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Matter at high density and low temperature is expected to be a color superconductor, which 
is a degenerate Fermi gas of quarks with a condensate of Cooper pairs near the Fermi surface 
that induces color Meissner effects. At the highest densities, where the QCD coupling is weak, 
rigorous calculations are possible, and the ground state is a particularly symmetric state, the 
color-flavor locked (CFL) phase. The CFL phase is a supcrfluid, an electromagnetic insulator, 
and breaks chiral symmetry. The effective theory of the low-energy excitations in the CFL phase 
is known and can be used, even at more moderate densities, to describe its physical properties. 
At lower densities the CFL phase may be disfavored by stresses that seek to separate the Fermi 
surfaces of the different flavors, and comparison with the competing alternative phases, which 
may break translation and/or rotation invariance, is done using phenomenological models. We 
review the calculations that underlie these results, and then discuss transport properties of several 
color-superconducting phases and their consequences for signatures of color superconductivity in 
neutron stars. 
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I. INTRODUCTION 
A. General outline 

The study of matter at ultra-high density is the "con- 
densed matter physics of quantum chromodynamics" . It 
builds on our understanding of the strong interaction, 
derived from experimental observation of few-body pro- 
cesses, to predict the behavior of macroscopic quanti- 
ties in many-body systems where the fundamental parti- 
cles of the standard model — quarks and leptons — become 
the relevant degrees of freedom. As in conventional 
condensed-matter physics, we seek to map the phase di- 
agram and calculate the properties of the phases. How- 
ever, we are in the unusual position of having a sector 
of the phase diagram where we can calculate many prop- 
erties of quark matter rigorously from first principles. 
This sector is the region of "asymptotically high" densi- 
ties, where quantum chromodynamics is weakly coupled. 
We will review those rigorous results and describe the 
progress that has been made in building on this solid 
foundation to extend our understanding to lower and 
more phenomenologically relevant densities. Quark mat- 
ter occurs in various forms, depending on the temper- 
ature T and quark chemical potential ^ (see Fig. [Ij. 
At high temperatures (T ^ ^) entropy precludes any 
pattern of order and there is only quark-gluon plasma 
(QGP), the phase of strongly interacting matter that has 
no spontaneous symmetry breaking, and which filled the 
universe for the first microseconds after the big bang. 
Quark-gluon plasma is also being created in small, very 
short-lived, droplets in ultrarelativistic heavy ion colli- 
sions at the Relativistic Heavy Ion Collider. 

In this review we concentrate on the regime of rela- 
tively low temperatures, T <^ fj,, where we find a rich 
variety of spontaneous symmetry breaking phases. To 
create such material in nature requires a piston that can 
compress matter to super-nuclear densities and hold it 
while it cools. The only known context where this might 
happen is in the interior of neutron stars, where gravity 
squeezes the star to an ultra-high density state where it 
remains for millions of years. This gives time for weak in- 
teractions to equilibrate, and for the temperature of the 
star to drop far below the quark chemical potential. We 
do not currently know whether quark matter exists in the 
cores of neutron stars. One of the reasons for studying 
color superconductivity is to improve our understanding 



of how a quark matter core would affect the observable 
behavior of a neutron star, and thereby resolve this un- 
certainty. 

When we speak of matter at the highest densities, we 
shall always take the high density limit with up, down 
and strange quarks only. We do so because neutron star 
cores are not dense enough (by more than an order of 
magnitude) to contain any charm or heavier quarks, and 
our ultimate goal is to gain insight into quark matter 
at densities that may be found in nature. For the same 
reason we focus on temperatures below about ten MeV, 
which are appropriate for neutron stars that are more 
than a few seconds old. 

As we will explain in some detail, at low temperatures 
and the highest densities we expect to find a degener- 
ate liquid of quarks, with Cooper pairing near the Fermi 
surface that spontaneously breaks the color gauge sym- 
metry ("color superconductivity"). Speculations about 
the existence of a quark matter phase at high density go 
back to the earliest day s of the quar k mod e l of hadrons 
(iBoccaletti et all Il966t ICarruthersl Il973l: lltohl . Il970t 



llvanenko and Kurdgelaidz^ 1965HPacinil 1966( ). and the 
possibility of quark Cooper pairing was noted even before 
there was a comprehensive theory of the stro ng interac- 
tion (jlvanenko and Kurdgelaid^ . ll969l[l970l) . After the 
development of quantum chromody namics (QCD), with 
its property of asymp totic freedom ( Gross and Wilczekl . 
Il973l : iPolitzeil . [ 19731 ). it became clear that a quark 
matter phase would exist at sufficiently high densit 
(iBavm and Chinl. ll976HChapline and Naueiibcr g. 197j 
19771 : ICollins and Perrvl.ll975l: iFreedman and McLerra 



19781 . Il977t iKislinger and Morlevl . Il976[) and the study 
of quark C ooper p a iring was pioneered by Barro is and 
Frautschi ( Barroi i 119771 . 119791 : iFrautschi I1978D . who 
first used the t erm "color superconductivity " , and by 
Bailin and Love (jBailin and Lovelll979lll984[) . who clas- 
sified many of the possible pairing patterns. Iw asaki 
and Iwado (jlwasakit . ,1995; .Iwasaki and Iwadol . Il995[ ) per- 
formed mean-field calculations of single-flavor pairing in 
a Nambu-Jona-Lasinio (NJL) model, but it was not until 
the predicti o n of l arge pairing gaps ( Alford et aLl . ll998t 



iRapD et~al. . 1998t) and th e color- flavor locked (CFL) 



phase ( Alford et all Il999lj ) that the phenomenology of 
color-superconducting quark matter became widely stud- 
ied. At present there are many reviews of t he sub- 
ject fr om various stages in it s dev el opment (lAlford 



2001 



1984 : LBuballa' '2005a'; Hong', '2 0011 IhsuI. 120001: iHuang 



Alford a.nd Raiagopall . 2006 ; Ba,ilin and Love , 



2002; Raiagopa,! and WilczeR I200C : 

a, 



20051: i Nardi ilh, ^..^ .^^^^ 

Reddvl.l200^: IrS] . l2004l : iRischkd . l2004l : ISchafeii . l2003b : 



ShovkovvT 2005 ) ■ and the reader may wish to consult 
them for alternative presentations with different em- 
phases. As these reviews make clear, the last decade 
has seen dramatic progress in our understanding of dense 
matter. We are now able to obtain, directly from QCD, 
rigorous and quantitative answers to the basic question: 
"What happens to matter if you squeeze it to arbitrar- 
ily high density?" . In Sec. IIVI we will show how QCD 
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becomes analytically tractable at arbitrarily high den- 
sity: the coupling is weak and the physics of confinement 
never arises, since long-wavelength magnetic interactions 
are cut off, both by Landau damping and by the Meiss- 
ner effect. As a result, matter at the highest densities 
is known to be in the CFL phase, whose properties (see 
Sec. Ull are understood rigorously from first principles. 
There is a well-developed effective field theory describing 
the low energy excitations of CFL matter (see Sec. W} . 
so that at any density at which the CFL phase occurs, 
even if this density is not high enough for a weak-coupling 
QCD calculation to be valid, many phenomena can nev- 
ertheless be described quantitatively in terms of a few 
parameters, via the effective field theory. 

It should be emphasized that QCD at arbitrarily high 
density is more fully understood than in any other 
context. High energy scattering, for example, can be 
treated by perturbative QCD, but making contact with 
observables brings in poorly understood nonperturba- 
tive physics via structure functions and/or fragmentation 
functions. Or, in quark-gluon plasma in the high tem- 
perature limit much of the physics is weakly-coupled but 
the lowest energy modes remain strongly coupled with 
nonperturbative physics arising in the nonabelian color- 
magnetic sector. We shall see that there are no analogous 
difficulties in the analysis of CFL matter at asymptotic 
densities. 

If the CFL phase persists all the way down to the tran- 
sition to nuclear matter then we have an exceptionally 
good theoretical understanding of the properties of quark 
matter in nature. However, less symmetrically paired 
phases of quark matter may well intervene in the inter- 
mediate density region between nuclear and CFL mat- 
ter (Sec. II. E[) . We enumerate some of the possibilities 
in Sec. IIIII In principle this region could also be un- 
derstood from first principles, using brute- force numer- 
ical methods (lattice QCD) to evaluate the QCD path 
integral, but unfortunately current lattice QCD algo- 
rithms are defeated by the fermion si gn problem in t he 
high-density low-temperature regime ISchmidtl . 120061) .^ 
This means we have to use models, or try to derive in- 
formation from astrophysical observations. In Sec. IVII 
we sketch an example of a (Nambu-Jona-Lasinio) model 
analysis within which one can compare some of the pos- 
sible intermediate-density phases suggested in Sec. II. El 
We finally discuss the observational approach, which in- 
volves elucidating the properties of the suggested phases 
of quark matter fSecs. IVi.CI and [VIip . and then finding 
astrophysical signatures by which their presence inside 



^ Condensation of Coope r pairs of quarks h a s been studied on the 
lattic e in 2-color Q C D jA Ues et aL|, [200^: I Fukushima and lidal. 



2007 



200C 



Hands ef aU 1200 6. 1999; iKogut et aLl l200ll . 11995 
20021: iNishida et aU i200}) , for high isospi n den- 



sity rather than baryon density llKogut and Sinclaiit l2002l : 
ISon and Stephanovl. [20011: ISplittorff ef al\ . |2001^ and in NJL- 
type models llHands and Waltersll2004ll 



neutron stars might be established or ruled out using as- 
tronomical observations fSec. IVIlT)) . 



B. Inevitability of color superconductivity 

At sufficiently high density and low temperature it 
is a good starting point to imagine that quarks form 
a degenerate Fermi liquid. Because QCD is asymp- 
totically free — the interaction becomes weaker as the 
momentum transferred grows — the quarks near the 
Fermi surface are almost free, with weak QCD interac- 
tions between them. (Small-angle quark-quark scatter- 
ing via a low-momentum gluon is no problem because 
it is cut off by Landau damping, which, together with 
Debye screening, keeps perturbation theory at high den- 
sity much better con t rolled t han at hig h temperature 
(jPisarski and Rischkd . Il999al: ISmH . Il999l ).) The quark- 
quark interaction is certainly attractive in some chan- 
nels, since we know that quarks bind together to form 
baryons. As we will now argue, these conditions are suf- 
ficient to guarantee color superconductivity at sufficiently 
high density. 

At zero temperature, the thermodynamic potential 
(which we will loosely refer to as the "free energy") is 

= i? — /iiV, where E is the total energy of the sys- 
tem, /i is the chemical potential, and N is the number 
of fermions. If there were no interactions then the en- 
ergy required to add a particle to the system would be 
the Fermi energy Ep = /x, so adding or subtracting par- 
ticles or holes near the Fermi surface would cost zero 
free energy. With a weak attractive interaction in any 
channel, if we add a pair of particles (or holes) with the 
quantum numbers of the attractive channel, the free en- 
ergy is lowered by the potential energy of their attrac- 
tion. Many such pairs will therefore be created in the 
modes near the Fermi surface, and these pairs, being 
bosonic, will form a condensate. The ground state will be 
a superposition of states with different numbers of pairs, 
breaking the fermion number symmetry. This argument, 
original ly developed by Bard een, Cooper, and Schrieffer 
(BCS) (jBardeen et all Il957[ ) is completely general, and 
can be applied to electrons in a metal, nucleons in nu- 
clear matter, '^He atoms, cold fermionic atoms in a trap, 
or quarks in quark matter. 

The application of the BCS mechanism to pairing in 
dense quark matter is in a sense more direct than in 
its original setting. The dominant interaction between 
electrons in a metal is the repulsive Coulomb interac- 
tion, and it is only because this interaction is screened 
that the attraction mediated by phonons comes into play. 
This means that the effective interactions that govern 
superconductivity in a metal depend on band structure 
and other complications and are very difficult to deter- 
mine accurately from first principles. In contrast, in 
QCD the "color Coulomb" interaction is attractive be- 
tween quarks whose color wave function is antisymmetric, 
meaning that superconductivity arises as a direct conse- 
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qucncc of the primary interaction in the theory. This has 
two important consequences. First, at asymptotic den- 
sities where the QCD interaction is weak we can derive 
the gap parameter and other properties of color super- 
conducting quark matter rigorously from the underlying 
microscopic theory. Second, at accessible densities where 
the QCD interaction is stronger the ratio of the gap pa- 
rameter to the Fermi energy will be much larger than 
in conventional BCS superconducting metals. Thus, su- 
perconductivity in QCD is more robust, both in the the- 
oretical sense and in the phenomenological sense, than 
superconductivity in metals. 

It has long been known that, in the absence 
of pairing, an unscreened static magn etic interac - 
tion results in a "non- Fermi- liquid " (IBavm et al. 



I990I: iBovanovskv and de Vega . l2001allb : Brown et al 



2000al: IChakravartv aZ.I. Il995t "ICerhold et al l 12004 : 
Holstein et all Il973l: IIdd et al\. |2"004. 2006: Manuel 



2000allbl: iNavak and Wilczekl . Il994l : IPolchinskil . Il994: 
Vanderhevden and OUitraultl 119971 )^ However, in QCD 
the magnetic interaction is screened at nonzero frequency 
(Landau damping) and this produces a particularly mild 
form of non-Fermi-liquid behavior, as we describe in 
Sec. IV. A. 21 In the absence of pairing but in the pres- 
ence of interactions, there are still quark quasiparticles 
and there is still a "Fermi surface", and the BCS argu- 
ment goes through. This argument is rigorous at high 
densities, where the QCD coupling g is small. The en- 
ergy scale below which non- Fermi liquid effects would be- 
come strong enough to modify the quasiparticle picture 
qualitatively is parametrically of order exp(— const /p^) 
whereas the BCS gap that results from pairing is para- 
metrically larger, of order exp(— const /g) as we shall see 
in Sec. IIVI This means that pairing occurs in a regime 
where the basic logic of the BCS argument remains valid. 

Since pairs of quarks cannot be color singlets, the 
Cooper pair condensate in quark matter will break the 
local color symmetry SU{3)c, hence the term "color su- 
perconductivity" . The quark pairs play the same role 
here as the Higgs particle does in the standard model: 
the color-superconducting phases can be thought of as 
Higgs phases of QCD. Here, the gauge bosons that ac- 
quire a mass through the process of spontaneous symme- 
try breaking are the gluons, giving rise to color Meissner 
effects. It is important to note that quarks, unlike elec- 
trons, have color and flavor as well as spin degrees of free- 
dom, so many different patterns of pairing are possible. 
This leads us to expect a panoply of different possible 
color superconducting phases. 

As we shall discuss in Sec. [TTl at the highest densities 
we can achieve an ab initio understanding of the proper- 
ties of dense matter, and we find that its preferred state 
is the CFL phase of three-flavor quark matter, which is 
unique in that all the quarks pair (all flavors, all colors, 
all spins, all momenta on the Fermi surfaces) and all the 
nonabelian gauge bosons are massive. The suppression 
of all of the infrared degrees of freedom of the types that 
typically indicate either instability toward further con- 



densation or strongly coupled phenomena ensures that, 
at sufhcicntly high density, the CFL ground state, whose 
only infrared degrees of freedom are Goldstone bosons 
and an abelian photon, is stable. In this regime, quantita- 
tive calculations of observable properties of CFL matter 
can be done from first principles; there are no remaining 
nonperturbative gaps in our understanding. 

As the density decreases, the effect of the strange 
quark mass becomes more noticeable, imposing stresses 
that may modify the Cooper pairing and the CFL phase 
may be replaced by other forms of color superconduct- 
ing quark matter. Furthermore, as the attractive inter- 
action between quarks becomes stronger at lower densi- 
ties, correlations beyond the two-body correlation that 
yields Cooper pairing may become important, and at 
some point the ground state will no longer be a Cooper- 
paired state of quark matter, but something quite differ- 
ent. Indeed, by the time we decrease the density down to 
that of nuclear matter, the average separation between 
quarks has increased to the point that the interactions 
arc strong enough to bind quarks into nucleons. It is 
worth noting that quark matter is in this regard differ- 
ent from Cooper-paired ultracold fermionic atoms (to be 
discussed in Sec. lIII.l| . For fermionic atoms, as the inter- 
action strength increases there is a crossover from BCS- 
paired fermions to a Bose-Einstein condensate (BEC) 
of tightly-bound, well-separated, weakly-interacting di- 
atoms (molecules). In QCD, however, the color charge 
of a diquark is the same as that of an antiquark, so di- 
quarks will interact with each other as strongly as quarks, 
and there will not be a literal analogue of the BCS/BEC 
crossover seen in fermionic atoms. In QCD, the neutral 
bound states at low density that are (by QCD standards) 
weakly interacting are nucleons, containing three quarks 
not two. 

We shall work with Nc — 3 colors throughout. In 
the limit Nc ^ 00 with fixed Aqcd (i-e fixed g'^Nc), 
Cooper pairing is not necessarily energetically preferred. 
A strong competitor for the large- TVc ground state is the 
chiral density wave (CDW), a condensa te of quark-hole 
pairs , each with total momentum 2pp ( Dervagin et all 
Il992f ). Quark-hole scattering is enhanced by a factor of 
Nc over quark-quark scattering, but, unlike Cooper pair- 
ing, it only uses a small fraction of the Fermi surface, and 
in the case of short range forces the CDW phase is ener- 
getically favored in one-dimensio nal systems, bu t not in 
two or more spatial dimensions ( Shankarl . 1 1994 ). How- 
ever, in QCD in the large Nc limit the equations govern- 
ing the CDW state become effectively one-dimensional 
because the gluon propagator is not modified by the 
medium, so the quark-hole interaction is dominated by 
almost coUinear scattering. Since pairing gaps are expo- 
nentially small in the coupling but medium effects only 
vanish as a power of the CDW state requires an 
exponentially large number of colors. It is estimated 
that for /i ~ 1 GeV, quark-hole pairin g becomes favored 
over Cooper pairing when Nc > 1000 (Shuster and Soiil . 
|2000( ). Recent work (jMcLerran and Pisarskil l2007t) dis- 
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cusses aspects of physics at large Nc at lower densities 
that may also be quite different from physics at Nc = 3. 

Before turning to a description of CFL pairing in Scc.HIl 
and Icss-symmetrically paired forms of color supercon- 
ducting quark matter in Sec. lIIH we discuss some generic 
topics that arise in the analysis of color-superconducting 
phases: the gap equations, neutrality constraints, the 
resultant stresses on Cooper pairing, and the expected 
overall form of the phase diagram. 



C. Quark Cooper pairing 



the condensate: 

We have written this gap equation and stability condi- 
tion somewhat schematically since for many patterns of 
pairing there will be gap parameters with different mag- 
nitudes in different channels. The free energy must then 
be minimized with respect to each of the gap parameters, 
yielding a coupled set of gap equations. The solution to 
with the lowest free energy that respects the neutral- 
ity constraints discussed below yields the favored phase. 



The quark pair condensate can be characterized in a 
gauge-variant way by the expectation value of the one- 
particle-irreducible quark-quark two-point function, also 
known as the "anomalous self-energy" , 



JO' ij ao 



(1) 



Here ip is the quark field operator, color indices a, (3 
range over red, green, and blue {r,g,b), flavor indices 
I, j range over up, down and strange (m, d, s), and a, b are 
the spinor Dirac indices. The angle brackets denote the 
one-particle-irreducible part of the quantum-mechanical 
ground-state expectation value. In general, both sides 
of this equation arc functions of momentum. The color- 
flavor-spin matrix P characterizes a particular pairing 
channel, and A is the gap parameter which gives the 
strength of the pairing in this channel. A standard BCS 
condensate is position-independent (so that in momen- 
tum space the pairing is between quarks with equal and 
opposite momentum) and a spin singlet (so that the gap 
is isotropic in momentum space). However, as we will see 
later, there is good reason to expect non-BCS conden- 
sates as well as BCS condensates in high-density quark 
matter. 

Although ID) defines a gauge-variant quantity, it is 
still of physical relevance. Just as electroweak symme- 
try breaking is most straightforwardly understood in the 
unitary gauge where the Higgs vacuum expectation value 
is uniform in space, so color superconductivity is typically 
analyzed in the unitary gauge where the quark pair op- 
erator has a uniform color orientation in space. We then 
relate the gap parameter A to the spectrum of the quark- 
like excitations above the ground state ( "quasiquarks" ) , 
which is gauge-invariant. 

In principle, a full analysis of the phase structure of 
quark matter in the fi-T plane would be performed by 
writing down the free energy fl, which is a function 
of the temperature, the chemical potentials for all con- 
served quantities, and the gap parameters for all possible 
condensates, including the quark pair condensates but 
also others such as chiral condensates of the form {qq). 
We impose neutrality with respect to gauge charges (see 
Sect. II. D] below) and then within the neutral subspace we 
minimize the free energy with respect to the strength of 



D. Chemical potentials and neutrality constraints 

Why do we describe "matter at high density" by intro- 
ducing a large chemical potential ^ for quark number but 
no chemical potentials for other quantities? The answer 
is that this reflects the physics of neutron stars, which 
arc the main physical arena that we consider. Firstly, 
on the long timescales relevant to neutron stars, the 
only global charges that are conserved in the standard 
model are quark number and lepton number, so only 
these can be coupled to chemical potentials (we shall 
discuss gauged charges below). Secondly, a neutron star 
is permeable to lepton number because neutrinos are so 
light and weakly-interacting that they can quickly escape 
from the star, so the chemical potential for lepton num- 
ber is zero. Electrons are present because they carry 
electric charge, for which there is a nonzero potential. 
In the first few seconds of the life of a neutron star the 
neutrino mean free path may be short enough to sus- 
tain a no nzero lepton numbe r chem i cal potential, see for 
instance (iBerdermann et all . 2004 Kaplan and Reddyl . 

" l2006^ ■ 



I2OO2I : iLaporta and Ruggieril . l2006l:TRuester al\ . 12006 
but we will not discuss that scenario. 

Stable bulk matter must be neutral under all gauged 
charges, whether they are spontaneously broken or 
not. Otherwise, the net charge density would create 
large electric fields, making the energy non-extensive. 
In the case of the electromagnetic gauge symmetry, 
this simply requires zero charge density, Q = 0. The 
correct formal requirement concerning the color charge 
of a large lump of matter is that it should be a color 
singlet, i.e., its wavefunction should be invariant under 
a general color gauge transformation. However, it is 
sufficient for us to impose color neutrality, meaning 
equality in the numbers of red, green, and blue quarks. 
This is a less stringent constraint (singlet neutral but 
neutral 7^ singlet) but the projection of a color neutral 
state onto a color singlet c osts no extra free e nergy in 
the t hermodynamic li r nit (|Amore et all |2002[ ). (See 
also (|Elze etall Il983l . [198^.) In general there are 8 
possible color charges, but because the Cartan subalge- 
bra of SU (3)c is two-dimensional it is always possible to 
transform to a gauge where all are zero except Q3 and 
Qs, the charges associated with the diagonal generators 
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diag(l,l,-2) 
(iBuballa and Shovkovvl l2005t 



T3 = idiag(l,-l,0) and Tf 
in (r,q,b) space 

iRaiagopal and Schmittl . 20061 ). In this review, we 
only discuss such gauges. So to impose color neutrality 
we just require Q3 = Qs ~ 0. 

In nature, electric and color neutrality arc enforced 
by the dynamics of the electromagnetic and QCD gauge 
fields, whose zeroth components serve as chemical po- 
tentials coupled to the charges Q,Q3,Qs, E^nd which 
are n aturally driven to values that set these charges to 
zero ( Alford and Raiagopall 20021 : jPietrich and Rischke . 
2004 iGerhold and RebhaJ 120031; [Rda and Bavml. 12001: 



Krvievskil . l2003f ) . In an NJL model with fermions but no 



gauge fields (sec Sec. IVI[) one has to introduce the chem- 
ical potentials /Ze, /is and fig, by hand in order to enforce 
color and electric neutrality. The neutrality conditions 
are then 



Q= ^ ^ 



dfle 

on 

dVL 



(3) 



(Note that we define an electrostatic potential fie that 
is coupled to the negative electric charge Q, so that in 
typical neutron star conditions, where there is a finite 
density of electrons rather than positrons, fie is positive.) 

Finally we should note that enforcing local neutral- 
ity is appropriate for uniform phases, but there are also 
non-uniform charge-separated phases ("mixed phases"), 
consisting of positively and negatively charged domains 
which are neutral on average. These are discussed further 
in Sec. UlLHl 



E. Stresses on BCS pairing 

The free energy argument that we gave in Sec. II.BI for 
the inevitability of BCS pairing in the presence of an 
attractive interaction relies on the assumption that the 
quarks that pair with equal and opposite momenta can 
each be arbitrarily close to their common Fermi surface. 
However, as wc will sec in Sec. |lTl the neutrality con- 
straint, combined with the mass of the strange quark and 
the requirement that matter be in beta equilibrium, tends 
to pull apart the Fermi momenta of the different flavors 
of quarks, imposing an extra energy cost ("stress") on 
the formation of Cooper pairs involving quarks of differ- 
ent flavors. This raises the possibility of non-BCS pairing 
in some regions of the phase diagram. 

To set the stage here, let us discuss a simplified exam- 
ple: consider two massless species of fermions, labeled 1 
and 2, with different chemical potentials fii and fi2, and 
an attractive interaction between them that favors cross- 
species BCS pairing with a gap parameter A. It will turn 
out that to a good approximation the color-flavor locked 



pairing pattern contains three such sectors, so this ex- 
ample captures the essential physics we will encounter in 
later sections. We define the average chemical potential 
and the stress parameter 



5fl = \{fli - fl2) 



(4) 



As long as the stress 5fi is small enough relative to 
A, BCS pairing between species 1 and 2 can occur, 
locking their Fermi surfaces together and ensuring that 
they occur in e qual numbers. At th e Chandrasekhar - 
Clogston point (jChandrasekhail Il962l : IClogstonl . Il962h . 
where 5fi = A/-\/2, the two-species model undergoes a 
first-order transition to the unpaired phase. At this point 
BCS pairing still exists as a locally stable state, with 
a completely gapped spectrum of quasiparticles. When 
5fi reaches A the spectrum becomes gapless at momen- 
tum p = fi, indicating that cross-s pecies BCS p airing is 
no longer favored at all momenta (| Alford et aZ.I . [20 04b). 
If the two species are part of a larger pairing pattern, 
the Chandrasekhar-Clogston transition can be shifted, 
and we shall see that in the two-species subsectors of 
the CFL pattern it is shifted to Sfi > A. The onset 
of gaplessness is therefore the relevant threshold for our 
purposes, and it always occurs at Sfi = A, independent 
of the larger context in which the two flavors pair. This 
follows from the fact that BCS pairing only occurs if the 
energy gained from turning a 1 quark into a 2 quark with 
the same momentum (namely fii — fi2) is less than the 
cost of breaki ng the Cooper pair formed by t hese quarks, 
which is 2 A (jUaiagopal and WilczeQ . 1200 ll ). Thus the 
1-2 Cooper pairs are energetically stable (or metastable) 
as long as < A. A more detailed treatment of this 
illust rative example can be found in ( Alford and Wari3 . 
|2005[) . 

This example uses massless quarks, but it can easily 
be modified to include the leading effect of a quark mass 
M . A difference in the masses of the pairing quarks also 
stresses the pairing, because it gives them different Fermi 
momenta at the same chemical potential, so the quarks 
in a 1-2 Cooper pair, which have equal and opposite mo- 
menta, will not both be close to their Fermi energies. The 
leading-order effect is easily calculated, since for a quark 
near its Fermi surface it acts like a shift in the quark 
chemical potential by — M^/(2/2) (given that Fermi mo- 
mentum pp p, to this order). 

Returning from our toy model to realistic quark mat- 
ter, the quark fiavors that are potentially relevant at 
neutron-star densities are the light up (u) and down 
(d) quarks, with current masses m„ and to^ that are 
< 5 MeV, and a medium-weight flavor, the strange (s) 
quark, with current mass ~ 90 MeV. Their effec- 
tive "constituent" masses in the vacuum are hundreds of 
MeV larger, but are expected to decrease with increasing 
quark density. We shall refer to the density-dependent 
constituent masses as Mu.d,s and shall typically neglect 
Mu and Aid- As our toy model has illustrated, however, 
the strange quark mass Mg will contribute to stresses on 
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discuss the factors that disfavor the CFL phase at inter- 
mediate densities, and survey the color superconducting 
phases that have been hypothesized to occur there. 

Various aspects of color superconductivity at high tem- 
peratures have been studied, including the phase struc- 
ture (see Sec. IVI.A[) . spectral functions, pair-forming 
and -breaking fluctuations, possible precursors to con- 
densation su ch as pseudo g aps, a nd various collective 
phenomena llAbuki et aL. 
20051 : iHatsuda et all \200^ 



■20021: iFukushima and lidal , 
Kitazawa et all l2002l 12004 . 



2005aibl . l2007l : IVoskresenskvl |2004 lYamamoto et aL . 



20071 ). However, this review centers on quark matter at 



neutron star temperatures, and throughout Sees. HTl and 
mil we restrict ourselves to the phases of quark matter 
at zero temperature. This is because most of the phases 
that we discuss are expected to persist up to critical tem- 
peratures that are well above the core temperature of a 
typical neutron star, which drops below 1 MeV within 
seconds of its birth before cooling down through the keV 
range over millions of years. 



FIG. 1 (Color online) A schematic outline for the phase dia- 
gram of matter at ultra-high density and temperature. The 
CFL phase is a superfluid (like cold nuclear matter) and has 
broken chiral symmetry (like the hadronic phase). 



II. MATTER AT THE HIGHEST DENSITIES 



A. Color-flavor locked (CFL) quark matter 



cross-flavor pairing, and those stresses will become more 
severe as the density (and hence ji) decreases. This will 
be a major theme of later sections. 



F. Overview of the quark matter phase diagram 

Fig. [1] shows a schematic phase diagram for QCD that 
is consistent with what is currently known. Along the 
horizontal axis the temperature is zero, and the density 
is zero up to the onset transition where it jumps to nu- 
clear density, and then rises with increasing ^. Neutron 
stars are in this region of the phase diagram, although 
it is not known whether their cores are dense enough to 
reach the quark matter phase. Along the vertical axis the 
temperature rises, taking us through the crossover from 
a hadronic gas to the quark-gluon plasma. This is the 
regime explored by high-energy heavy-ion colliders. 

At the highest densities we find the color-flavor locked 
color-superconducting phase, ^ in which the strange quark 
participates symmetrically with the up and down quarks 
in Cooper pairing. This is described in more detail in 
Sees, ini IIVI and |Vl It is not yet clear what happens 
at intermediate density, and in Sees. |TTT] and |Vl] we will 



^ As explained in Sec. lI.Al we fix Nf = 3 at all densities, to main- 
tain relevance to ne utron star inter iors. Pairing with arbitrary 
Nf has been studied llSchaferl . l2000al ). For Nf a multiple of three 
one finds multiple copies of the CFL pattern; for Nf = 4, 5 the 
pattern is more complicated. 



Given that quarks form Cooper pairs, the next ques- 
tion is who pairs with whom? In quark matter at suf- 
ficiently high densities, where the up, down and strange 
quarks can be treated on an equal footing and the disrup- 
tive effects of the strange quark mass can be neglected, 
the most symmetric and most attractive option is the 
color-flavor locked phase, where quarks of all three colors 
and all three flavors form conventional zero-momentum 
spinless Cooper pairs. This pattern, anticipated in early 
studies of alternative condensa tes for zero-dens i ty ch i- 
ral symmetry breaking (Sredni cki and Susskind . 198ll) . 
is enco ded in the quark-quark self-energy ( Alford et all . 
ll999bD 



The symmetry breaking pattern is 



[SU{3)c] X U{1)b 
X SUi3)L X SU{3)b. ^ SU{3) c+L+R 



X Z2 



(5) 



(6) 



Color indices a,/? and flavor indices i,j run from 1 to 3, 
Dirac indices are suppressed, and C is the Dirac charge- 
conjugation matrix. Gauge symmetries are in square 
brackets. Acfl is the CFL gap parameter. The Dirac 
structure C75 is a Lorentz singlet, and corresponds to 
parity-even spin-singlet pairing, so it is antisymmetric in 
the Dirac indices. The two quarks in the Cooper pair are 
identical fermions, so the remaining color-|-flavor struc- 
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turc must be symmetric. The dominant color-flavor com- 
ponent in (O transforms as (3^,3^1), antisymmetric in 
both. The subdominant term, muhiphed by k, trans- 
forms as (6s, 6s). It is almost certainly not energetically 
favored on its own (all the arguments in Sec. III. A. 51 for 
the color triplet imply repulsion for the sextet), but in 
the presence of the dominant pairing it breaks no ad- 
ditio nal symmetries, so k is in general small but not 
zero (lAlford et all Il999bt iPisarski and Rischkl. Il999ct 



ISchafeil . l2000al: IShovkovv and Wiiewardhanal . Il999t) . 



1. Color-flavor locking and chiral symmetry breaking 

A particularly striking feature of the CFL pairing pat- 
tern is that it breaks chiral symmetry. Because of color- 
flavor locking, chiral symmetry remains broken up to ar- 
bitrarily high densities in three-flavor quark matter. The 
mechanism is quite different from the formation of the 
(i/'V') condensate that breaks chiral symmetry in the vac- 
uum by pairing left-handed (L) quarks with right-handed 
(R) antiquarks. The CFL condensate pairs L quarks with 
each other and R quarks with each other (quarks in a 
Cooper pair have opposite momentum, and zero net spin, 
hence the same ehirality) and so it might naively appear 
chirally symmetric. However, the Kronecker deltas in 
([5|) connect color indices with flavor indices, so that the 
condensate is not invariant under color rotations, nor un- 
der flavor rotations, but only under simultaneous, equal 
and opposite, color and flavor rotations. Color is a vec- 
tor symmetry, so the compensating flavor rotation must 
be the same for L and R quarks, so the axial part of 
the flavor group, which is the chiral symmetry, is broken 
by the locking of colo r and flavor rotations to each other 
(|Alford et aLi . il999bf ). Such locking is familiar from other 
contexts, including the QCD vacuum, where a conden- 
sate of quark-antiquark pairs locks S'C/(3)l to SU{3)r 
breaking chiral symmetry "directly" , and the B phase 
of superfluid ■^He, where the condensate transforms non- 
trivially under rotations of spin and orbital angular mo- 
mentum, but is invariant under simultaneous rotations 
of both. 

The breaking of the chiral symmetry is associated with 
an expectation value for a gauge-invariant order param- 
eter with the structure i/j'ipipip (see Sec.|V|. There is also 
a subdomi nant "conventio nal" chiral condensate {'ipip) ^ 
(V'C75V') (ISchafed . l2000aD . These gauge-invariant ob- 
servables distinguish the CFL phase from the QGP, and 
if a lattice QCD algorithm applicable at high density ever 
becomes available, they could be used to map the pres- 
ence of color-flavor locking in the phase diagram. 

We also expect massless Goldstone modes associated 
with chiral symmetry breaking (see Sees. in.A.4l and IV)) . 
In the real world there is small explicit breaking of chi- 
ral symmetry from the current quark masses, so the order 
parameters will not go to zero in the QGP, and the Gold- 
stone bosons will be light but not massless. 



2. Superfluidity 

The CFL pairing pattern spontaneously breaks the 
exact global baryon number symmetry U{1)b, leaving 
only a discrete Z2 symmetry under which all quark fields 
are multiplied by —1. There is an associated gauge- 
invariant 6-quark order parameter with the flavor and 
color structure of two Lambda baryons, (AA) where 
^ _ gobcg^^^a^&^c^ This order parameter distinguishes 
the CFL phase from the QGP, and there is an associ- 
ated massless Goldstone boson that makes the CFL phase 
a superfluid, see Sec. IV.C.2I The vortices that result 
when CFL quark matter is rotated have been studied in 
tealaehandran et al\ '2006"; 'Forbes and Zhitmtskyl, 120021 : 
llida and Bavm . ,2002^ ,Nakano et aL . ,2007l) . 



3. Gauge symmetry breaking and electromagnetism 

As explained above, the CFL condensate breaks the 
SU{3)c X SU{3)l X SU{3)r symmetry down to the diag- 
onal group SU{3)c+L+R of simultaneous color and flavor 
rotations. Color is a gauge symmetry, and one of the gen- 
erators of SU{3)l+b. is the electric charge, which gener- 
ates the U{1)q gauge symmetry. This means that the 
unbroken SU{3)c+l+fi. contains one gauged generator, 
corresponding to an unbroken U{1)q which consists of 
a simultaneous electromagnetic and color rotation. The 
rest of the color group is broken, so by the Higgs mech- 
anism seven gluons and one gluon-photon linear com- 
bination become massive via the Meissner effect. The 
orthogonal gluon-photon generator Q remains unbroken, 
because every diquark in the condensate has Q = 0. The 
mixing angle is cos 9 = g/ \/ g'^ + 4e^/3 where e and g are 
the QED and QCD couplings. Because e g the angle 
is close to zero, meaning that the Q photon is mostly the 
original photon with a small admixture of gluon. 

The Q photon is massless. Given small but nonzero 
quark masses, there are no gapless Q-charged excitations; 
the lightest ones are the pseudoscalar pseudo-Goldstone 
bosons TT* and (see Sees. III. A. 41 and |V | . so for tem- 
peratures well below their m asses (and well below the 
electron mass (jShovkovv and E llis. 2003il) the CFL phase 
is a transparent insulator, in which Q-electrie and mag- 
netic fields satisfy Maxwell's equations with a dielectric 
constant and index of refraction that can be calculated 
directly from QCD (|Litim and Manuell . [200l[) . 



n = l 



2^ eos^ 9 



97r2 



^CFL 



(7) 



(This result is valid as long as n — 1 <§; 1.) Apart from 
the fact that n 7^ 1, the emergence of the Q photon is 
an exact QCD-scale analogue of the TeV-scale sponta- 
neous symmetry breaking that gave rise to the photon as 
a linear combination of the and hypercharge gauge 
bosons, with the diquark condensate at the QCD scale 
playing the role of the Higgs condensate at the TeV scale. 
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If one could shine a beam of ordinary light on a 
lump of CFL matter in vacuum, some would be re- 
flected and some would enter, refracted, as a beam 
of Q-light. The reflection and ref r action coefficients 
are known (iManuel and Raiagopal |2002| ) (see also 
(| Alford and Goodl . 120041 )'). The static limit of this aca- 
demic result is relevant: if a volume of CFL matter finds 
itself in a static magnetic field as within a neutron star, 
surface currents are induced such that a fraction of this 
field is expelled via the Meissner effect for the non-Q com- 
pone nt of Q, while a fracti on is admitted as Q-magnetic 
field (|Alford et all . l2000bl ). The magnetic field within 
the CFL volume is not confined to fiux tubes, and is not 
frozen as in a conducting plasma: CFL quark matter is 
a color superconductor but it is an electromagnetic insu- 
lator. 

All Cooper pairs have zero net Q-charge, but some 
have neutral constituents (both quarks Q-neutral) and 
some have charged constituents (the two quarks have 
opposite Q-charge). The Q-componcnt of an external 
magnetic field will not affect the first type, but it will 
affect the pairing of the second type, so external mag- 
netic fields can modify the CFL phase to the so-called 
magnetic CFL ("MCFL") p hase. The MCFL phase 
has a different gap structure ( Ferrer et al. . 20051 200 
and a different effective theory (jFerrer and de la Incer 
l2007bD . The original analyses of the MCFL phase were 
done for rotated magnetic fields B large enough that 
all quarks are in the lowest Landau level; solving the 
gap equations at lower B shows that the gap parameters 
in the MCFL phase ex hibit de Haas-van Alphen oscilla- 
tions, periodic in 1/B ( Fukushima and Warringal l2007t 



iNoronha and Shovkovvl l2007t) 



4. Low-energy excitations 

The low-energy excitations in the CFL phase are: the 8 
light pseudoscalars arising from broken chiral symmetry, 
the massless Goldstone boson associated with superfiuid- 
ity, and the Q-photon. The pseudoscalars form an octet 
under the unbroken SU{S) color-|-flavor symmetry, and 
can naturally be labeled according to their Q-charges as 
pious, kaons, and an 77. The effective Lagrangian that 
describes their interactions, and the QCD calculation of 
their masses and decay constants will be discussed in 
Sec. |Vl "We shall find, in particular, that even though 
the quark-antiquark condensate is small, the pion decay 
constant is large, /,r ~ 

The symmetry breaking pattern ^ does not include 
the spontaneous breaking of the U{1)a "symmetry" be- 
cause it is explicitly broken by instanton effects. How- 
ever, at large densities these effects become arbitrar- 
ily small, and the spontaneous breaking of U{1)a will 
have an associated order parameter and a ninth pseudo- 
Goldstone boson with the quantum numbers of the 77'. 
This introduces the possibility of a second type of vor- 
tices ( Forbes and Zhitnitskvl . [2OO2I: ISon et aZ.I . l200l[) . 



Among the gapped excitations, we find the quark- 
quasiparticles which fall into an 8 ® 1 of the unbro- 
ken global SU(3)c+L+R,i so there are two gap parame- 
ters Al and Ag. The singlet has the larger gap Ai = 
(2 + 0{k))As. We also find an octet of massive vector 
mesons, which are the gluons that have acquired mass 
via the Higgs mechanism. The symmetries of the 3- 
fiavor CFL phase are the same as those one would ex- 
pect for 3-fiavor hypernuclear matter, and even the pat- 
tern of gapped excitations is remarkably similar, differ- 
ing only in the absence of a ninth massive vector me- 
son. It is therefore possible that there is no phase transi- 
tion between hypernu c lear m atter and CFL quark matter 
( Schafer and Wilczekl . [19993) . This hadron-quark conti- 
nuity can arise in nature only if the strange quark is 
so light that there is a hypernuclear phase, and this 
phase is characterized by proton-S~, neutron-S*^ and 
E^-E" pairing, which can then continuously evolve into 
CFL quark matter u pon further increasing the density 
(lAlford aLl . ll999aD . 



5. Why CFL is favored 

The dominant component of the CFL pairing pattern is 
the color 3^1, flavor 3^, and Dirac C75 (Lorentz scalar). 
There are many reasons to expect the color 3a to be 
favored. First, this is the most attractive channel for 
quarks interacting via single-gluon exchange which is the 
dominant interaction at high densities where the QCD 
coupling is weak; second, it is also the most attractive 
channel for quarks interacting via the instanton-induced 
't Hooft interaction, which is important at lower den- 
sities; third, qualitatively, combining two quarks that 
are each separately in the color-3 representation to ob- 
tain a diquark that is a color-3yi lowers the color- 
flux at large distances; and, fourth, phenomenologically, 
the idea that baryons can be modeled as bound states 
of a quark and a color-antisymmetric diquark, taking 
advantage of the attraction in this diquark channel, 
has a long histo ry and h as had a recent ren ai ssance 
(lAnselmino et all. [l993t IClose and Tornqvistl. 12002 : 



Jaffe l . 119771: 1 Jaffe and Wilczekl . l2003HSelem and Wilraek 
20061) . 

It is also easy to understand why pairing in the 
Lorentz-scalar channel is favorable: it leaves rotational 
invariance unbroken, allowing for quarks at all angles 
on the entire Fermi-sphere to participate coherently in 
the pairing. Many calculations have shown that pair- 
ing is weaker in ch a nnels that b reak rotationa l sym- 
metrv (lAlford et all. l2003l. Il998l ; [ Buballa et dl. \200c 
Iwasaki and IwadoT ll995l: ISchafe'R [2000b: Sc hmitt et al 
2OO2D . There is also a rotationally invariant pairing 



channel with negative parity described by the order 
parameter {ipCip). Perturbative gluon exchange inter- 
actions do not distinguish between positive and nega- 
tive parity diquarks, but non-perturbativc instanton in- 
duced interactions do, favoring the positive parity chan- 
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nel (|Alford et all Il998l : iRapp et al l. ll998ll2000D . 

Once we have antisymmetry in color and in Dirac in- 
dices, we are forced to antisymmetrize in flavor indices, 
and the most general color-flavor structure that the ar- 
guments above imply should be energetically favored is 



Unpaired 



2SC pairing 



CFL pairing 



A 

B ■ 



(8) 



CFL pairing corresponds to (j)-^ = and this is the 
only pattern that pairs all the quarks and leaves an en- 
tire SU (3) global symmetry unbroken. The 2SC pattern 
is (jijj = 6^ 6 n, in which only m and d quarks of two 
colors pair (lAlford et all 1998 : iBailin and Lovd , Il984t 
lBarroisl . [l9"79HRapp et an . ll998h . see SecIlLAl As long 
as the strange quark mass can be neglected (the para- 
metric criterion turns out to be Acfl ^ -^s/M; see 
Sec. nil. Bp calculations comparing patterns of the struc- 
ture ([5]) always find the CFL phase to have the highest 
condensation energy, making it the favored pattern. This 
has been confirmed i n weak-cou pling Q CD calculations 

valid at high density jEvans et al. 200Q : ISchafe^, |2000at 

IShovkovv and Wijewardhana . 1999), in the G inzburg- 
Landau approximation ( lida and Bavml . l200ll ). and in 
many calcula t ions u si ng Nambu-J o na-Lasinio models 
(lAlford et all Il999allbl: iMalekzadehl . l2006l : iRapp et all 
I2OOOI : ISchafer and Wilczekl . Il999ct) . In the high-density 
limit where A ^ M^/fi and A <|; /i we can expand in 
powers of A//i and explicitly compare CFL to 2SC pair- 
ing. The CFL condensation energy is (8A|-|-Aj)^^/(47r^) 
which is 12AcFLA*^/(47r^) when k < 1 (see Sec. IILA.4|) 
whereas the condensation energy in the 2SC phase is only 
4A2scA^^/(47r^). We shaU see later that the 2SC gap pa- 
rameter turns out to be larger than the CFL gap param- 
eter by a factor of 2^^^, so up to corrections of order n 
the CFL condensation energy is larger than that in the 
2SC phase by a factor of 3 x 2~^/^. At lower densities 
the condensation energies become smaller, and we can- 
not neglect negative Mf terms which are energy penalties 
induced by the neutrality requirement. Their coefRcient 
is larger for CFL than for 2SC, partly (but usually not 
completely) cancelling the extra condensation energy — 
see Fig. 13] and Sec.|llIS 



B. Intermediate density: stresses on the CFL phase 



As we noted in section lI.El BCS pairing between 
two species is suppressed if their chemical potentials are 
sufficiently different. In real- world quark matter such 
stresses arise from the strange quark mass, which gives 
the strange quark a lower Fermi momentum than the 
down quark at the same chemical potentials fi and fi^, 
and from the neutrality requirement, which gives the up 
quark a different chemical potential from the down and 
strange quarks at the same /x and He- Once flavor equilib- 
rium under the weak interactions is reached, we find that 
all three flavors prefer to have different Fermi momenta 
at the same chemical potentials. This is illustrated in 



red green blue 



red green blue 



red green blue 



FIG. 2 (Color online) Illustration of the splitting apart of the 
Fermi momenta of the various colors and flavors of quarks 
(exaggerated for easy visibility). In the unpaired phase, re- 
quirements of neutrality and weak interaction equilibration 
cause separation of the Fermi momenta of the various flavors. 
The splittings increase with decreasing density, as fi decreases 
and Ms{fi) increases. In the 2SC phase, up and down quarks 
of two colors pair, locking their Fermi momenta together. In 
the CFL phase, all colors and flavors pair and have a common 
Fermi momentum. 



Fig. [21 which shows the Fermi momenta of the different 
species of quarks. 

In the unpaired phase (Fig. [2l left panel) , the strange 
quarks have a lower Fermi momentum because they are 
heavier, and to maintain electrical neutrality the number 
of down quarks is correspondingly increased. To lowest 
order in the strange quark mass, the separation between 
the Fermi momenta is Spp = M'^ / (4/j,) , so the splitting 
becomes larger as the density is reduced, and smaller as 
the density is increased. The phase space at the Fermi 
surface is proportional to fi^, so the resultant difference 
in quark number densities is Ud—riu = n^—ns (x fi^Spp ^ 
fJ^M^. Electrons are also present in weak equilibrium, 
with fie = M's/(4:iJ.), so their charge density is paramet- 
rically of order fj.^ ~ / fi^ <C fJ^M'^, meaning that they 
are unimportant in maintaining neutrality. 

In the CFL phase all the colors and flavors pair with 
each other, locking all their Fermi momenta together at 
a common value (Fig. [21 right panel). This is possi- 
ble as long as the energy cost of forcing all species to 
have the same Fermi momentum is compensated by the 
pairing energy that is released by the formation of the 
Cooper pairs. Still working to lowest order in , we 
can say that parametrically the cost is fJ,^Sp^ ~ Mf, 
and the pairing energy is ^^A^p^, so we expect CFL 
pairing to become disfavored when Acfl ^ -^s/m- I^i 
fact, the CFL phase remains fav ored over the unpaired 
phase as long as Acfl > / 4/i (jAlford and Raiagopall 
|2002| ). but already becomes unstable against unpair- 



model calculations (lAbuki et al. 


. 20051 Alford et al. 


2005cl: Alford and Raiagopall 2002 


; Blaschke et al.. 2005; 


Fukushima et all 2005; iRiister et al.. 20051) find that if 



the attractive interaction were strong enough to induce 
a 100 MeV CFL gap when = then the CFL phase 
would survive all the way down to the transition to nu- 



11 



clear matter. Otherwise, there must be a transition 
to some other quark matter phase: this is the "non- 
CFL" region shown schematicahy in Fig. [TJ When the 
stress is smaU, the CFL pairing can bend rather than 
break, developing a condensate of K° mesons, described 
in Sec. III. CI below. When the stress is larger, however, 
CFL pairing becomes disfavored. A comprehensive sur- 
vey of possible BCS pairing patterns shows that all of 
them suffer from the stress o f Fermi surface splitting 
( Raiagopal and Schmittl . l2006l ). so in the intermediate- 
density "non-CFL" region we expect more exotic non- 
BCS pairing patterns. In Sec. IIIII we give a survey of 
possibilities that have been explored. 



C. Kaon condensation: the CFL-7<'° phase 



Bedaque and Schafer ( Bedaque and Schafed . l2002t ) 
showed that when the stress is not too large (high den- 
sity), it may simply modify the CFL pairing pattern 
by inducing a flavor rotation of the condensate. This 
modification can be interpreted as a condensate of "JC^" 
mesons. The K'^ meson carries negative strangeness (it 
has the same strangeness as a s quark), so forming a K'^ 
condensate relieves the stress on the CFL phase by re- 
ducing its strangeness content. At large density kaon 
condensation occurs for Ms > m^^'^A^/'^, where m is 
mass of the light (u and d) quarks. At moderate density 
the critical strange quark mass i s increased by in stanton 
contribution to the kaon mass (ISchafed l2002aD . Kaon 
condensation was initially demonstrated using an effec- 
tive theory of the Goldstone bosons, but with some effort 
can also be seen in an NJL calculation (iBuballal . [2005bt 
lForbesl . [2005h . The CFL-if° phase is a superfluid; it is 
a neutral insulator; all its quark modes are gapped (as 
long as Ms/{2fj.) < A); it breaks chiral symmetry. In all 
these respects it is similar to the CFL phase. Once we 
turn on small quark masses, different for all flavors, the 
SU{3)c+L+R symmetry of the CFL phase is reduced by 
explicit symmetry breaking to just U{1)q x C/(l)y, with 

Y a linear combination of a diagonal color generator and 
hypercharge. In the CFL-A'° phase, the kaon condensate 
breaks J7(l)y spontaneously. This modifies the spectrum 
of both quarks and Goldstone modes, and thus can affect 
transport properties. 



III. BELOW CFL DENSITIES 

As we discussed in the introduction (end of Sec. II. A[) 
and above fSec. lILBj) . at intermediate densities the CFL 
phase suffers from stresses induced by the strange quark 
mass, combined with beta-equilibration and neutrality 
requirements. It can only survive down to the tran- 
sition to nuclear matter (occurring at quark chemical 
potential /i = /^nuc) if the pairing is strong enough: 
roughly Acfl > Ms(/i„uc)^/2Ainuc, ignoring strong in- 
teraction corrections, which are presumably important 
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FIG. 3 (Color online) Free energy of various phases of dense 
3-flavor quark matter, assuming Acfl = 25 MeV. The ho- 
mogeneous phases are CFL and 2SC, their gapless analogs 
gCFL and g2SC, and the kaon-condensed phase CFL-K°. The 
true ground state must have a free energy below that of the 
gCFL phase, which is known to be unstable. The inhomoge- 
neous phases are curCFL-Jf" , which is CFL-if" with meson 
supercurrents, and 2PW, CubeX, and 2Cube45z, which are 
crystalline color superconducting phases. The transition from 
CFL-K° to curCFL-A'" is marked with a dot. In 2PW the 
condensate is a sum of only two plane waves. CubeX and 
2Cube45z involve more plane waves, their condensation en- 
ergies are larger but less reliably determined, so their curves 
should be assumed to have error bands comparable in size to 
the difference between them. 



in this regime. It is therefore quite possible that other 
pairing patterns occur at intermediate densities, and in 
this section we survey some of the possibilities that have 
been suggested. 

Fig. [3] shows a comparison of the free energies of some 
of these phases. We have chosen Acfl = 25 MeV, so 
there is a window of non-CFL pairing between nuclear 
density and the region where the CFL phase becomes 
stable. (For stronger pairing, Acfl ^ 100 MeV, there 
would be no such window.) The curves for the CFL, 2SC, 
gCFL, g2SC, and crystalline phases (2PW, CubeX and 
2Cube45z) are obtained from an NJL model as described 
in Sec. IVII The curves for the CFL-K'^ and meson su- 
percurrent (curCFL-X^) phases are calculated using the 
CFL effective theory with parameters chosen by match- 
ing to weak-coupling QCD, as described in Sec.|Vl except 
that the gap was chosen to match Acfl — 25 MeV. The 
phases displayed in Fig. [3] are discussed in the following 
sections. 



A. Two-flavor pairing: the 2SC phase 

After CFL, 2SC is the most straightforward less- 
symmetrically paired form of quark matter, and was one 
of the first patterns to be analyzed ( Alford et all . [l998l : 
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Bailin and Lov3 . Il979l . Il984 iBarroisl Il979t iRapp et all 
19981 ). In the 2SC phase, quarks with two out of three 



B. The unstable gapless phases 



colors (red and green, say) and two out of three flavors, 
pair in the standard BCS fashion. The flavors with the 
most phase space near their Fermi surfaces, namely u 
and d, are the ones that pair, leaving the strange and 
blue quarks unpaired (middle panel of Fig. [2|). Accord- 
ing to NJL models, if the coupling is we ak then there 
is no 2SC region in the phase diagram (jSteiner et all . 
I2OO2D . This can be understood by an expansion in pow- 
ers of Ms, which finds that the CFL^2SC transition 
occurs at the same point as th e 2SC— ^unpaired tran- 
sition , leaving no 2SC window ( Alford and Raiagopall 
|2002[ ) (this is the situation in Fig. [3]). However, NJL 
models with stronger coupling leave open the possibility 
of a 2S C window in the "non-CFL" r egion of the TJhase d i- 
agram (|Abuki and Kunihirol . l2006t iRiister et all l2005f) . 
(These calculations have to date not included the pos- 
sibility of meson current or crystalline color supercon- 
ducting phases, discussed below, that may prove more 
favorable.) 

The 2SC pairing pattern, corresponding to (j)^ = S^Sg 

in dHl), is {tpiCj5'4'j) oc A2scey3e"'^^ where the sym- 
metry breaking pattern, assuming massless up and down 
quarks, is 



[5f/(3)e] X SU{2)l X SU{2)b, x U{1)b x C/(1). 



3 [C/(1)q] 

[SUi2)rg] X SU{2)l X SU{2)r x [/(l)^ x U{1). 



3 [C^(1)q] 

.(9) 

using the same notation as in Eq. The unpaired 

massive strange quarks introduce a U{l)s symmetry. 
The color SU(3)c gauge symmetry is broken down 
to an SU{2)rg red-green gauge symmetry, whose con- 
finement distance rises exponentially with density, as 
expfconst /x/f q A)') dRischke et all 200ll) ( see al so 
(|Casalbuoni et all l2002bl : lOuved and Sanninol . 12001^ . 
An interesting feature of 2SC pairing is that no global 
symmetries are broken. The condensate is a singlet of 
the SU{2)l X SU{2)fj flavor symmetry, and baryon num- 
ber survives as _B, a linear combination of the original 
baryon number and the broken diagonal Tg color gener- 
ator. Elcctromagnetism, originally a linear combination 
of B, S, and I3 (isospin), survives as an unbroken linear 
combination Q of B, S, and I3. 2SC quark matter is 
therefore a color superconductor but is neither a super- 
fluid nor an electromagnetic superconductor, and there 
is no order parameter tha t distinguishes it fro m the un- 
paired phase or the QGP (jAlford ad Il998[) . With re- 
spect to the unbroken U{1)q gauge symmetry, the 2SC 
phase is a conductor not an insulator because some of 
the ungapped blue and strange quarks are Q-charged. 



As was noted in Sec. III.Bl and can be seen in Fig. [31 
the CFL phase becomes unstable when /i w ^AI^ / Aqfl- 
At this point the pairing in the gs-bd sector suffers the 
instability discussed in Sec. lI.El and it becomes energeti- 
cally favorable to convert gs quarks into bd quarks (both 
near their common Fermi momentum).'^ If we restrict 
ourselves to diquark condensates that are spatially ho- 
mogeneous, the result is a modification of the pairing in 
which there is still pairing in all the color-flavor channels 
that characterize CFL, but gs-bd Cooper pairing ceases 
to occur in a range of momenta near the Fermi sur face 
(lAlford et aLl . l2n04bL [20053 : iFukushima et aLl . l20n,^ . In 
this range of momenta there are bd quarks but no gs 
quarks, and quark modes at the edges of this range are 
ungapped, hence this is called a gapless phase ("gCFL"). 
Such a phenomenon w as first proposed for two fl avor 
quar k matter ("g2SC") ( Shovkovv and IIuan^ . l2003f ). see 



also ( Gubankova et all 12003 ). It has been confirmed 



in NJL analyses such as those in (lAbuki et all 2005 



Abuki and Kunihiro '. 2006: Alfor d et all l2005bl l2004allb 



2005cl:lFukushima et aL ,2005. : .Riister et aLl . |2004 12005^ . 



which predict that at densities too low for CFL pairing 
there will be gapless phases. 

In Fig. [3l where Acfl = 25 MeV, we see the transi- 
tion from CFL to gCFL at M^/n « 2Acfl = 50 MeV. 
(It is interesting to note that, whereas the CFL phase 
is a Q-insulator, the gCFL phase is a Q-conductor, be- 
cause it has a small electron density, balanced by un- 
paired bu quarks from a very thin momentum shell of 
broken bu-rs pairing; the CFL— *-gCFL transition is the 
analogue of an insulator to metal transition at which a 
"band" that was unfilled in the insulating phase drops 
below the Fermi energy, making the material a metal.) 
The gCFL phase then remains favored beyond the value 
M'^/H « 4AcFL = 100 MeV at which the CFL phase 
would become unfavored relative to comple tely unpaired 
quark matter ( Alford and Raiagopall . |2002| ) . 

However, it turns out that in QCD the gap- 
less p hases, both g2SC (ICiannakis and Renl . 
l2005al: ItRiang and Shovkovvl l2004bD and gCFL 

[2OOI), are 



(jCasalbuoni et all 



2005bl : iFukushimal . 
unstable at zero temperature. (Increasing the tempera- 
ture above a critical value removes the instability; the 
critical value varies dramatically between phases, from 



a fra ction of an MeV to of order 10 MeV ( Fukushimal . 
I2OO5I ).) The instability manifests itself in an imaginary 
Meissner mass rrij^j for some of the gluons. is the 

low-momentum current-current two-point function, and 



^ The onset of gaplessness occurs at the ^ at which ^{^ibd~ t^gn) = 
AcFL, as explained in Sec. II. El Note that in the CFL phase 
(i^bd ~ Mgs) = ^s/Mi twice its value in unpaired quark mat- 
ter because of the nonzero color chemical potential /ig oc M^//i 
required by color neutrality in the presence o f CFL pairing 
l lAlford and Raiagopall |2002| : ISteiner et al\ . |2002|) . 
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TO^^/(g^A'^) (where the strong interaction coupling is 
g) is the coefficient of the gradient term in the effective 
theory of small fluctuations around the ground-state 
condensate, so a negative value indicates an insta- 
bility towa rds sponta n eous breaking of tran s lationa l 
invar ia nce llFu kushima'. '2006t iHashimotol. 20061 'Huang , 
2006[ llida and Fukushimaj, l2006l : iReddva nd Rupalj, 
20051 ). Calculati o ns in a simple two-species model 



(lAlford and Wand . I2005D show that gapless charged 
fermionic modes generically lead to imaginary mj^j . 

The instability of the gapless phases indicates that 
there must be other phases of even lower free energy, that 
occur in their place in the phase diagram. The nature of 
those phases is not reliably determined at present; likely 
candidates are discussed below. 



C. Crystalline color superconductivity 

The Meissner instability of the gCFL phase points 
to a breaking of translational invariance, and crys- 
talline color superconductivity represents a possible res- 
olution of that instability. The bas ic idea, first pro- 
posed in condensed matter ph y sics (Fulde and FerreUl . 
Il964l : iLarkin and Ovchinnikovl Il965 ) and more re- 
cently analyz ed i n the context of color supercond uctiv- 
ity jAlford erall. l2001al : iBowers and Ra iagopal l2002t 
ICasalbuoni and Nardullif . 12004) . is to allow the different 
quark flavors to have different Fermi momenta, thus ac- 
commodating the stress of the strange quark mass, and to 
form Cooper pairs with nonzero momentum, each quark 
lying close to its respective Fermi surface. The price 
one must pay for this arrangement is that only fermions 
in certain regions on the Fermi surface can pair. Pairs 
with nonzero momenta chosen from some set of wave vec- 
tors Qd yield condensates that vary in position space like 
exp(iqo -x), forming a crystalline pattern whose Bra- 
vais lattice is the set of qa. 

Analyses to date have focused on u-d and u-s pair- 
ing, neglecting pairing of d and s because the separa- 
tion of their Fermi momenta is twice as large (Fig. 
If the (ud) condensate includes only pairs with a single 
nonzero momentum q, this means that in position space 
the condensate is a single plane-wave and means that in 
momentum space pairing is allowed on a single ring on 
the u Fermi surface and a single ring on the opposite 
side of the d Fermi surface. The simplest "crystalline" 
phase of three- fiay o r quark matter that has been anal yzed 
(jCasalbuoni et all l2005at iMannarelli et all l2006bt ) in- 
cludes two such single-plane wave condensates ("2PW"), 
one (ud) and one (us). The favored orientation of the 
two q's is parallel, keeping the two "pairing rings" on 
the u Fermi surface (fro m the (us) and (ud ) conde nsates) 
as far apart as possible ( Mannarelli et all l2006b[ ) . This 
simple pattern of pairing leaves much of the Fermi sur- 
faces unpaired, and it is much more favorable to choose 
a pattern in which the (us) and (ud) condensates each 
include pairs with more than one q-vector, thus more 



than one ring and more than one plane wave. Among 
such more realistic pairing patterns, the two that ap- 
pear most favorable have either four q's per condensate 
that together point at the eight corners of a cube in mo- 
mentum space ("CubeX") or eight q's per condensate 
that each point at the corners of separate cubes, ro- 
tated relative to each other by 45 degrees ("2Cube45z") 
( Raiagopal and Sharmal . [2"006bl ) . It has been shown that 
the chromom agnetic instabili t y is n o longer present in 



20061 ). The free energies 



these phases ([Ciminale et all 
of the 2PW, CubeX and 2Cube45z phases as calculated 
within an NJL model (see Sec. IVI[) are shown in Fig. [31 
The calculation is an expansion in powers of (A/Spp)'^ 
which in the CubeX and 2Cube45z phases turns out to 
be of order a tenth to a quarter. According to results 
obtained in a calculation done to third order in this ex- 
pansion parameter, the CubeX and 2Cube45z condensa- 
tion energies are large enough that one or other of them 
is favored over a wide range of M"^ / ^ as illustrated in 
Fig. [21 The uncertainty in each is of the same order as 
the difference between them, so one cannot yet say which 
is favored, but the overall scale is plausible (one would 
expect condensation energies an order of magnitude big- 
ger than that of the 2PW state). We discuss crystalline 
color superconductivity in greater detail in Sec. IVII 



D. Meson supercurrent ("curCFL-/('°") 

Kaon condensation alone does not remove the gap- 
less modes that occur in the CFL phase when Ms be- 
comes large enough, but it does affect the number of 
gapless modes and the onset value of Mg. In the CFL- 
phase, the electrically charged [hs) mode becomes 
gapless at / fJ. ~ 8 A/3 (compared to 2 A in the CFL 
phase), and the electricall y neutral (bd) mode becomes 



pnasej, and tne electricall y neutral (od) mode becomes 
gapless for M,^/m ~ 4 A (iKrvievski and Schafe^ . l2005t 
iKrvievski and Yamadal . 120051 ) . fl n an NJL model anal- 
vsis (jForbesl l2005l ) ."the charged mode in the CFL-A'° 
phase becomes gapless at M^/fJ. ~ 2.44A for A = 25 
MeV as in Fig. ^. The gapless CFL-K° phase has 
an instability which is similar to the instability of the 
gCFL phase. This instability can be viewed as a ten- 
dency towards spon taneous generation of Goldston e bo- 
son (kaon) currents (|Krvievskil[2005l : ISchafeil . l2006[ ). The 
currents correspond to a spatial modulation of the kaon 
condensate. There is no net transfer of any charge be- 
cause the Goldstone boson current is counterbalanced 
by a backflow of ungapped fermions. The meson super- 
current ground state is lower in energy than the CFL- 
K'^ state and t he m agnetic screening masses are real 
( Gerhold et al I I2007D . Because the ungapped fermion 
mode is electrically charged, both the magnitude of the 
Goldstone boson current needed to stabilize the phase 
and the magnitude of the resulting energy gain relative 
to the phase without a current are very small. Gold- 
stone boson currents can also be generated in the gCFL 
phase without K'^ condensation. In this case gauge in- 
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variance implies that the supercurrent state is equivalent 
to a single plane-wave LOFF state, but the analyses can 
be carried out in the limit that the gap is large compared 
to th e magnitude of the current ( Gerhold and Schafeij . 
I2OO6I ) . This analysis is valid near the onset of the gCFL 
phase, but not for larger mismatches, where states with 
multiple currents are favored. 



E. Single-flavor pairing 

If the stress due to the strange quark mass is large 
enough then there may be a range of quark matter 
densities where no pairing between different flavors is 
possible, whether spatially uniform or inhomogeneous. 
From Fig. [3] we can estimate that this will occur when 
M^/(/iAcFL) ^ 10, so it requires a large effective strange 
quark mass and/or small CFL pairing gap. The best 
available option in this case is Cooper pairing of each fla- 
vor with itself. Single-flavor pairing may also arise among 
the strange quarks in a 2SC phase, since they are not in- 
volved in two-flavor pairing. We will discuss these cases 
separately below. 

To maintain fermionic antisymmetry of the Cooper 
pair wavefunction, single-flavor pairing phases have to 
either be symmetric in color, which greatly weakens or 
eliminates the attractive interaction, or symmetric in 
Dirac indices, which compromises the uniform partic- 
ipation of the whole Fermi sphere. As a result, they 
have much lower critical temperatures than multi-flavor 
phases such as the CFL or 2SC phases, perhaps as 



large as a few MeV, more typically in the eV to many 
keV range fAlford et aLl.l2003lll998l:lBuballa aD. l 2003l : 
ISchafeil . 12000^ : ISchmittI I2OO5I: ISchmitt et all 12002^ 



Matter in which each flavor only pairs with itself 
has been studied using NJL models and weakly-coupled 
QCD. These calculations agree that the energetically fa- 
vored state is color-spi n-locke d (CSL) pairing for each 
flavor (iBailin and Lovd . Il979l : ISchafeii . l2000bl : ISchmltO . 
I2005D . CSL pairing involves all 3 colors, with the color 
direction of each Cooper pair correlated with its spin 
direction, breaking SU{3)c x S'0(3)rot S'0(3)c-i-rot- 
The phase is isotropic, with rotational symmetry sur- 
viving as a group of simultaneous spatial and color rota- 
tions. Other possible phases exhibiting spin-one, single- 
flavor, pairin g include the polar, planar, a nd A phases 



navor, pairin g include tne polar, planar, a nd A pnases 
described in ( Schafed . l2000bl:ISchmittl.l2005l) (for an NJL 



model treatment see () Alford et aZ.1 . 120031 )). Some of these 
phases exhibit point or line nodes in the energy gap at the 
Fermi surface, and hence do break rotational symmetry. 

If 2SC pairing occurs with strange quarks present, one 
might expect the strange quarks of all three colors to un- 
dergo CSL self-pairing, yielding an isotropic "2SC-I-CSL" 
pattern. However, the 2SC pattern breaks the color sym- 
metry, and in order to maintain color neutrality, a color 
chemical potential is generated, which splits the Fermi 
momentum of the blue strange quarks away from that 
of the red and green strange quarks. This is a small 



effect, but so is the CSL pairing gap, and NJL model 
calculations indicate that the color chemical potential 
typically destroys CSL p airing of the strange quarks 
(jAlford and Cowanl . [2001 . The system falls back on the 
next best alternative, which is spin-one pairing of the red 
and green strange quarks. 

Because their gaps and critical temperatures can range 
as low as the eV scale, single-flavor pairing phases in 
compact stars would appear relatively late in the life of 
the star, and might cause dramatic changes in its behav- 
ior. For example, unlike the CFL and 2SC phases, many 
single-flavor-paired p hases are electrical superconductors 
(jSchmitt et al so their appearance could signifi- 

cantly affect the magnetic field dynamics of the star. 



F. Gluon condensation 

In the 2SC phase (unlike in the CFL phase) the mag- 
netic instability arises at a lower value of the stress on 
the BCS pairing than that at which the onset of gap- 
less pairing occurs. In this 2SC regime, analyses done 
using a Ginzburg-Landau approach indicate that the in- 
stability can be cur ed by the appearance of a c hromo- 
electric condensate (iGorbar eit a/.l . l2006al l2007l l2006bt 
[Hashimoto and Miranskv . 2007 ). The 2SC condensate 
breaks the color group down to the SU{2)rg red-green 
subgroup, and five of the gluons become massive vector 
bosons via the Higgs mechanism. The new condensate in- 
volves some of these massive vector bosons, and because 
they transform non-trivially under SU {2)rg it now breaks 
that gauge symmetry. Because they are electrically 
charged vector particles, rotational symmetry is also bro- 
ken, and the phase is an electrica l superconductor. Alter- 
native ly, it has been suggested ( Ferrer and de la Inceral . 
l2007aD that the gluon condensate may be inhomogeneous 
with a large spontaneously-induced Q magnetic field. 



G. Secondary pairing 

Since the Meissner instability is generically associated 
with the presence of gapless fermionic modes, and the 
BCS mechanism implies that any gapless fermionic mode 
is unstable to Cooper pairing in the most attractive chan- 
nel, one may ask whether the instability could be re- 
solved without introducing spatial inhomogeneity sim- 
ply by "secondary pairing" of the gapless quasiparticles. 



which would then acq u ire th eir own gap ((Hong . .2005 
iHuang and Shovkovvl . |2003( ). Furthermore, there is a 
mode in the gCFL phase whose dispersion relation is well 
approximated as quadratic, e oc (fc — const)^, yielding a 
greatly increased density of states at low energy (diverg- 
ing as e~^/^), so its secondary pairing is much stronger 
than would be predicted by BCS theory: oc for 
an effective four-fermion coupling strength G, as com- 
pared witShthe standard BCS result A cx exp(— const /G) 
(lHong| . [2005h . This result is confirmed by an NJL study 
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in a two-spccics model ( Alford and Wan^l2006l ). but the 
secondary gap A^. was found to be still much smaller 
than the primary gap Ap, so it does not generically re- 
solve the magnetic instability (in the temperature range 
As < T < Ap, for example). 



H. Mixed phases 

Another way for a system to deal with a stress on 
its pairing pattern is to form a mixed phase, which is 
a charge-separated state consisting of positively and neg- 
atively charged domains which are neutral on average. 
The coexisting phases have a common pressure and a 
common value of the charge chemical potential which 
i s not equal to th e neutrality value for either phase 
(|Glendenninej . Il992l : iRavenhaU et al. I, [mi). The size of 
the domains is determined by a balance between surface 
tension (which favors large domains) and electric field 
energy (which favors small domains). Separation of color 
charge is expected to be suppressed by the very high 
energy cost of color electric fields, but electric charge 
separation is quite possible, and may occur at the in- 
terface between color-superconducti ng quark matter and 
nuclear matter ( Alford et aZ.I . [20"oibf ) and an interface be- 
tween quark matter an d the vacuum (|Alford et aZ.1 . 120061 : 
iJaikumar et aLl . [2006al ) . iust as it occur s at interfaces be- 
tween nuclear matter and a nucleon gas ( Ravenhall et all . 
Il983l ). Mixed phases are a generic phenomenon, since, 
in the approximation where Coulomb energy costs are 
neglected, any phase can always lower its free energy 
density by becoming charged (this follows from the fact 
that free energies are concave functions of chemical po- 
tentials). In this approximation, if two phases A and B 
can coexist at the same pressure with opposite charge 
densities then such a mixture will always be favored over 
a uniform neutral ph ase of either A or B. For a pedagog- 
ical discussion, see (lAlford et all l2005d) . Surface and 
Coulomb energy costs can cancel this energy advantage, 
however, and have to be calculated on a case-by-case ba- 
sis. 

In quark matter it has been found that as long as 
we require local color neutrality such mixed phases are 



not the favored respo nse to the st r ess impo sed by the 
strange quark mass (lAlford et all l2004allbl ). Phases 



i nvolving color charge separation have been studied 
( Neumann et all l2003l ) but it seems likely that the en- 
ergy cost of the color-electric fields will disfavor them. 



I. Relation to cold atomic gases 

An interesting class of systems in which stressed 
superconductivity can be studied experimentally is 
trapped atomic gases in which two different hyperfine 
states ("specie s ") of the atom pair with each other 
(jGiorgini et al I l200it ). This is a useful experimental 
model because the stress and interaction strength are 



both under experimental control, unlike quark matter 
where one physical variable (fi) controls both the cou- 
pling strength and the stress. The atomic pairing stress 
can be adjusted by changing the relative number of atoms 
of the two species ( "polarization" ). The scattering length 
of the atoms can be controlled using Feshbach reso- 
nances, making it possible to vary the strength of the 
inter-atomic attraction from weak (where BCS pairing 
occurs) through the unitarity limit (where a bound state 
forms) to strong (Bose- Einstein condensation of diatomic 
molecules) . 



The theoretical expectation is that, in the weak cou- 
pling limit, there will be BCS pairing as long as Sfi, 
the chemical potential difference between the species, is 
small enough. The BCS phase is unpolarized because 
the Fermi surfaces are locked together. A first-order 
transition from BCS to crystalline (LOFF) pairing is 
expected at 6iJ = A0/V2, where Ap is the BCS gap 
at (5/i = 0; then at Sur. a continuous t ransition to the 
unpaired phase ( Chandrasckhar*. '1962"; 'Clogston, 'l96 
iFulde and Ferrell lT964t Xarkin and Ovchinnikov, .196 
For the single plane wave LOFF state 6fic — 0.754Ao, 
but for multiple plane wave states Sfic may be larger. 



Experiments with cold trapped atoms near the uni- 
tary limit (strong coupling) have seen phase separation 
between an unp olarized su pe rfluid a nd a polarized 
normal state (IChin et al. I, I2OO6I : [Partridge et~d\ . 
I2OO6I : IZwierlein et all l2006^ . If one ignores the crys- 
talline phase (pe rha ps onl y favored at weak coupling 
(iMannarem et all l2006at ISheehv and Radzihovskvl 



l2006l . T2007lV r this is consistent with the theoretical 
expectation for the BCS regime: the net polarization 
forces the system to phase separate, yielding a mixture 
of BCS and unpaired phases wit h 6ii fixed at t h e first 
order transition be t ween them ( Bedaque et all 120031 : 
ICarlson and Reddvl . [200I). It remains an exciting 
possibility that crystalline superconducting (LOFF) 
phases of cold atoms may be observed: this may require 
experiments closer to the BCS regime. 



In the strong coupling limit the superfluid consists of 
tightly bound molecules. Adding an extra atom requires 
energy A. For \6pi\ > A the atomic gas is a homogeneous 
mixture of an unpolarized superfluid and a fully polarized 
Fermi gas, so the system is a stable gapless superfluid. 
This means that in strong coupling polarization can be 
carried by a gapless superfluid, whereas in weak coupling 
even a small amount of polarization leads to the appear- 
ance of a mixed BCS/LOFF phase. It is not known what 
happens at intermediate coupling, but one possibility is 
that the gapless superfluid a nd the LOFF phase are con - 
nected by a phase transition ( Son and Stephanovl .l2006f). 
This transition would correspond to a magnetic instabil- 
ity of the gapless superfluid. 
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IV. WEAK-COUPLING QCD CALCULATIONS 

We have asserted in Sees. U and |lT] that at sufficiently 
high densities it is possible to do controlled calculations 
of properties of CFL quark matter directly from the QCD 
Lagrangian. We describe how to do such calculations in 
this section. We shall focus on the calculation of the gap 
parameter, but shall also discuss the critical temperature 
Tc for the transition from the CFL phase to the quark- 
giuon plasma and the Meissner and Debye masses that 
control color-magnetic and color-electric effects in the 
CFL phase. Phenomena that are governed by the mass- 
less Goldstone bosons and/or the light pseudo-Goldstonc 
bosons are most naturally described by first constructing 
the appropriate effective theory and then, if at sufficiently 
high densities, calculating its parameters directly from 
the QCD Lagrangian. We defer these analyses to Sec.lVl 

Although the weak-coupling calculations that we de- 
scribe in this section are only directly applicable in the 
CFL phase, we shall present them in a sufficiently gen- 
eral formalism that they can be applied to other spatially 
homogeneous phases also, including for example the 2SC 
and CSL phases. These phases can be analyzed at weak- 
coupling either just by ansatz, or by introducing such a 
large strange quark mass that CFL pairing is disfavored 
even at enormous densities. Such calculations provide 
insights into the properties of these phases, even though 
they do not occur in the QCD phase diagram at high 
enough densities for a weak-coupling approach to be ap- 
plicable. To keep our notation general, we shall refer to 
the gap parameter as A; in the CFL phase, Acfl = A. 

We shall see that at weak coupling the expansion pa- 
rameter that controls the calculation of log(A/^) is at 
best g, certainly not g^. (The leading term is of order 
1/g; the \ogg and 5° terms have also been calculated. 
The O{glogg) and 0{g) terms arc nonzero, and have 
not yet been calculated. Beyond 0{g), it is possible that 
fractional powers of g may arise in the series.) We there- 
fore expect the weak-coupling calculations to be quan- 
titatively reliable only at densities for which g{fi) < 1, 
which corresponds to densities many orders of magni- 
tude greater than that at the centers of neutron stars . 
Indeed, it has been shown ( Raiagopal and Shusteij . [2000l ) 
that some of the 0{g) contributions start to decline in 
magnitude relative to the 17° term only for g{^) < 0.8 
which corresponds, via the two-loop QCD beta function, 
to /i > 10® MeV meaning densities 15-16 orders of mag- 
nitude greater than those at the centers of compact stars. 
The reader may therefore be tempted to sec this section 
as academic. From a theoretical point of view, it is ex- 
ceptional to have an instance where the properties of a 
superconducting phase can be calculated rigorously from 
a fundamental short-distance theory, making this explo- 
ration a worthy pursuit even if academic. From a prac- 
tical point of view, the quantitative understanding that 
we derive from calculations reviewed in this section pro- 
vides a completely solid foundation from which we can 
extrapolate downwards in fi. The effective field theory 



described in Sec. |V] gives us a well-defined way of doing 
so as long as we stay within the CFL phase, meaning that 
we can come down from /i > 10® MeV all the way down 
to ^ ^ Mg/(2AcFL)- Finally, we shall gain qualitative 
insights into the CFL phase and other color supercon- 
ducting phases, insights that guide our thinking at lower 
densities. 

The QCD Lagrangian is given by 



C = -tpii-ff^Dfj, + fl'jo - m)^ - 



(10) 



Here, ij) is the quark spinor in Dirac, color, and flavor 
space, i.e., a 4A^c./V/-component spinor, and -0 = V'Wo- 
The covariant derivative acting on the fermion field is 
= dfj, + igTaAlj^, where g is the strong coupling con- 
stant, A^;^ are the gauge fields, T'' = A°/2 (a = 1, . . . , 8) 
are the generators of the gauge group SU{3)c, and A'' 
are the Cell-Mann matrices. The field strength tensor is 
Gi^" = d^,Al-d^Al + gr'"'AlAl with the 5[/(3)c struc- 
ture constants f"-^'^. The chemical potential jl and the 
quark mass rh = diag(m„, m^, nis) are diagonal matrices 
in flavor space. If weak interactions are taken into ac- 
count flavor is no longer conserved and there are only two 
chemical potentials, one for quark (baryon) number, ^, 
and one for electric charge, //g. At the very high densities 
of interest in this section, the constituent quark masses 
are essentially the same as the current quark masses m„ , 
rud and ms meaning that we need not distinguish be- 
tween them. Furthermore, at asymptotic densities we 
can neglect even the strange quark mass, so throughout 
most of this section we shall set m„ = rud ~ iris = 0. 

If the coupling is small then the natural starting point 
is a free Fermi gas of quarks. In a degenerate quark gas 
all states with momenta p < pp = (/i^ — m^)^/^ are oc- 
cupied, and all states with p > pp are empty. Because of 
Pauli-blocking, interactions mainly modify states in the 
vicinity of the Fermi surface. Since the Fermi momen- 
tum is large, typical interactions between quarks near 
the Fermi surface involve large momentum transfer and 
are governed by the weak coupling ^(/i). Interactions in 
which quarks scatter by only a small angle involve only 
a small momentum transfer and are therefore potentially 
dangerous. However, small momenta correspond to large 
distances, and medium modifications of the exchanged 
gluons are therefore important. In a dense medium, elec- 
tric gluons are Debye screened at momenta q ^ gfx. The 
dominant interaction for momenta below the screening 
scale is due to unscreened, almost static, magnetic glu- 
ons. In a hot quark-gluon gas, interactions between mag- 
netic gluons become nonperturbative for momenta less 
than g^T. This phenomenon does not take place in a 
very dense quark liquid, and gluon exchanges with arbi- 
trarily small momenta remain perturbative. On a quali- 
tative level this can be attributed to the absence of Bose 
enhancement factors in soft gluon propagators. A more 
detailed explanation will be given in Sec. IV. A. 21 The 
unscreened magnetic interactions nevertheless make the 
fluid a "non- Fermi liquid" at temperatures above the crit- 
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ical temperature for cofor superconductivity. We shall 
discuss this also in Sec. IV. A. 21 where we shall see that 
these non-Fermi liquid effects do not spoil the basic logic 
of the BCS argument that diquark condensation must oc- 
cur in the presence of an attractive interaction, but are 
crucial in the calculation of the gap that results. 



A. The gap equation 

As discussed in See. II. B| any attractive interaction 
in a many-fermion system leads to Cooper pairing. 
QCD at high density provides an attractive interaction 
via one-gluon exchange. In terms of quark represen- 
tations of SU{3)c, the attractive channel is the anti- 
symmetric anti-triplet 3a-, appearing by "pairing" two 
color triplets: 3 (8) 3 = 3^ © ds- Consequently, only 
quarks of different colors form Cooper pairs. There 
is an induced pairing in the symmetric s extet channe l 
6,9- H ow ever, this pairing is ni u ch wea ker (lAlford et al\ . 
IQQQaHbl: iPisarski and Rischkd. Il999d : ISchafed . l2000a : 
Shovkovv and Wiiewardhanal . Il999l ) , and we shall largely 
neglect it in the following. As in an electronic supercon- 
ductor, Cooper pairing results in an energy gap in the 
quasiparticle excitation spectrum. Its magnitude at zero 
temperature A is crucial for the phenomenology of a su- 
perconductor. In addition, it also sets the scale for the 
critical temperature Tc of the phase transition which can 
be expected to be of the same order as A (in BCS theory, 
Tc = 0.57A). Over the course of the next five subsections, 
we shall discuss the QCD gap equation, which is used to 
determine both A and Tc. 

Our starting point is the partition function 



(11) 



with the action S — J d'^x C and the Lagrangian (fTO)l . In 
the following we shall only sketch the derivation of the 
gap equat ion. De t ails fol l owing the same lines can be 
found in jManue ll. [2"000bl: iPisarski a nd Risc hki. Il999b . 
2000al :[Risehke. 200"4l: ISchmittl . l2004 i2005: Sc hmitt et al 



2002D . 



We begin by introducing Nambu-Gorkov spinors. This 
additional two-dimensional structure proves convenient 
in the theoretical description of a superconductor or 
a superfluid, see for ins tance ( Abrikosov et all Il963t 
iFetter and Waleckal . fl97l[) . It allows for the introduction 
of a source that couples to quark bilinears (as opposed 
to quark-anti-quark bilinears). Spontaneous symmetry 
breaking is realized by taking the limit of a vanishing 
source. The Nambu-Gorkov basis is given by 



* = 



* = (V',V'c) 



(12) 



where ipc = is the charge-conjugate spinor, ob- 

tained by multiplication with the charge conjugation ma- 
trix C = 17^7*^. In a free fermion system, the new basis 



is a pure doubling of degrees of freedom with the inverse 
fermion propagator consisting of the original free propa- 
gators. 



nan 





[Gal 



(13) 



where [Gt]-HX,Y) = (z7'^9^ ± ^70) (5W(X - y). 
Here and in the following capital letters denote four- 
vectors, e.g., X = {xo,x.). The effect of a nonzero 
diquark condensate can now be taken into account 
through adding a suitable source term to the ac- 
tion and computing the effective action F as a func- 
tio nal of the ghio n and fermio n pr opagators D and 
S (lAbukil. ,2003: MTra^sky aU [2OOI; Risclikg 
iRiister and Rischkd . |2004 ISchmitll . |2004 lTakagi[|2003D : 



-^TTlogD-'-^Tr{D^'D-l) 
-iTrlog^-i+ iTr(5o-i5-l) 



(14) 



This functional is called the "2PI effective action" since 
the contribution T2\ p,S] consist s of all two-part i cle ir- 
reducible diagrams (iBavml . Il962l : ICornwall et all . Il974l : 
iLuttinger and Wardl . Il96df i. This formalism is particu- 
larly suitable for studying spontaneous symmetry break- 
ing in a self-consistent way. The ground state of the 
system is obtained by finding the stationary point of the 
effective action. The stationarity conditions yield Dyson- 
Schwinger equations for the gauge boson and fermion 
propagators. 



n. 



(15a) 
(15b) 



where the gluon and fermion self-energies are the func- 
tional derivatives of F2 at the stationary point. 



.5T 



(5F2 



n=-2TT^' S = 2-^- (16) 



5D 



5S 



Writing the second of these equations as F2[5] = 
(1/4) Tr(S S'), we can then use the Dyson-Schwinger 
equation (|15b|) to evaluate the fermionic part of the effec- 
tive action at the stationary point, obtaining the pressure 

P = ^Tr log - - S^' S) . (17) 

We shall return to this expression for the pressure in 
Sec. lITO 

Here, we proceed to analyze the Dyson- Sehwinger 
equation (|15bp for the fermion propagator. We denote 
the entries of the 2x2 matrix S in Nambu-Gorkov space 
as 



(18) 
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where the off-diagonal elements are related via $^ ~ 
7o(tf>^)^7o- One can invert the Dyson-Schwinger equa- 
tion formally to obtain the full fermion propagator in 
the form 



S = 



G+ F- 
F+ G- 



(19) 



where the fermion propagators for quasiparticles and 
charge-conjugate quasiparticles are 



(20) 

and the so-called anomalous propagators, typical for a 
superconducting system, are given by 



(21) 



They can be thought of as describing the propagation of 
a charge-conjugate particle (i.e., a hole) with propagator 
([Gg"]~^ + E~)~^ that is converted into a particle with 
propagator G^ , via the condensate (Or, a particle 
that is converted into a hole via the condensate.) The 
essence of superconductivity or superfluidity is the ex- 
istence of a difermion condensate that makes the quasi- 
particle excitations superpositions of elementary states 
with fermion-numbcr ...,—5,-3,-1,1,3,5,...; we see 
the formalism accommodating this phenomenon here. 

We shall approximate r2 by only taking into account 
two-loop diagrams. Upon taking the functional deriva- 
tive with respect to S, this corresponds to a one-loop 
self-energy E. We show E diagrammatically in Fig. ID 
We shall argue later that this approximation is sufficient 
to calculate log(A/yu) up to terms of order g'^. Upon 
making this approximation, the gap equation takes the 
form shown in the lower panel of Fig. |4l namely 



$+(X) 



Q 



F+{Q)YnDf^{K-Q), (22 



in momentum space, where D'^'l{K — Q) is the gluon 
propagator. 

Note that in the derivation of the gap equation we have 
assumed the system to be translationally invariant. This 
assumption fails for crystalline color superconductors, see 
Sec. lVIl There has been some work on analyzing a partic- 
ularly simple crystalline phase in QCD at asymptoticall y 
high densities and weak coupling ( Leibovich et allhoOlh . 
but the formalism we are employing does not allow us to 
incorporate it into our presentation and, anyway, this 
subject remains to date largely unexplored. 



B. Quasiparticle excitations 

Before we proceed with solving the gap equation, it 
is worthwhile to derive the dispersion relations for the 
fermionic quasiparticle excitations in a color supercon- 
ductor. That is, we suppose that the gap parameter(s) 



E = 






FIG. 4 Upper panel: Diagrammatic representation of the 
quark self-energy in Nambu-Gorkov space. Curly lines cor- 
respond to the gluon propagator D. The quasiparticle propa- 
gators and G~ are denoted by double lines with an arrow 
pointing to the left and right, respectively. The anomalous 
propagators in the off-diagonal entries are drawn accord- 
ing to their structure given in Eq. (|2ip : thin lines correspond 
to the term ([Gj]~^ -I- E^)~^, while the cross-hatched and 
hatched circles denote the gap matrices and $~, respec- 
tively. Lower panel: The QCD gap equation p2|l is obtained 
by equating with the lower left entry of the self-energy de- 
picted in the upper panel (the other off-diagonal component 
yields an equivalent equation for $^). 



A have been obtained in the manner that we shall de- 
scribe below and ask what are the consequences for the 
quasiparticle dispersion relations. Based on experience 
with ordinary superconductors or superfluids, we expect 
(and shall find) gaps in the dispersion relations for the 
fermionic quasiparticles. We may also expect that in 
some color superconducting phases, quasiparticles with 
different colors and flavors, or different linear combina- 
tions of color and flavor, differ in their gaps and disper- 
sion relations. Indeed, some gaps may vanish or may be 
nonzero only in certain directions in momentum space. 

The quasiparticle dispersion relations are encoded 
within the anomalous self-energy defined in (|18p . 
which satisfies the gap equation ([22|) . We shall assume 
that can be written in the form 



(23) 



where is a matrix in color, flavor and Dirac space, 

and A^^"^^ = (1 + 6707 • k)/2 are projectors onto states of 
positive (e = +) or negative (e = — ) energy. The corre- 
sponding gap functions are denoted as A*^^) {K) and will 
be determined by the gap equation. Here and in the fol- 
lowing the energy superscript is denoted in parentheses 
to distinguish it from the superscript that denotes com- 
ponents in Nambu-Gorkov space. In our presentation we 
shall assume that A4 is momentum- independent, corre- 
sponding to a condensate of Cooper pairs with angular 
momentum J = 0, but the formalism can easily be ex- 
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tended to allow a momentum-dependent A^t as required 
for example in the analysis of the CSL phase and we shall 
quote results for this case also. Note that in Eq. ([23|) we 
are assuming that every nonzero entry in M is associated 
with the same gap functions A^*^^; the formalism would 
have to be generalized to analyze phases in which there is 
more than one independent gap function, as for example 
in the gCFL phase. 

We shall analyze color superconducting phases whose 
color, flavor and Dirac structure takes the form 



M 



a/3 



(24) 



where the 75 Dirac structure selects a positive parity con- 
densate, where, as described in Sees. |I] and [TTl the anti- 
symmetric color matrix is favored since QCD is attrac- 
tive in this channel and the antisymmetric flavor matrix 
is then required, and where ^ is a 3 x 3 matrix. We 
note that because the full flavor symmetry is the chiral 
SU{3)l X SU{3)r symmetry, the matrix (j) is actually a 
pair In this section we shall assume (j)^ = (j)^ . 

The case (f)^' 7^ (f)^ , which corresponds to a meson con- 
densate in the CFL phase, is discussed in Sec. IV. CI 

The excitation spectrum is given by the poles of the 
propagator S in p^ . (We shall see that the diagonal 
and the off-diagonal entries in S have the same poles.) 
It will turn out that the Hermitian matrix A4Ai^ deter- 
mines which quasiparticles are gapped and determines 
the ratios among the magnitudes of (possibly) different 
gaps. It is convenient to write this matrix via its spectral 
representation 



(25) 



where are the eigenvalues and Vr the projectors onto 
the corresponding eigenstates. 

The final preparation that we must discuss prior to 
computing the propagator is that we appro ximate the di- 
agona l elements of the quark self-energy as (Brown et all . 
I2000d: iGerhold and Rebhli^ . l2005t Im anue 1 l2000bD 



g 



ASe'^ml 

'log _2,9 



(26) 



where = N jg^ / {Qii^) is the square of the effective 
gluon mass at finite density, and e is the Euler constant. 
The expression (j26p is the low energy approximation to 
the one- loop self-energy, valid for /cq ~ A <C m^, for the 
positive energy (e = 4-) states. Taking the low energy 
approximation to and neglecting the self-energy cor- 
rection to the negative energy states will prove sufficient 
to determine log(A/^) up to order (7°. 

With all the groundwork in place, we now insert Go, 
S from (US]), and $ from ^ into ^ and and 
hence (fTg]). and use ([25|l to simplify the result. We find 
that the diagonal entries in the fermion propagator S are 



given by 



p. A 



[k^/zi-){k^)Y -[ 

while the anomalous propagators are 



2 r.Wi2 



(27) 



, (28a) 



t::. [fco/z(-)(fco)] - [ 



In writing these expressions, we have defined the wave 
function renormalization factor 



(29) 



for the positive energy e = + components, originating 
from the self-energy ((26|) . (By neglecting the negative en- 
ergy contribution to in (|26p . we are setting the nega- 
tive energy wave function renormalization Z^^\kQ) = 1.) 
We have furthermore defined 



y'(efc-M)' + A.|A(-)|2. (30) 



The r'th quasipartiele and antiquasiparticle energies are 
then given by solving fcp = ^^'^H^o)^^'^^ ^^i' fco- To lead- 



k,r 

ing order in g, wave function renormalization can be ne- 
glected and the quasipartiele and antiquasiparticle ener- 
gies are given by the e'j^l themselves. We see from ([30)) 
that the antiparticles have e > ^ — in fact, for k near 
fi they have e ~ 2fi. They therefore never play an im- 
portant role at high density. This justifies our neglect of 
the negative energy and hence of Z'^^-'. And, it jus- 
tifies the further simplification that we shall henceforth 
employ, setting the antiparticle gap to zero, A^"-* = 0, 
and denoting A = A'^+'. We shall also use the notation 

Z{ko) = Z^+\ko) and ek,r = f-i^J ■ We then see that the 
minimum value of Ck^r occurs at the Fermi surface, where 
k = fj,, and is given by -^/A^A which is conventionally 
referred to as the gap, again neglecting wave function 
renormalization. We see that although we must solve 
the gap equation in order to determine the magnitude of 
the gap parameter A, as we will do in Sees. IIV.DI and 
IIV.EI the ratios among the actual gaps in the quasipar- 
tiele spectra that result are determined entirely by the 
Ar's, namely the eigenvalues of MM"^ . 

We close this subsection by evaluating the pattern of 
quasipartiele gaps explicitly for the CFL and 2SC phases. 
We list the order parameters (f)^ , eigenvalues A^ , and cor- 
responding projectors Vr for these two phases in Table 
m In the CFL phase, one finds the eigenvalues Ai = 4 
with degeneracy Tr[7^i] = 1 and A2 = 1 with degeneracy 
Tr[7'2] = 8. This means that all nine quasiparticles are 
gapped. There is an octet with gap A. and a singlet with 
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gap 2A. The octet Cooper pairs are gu-rd, bd-gs, bu-rs, 
as well as two linear combinations of the three quarks 
ru-gd-bs. (Here, gu refers to a green up quark, etc.) The 
singlet Cooper pair with twice the gap is the remaining 
orthogonal combination of ru-gd-bs. In the 2SC phase, 
on the other hand, we find four quasiparticles with Ai = 1 
and hence gap A and 5 quasiparticles with A2 = that 
are unpaired. The gapped quasiparticles involve the first 
two colors, red and green, and the first two flavors, up 
and down. The Cooper pairs have color-flavor structure 
ru-gd and gu-rd. (Note that all these color-flavor combi- 
nations depend on the chosen basis of the color and fla- 
vor (anti)triplets. This basis is fixed by Eq. ([24]); apply- 
ing color (flavor) rotations to e°'P^ i^ijB) would change 
the basis and yield different, physically equivalent, color- 
flavor combinations for the 2SC and CFL phases.) 

The formalism of this section can easily be applied to 
patterns of pairing in which depends on the direc- 
tion of the quark momentum k. Such phases arise if the 
Cooper pairs carry total angular momentum J = 1. This 
allows for pairing between quarks of the same flavor, as 
discussed in Sec. IIII.EI Depending on the specific struc- 
ture of Mk, the eigenvalues may become momentum 
dependent and lead to nodes in the gap function along 
certain directions in momentum space. 



C. Pressure and condensation energy 

We can now return to our expression pT|) for the pres- 
sure P (equivalently, the thermodynamic potential since 
51 = —P) for a color superconductor and use the re- 
sults of Sec. IIV.BI to evaluate it for a superconducting 
phase of the form (|24p . We first substitute the expres- 
sions (f27| and ((28|) for the fermion propagator (fT9|) in the 
pressure ([T7| . In order to obtain a result that is valid 
at both nonzero and zero temperature, it is then conve- 
nient to switch to Euclidean space, and perform the sum 
over Matsubara frequencies. Upon doing the trace over 
Nambu-Gorkov, color, flavor, and Dirac space we find 



P 



Xr |A 



"k.r 



2e 



, , tanh , 

k.r \ / 



(31) 



Including the effects of wave function renormalization 
would modify this expression at order g. In most con- 
texts, we shall only consider the pressure (pij) at zero 
temperature. In this case, with '^'"[■'^'■] ~ ^cNf, 



P = NrN 



127r2 



SP 



(32) 



where we denote the pressure difference of the color- 
superconducting phase compared to the unpaired phase 
by SP. If we make the simplifying assumption (corrected 



in the next subsection) that the gap function is a con- 
stant in momentum space in the vicinity of the Fermi 
surface, we find the easily interpretable result 



(33) 



At T = this quantity is the condensation energy density 
of the color-superconducting state. The fact that 6P > 
implies that the superconducting state is favored relative 
to the normal phase. We observe that SP is proportional 
to the sum of the energy gap squared of the r-th branch, 
multiplied by the corresponding degeneracy Tr[Pr]- 

We can use the result ([55]) to understand how to com- 
pare the favorability of different patterns of color super- 
conducting pairing: the phase with lowest free energy 
(highest SP) is favored. As an example, in the CFL 
phase SP = {^i'^ /{47:^)){8 • 1 + 1 • 4)A^pL while in the 
2SC phase SP = {fi^ /{4,tt^)){4: -1 + 0- 5)Alsc suggesting 
that the CFL phase is favored. (We shall make this con- 
clusion firm in Sec. IIV.EI where we shall find that Acfl 
is smaller than A2SC but only by a factor of 2^^^. This 
factor will also turn out to be determined entirely by the 
Ar's and Tr[Pr]'s.) 

In principle, in order to generalize the conclusion that 
the CFL phase is favored one has to compare the con- 
densation energies of all possible phases described by the 
order parameter M in Eq. ([M]) . This is difhcult because 
(j) is an arbitrary complex 3x3 matrix. At asymptotic 
densities, however, we can neglect the strange quark mass 
and treat the quarks as degenerate in mass. The result- 
ing SU{3)c X SU{3)f symmetry simplifies the task. (/ is 
L ov R for 0i or 0ij.) The matrix (j) transforms under 
color-flavor rotations as ^ U'^ipV with U S SU{3)c, 
V £ SU{3)f. This means that two order parameters 
and (pV describe the same physics. Now note that 
for any (j) there exists a transformation ([/, V) such that 
is diagonal. Therefore, we need consider only di- 
agonal matrices </>. Choosing all diagonal elements to be 
nonzero corresponds to the maximum number of gapped 
quasiparticles. Hence, once wc show (below) that A2SC is 
not much larger than Acfl it is easy to understand that 
the CFL phase with = 1, yielding an order parameter 
that is invariant under the largest possible subgroup of 
the original symmetries, is the ground state at asymptot- 
ically large densities. 

At lower densities, the flavor symmetry is explicitly 
broken by the mass of the strange quark (the symme- 
tries are further broken by different chemical potentials 
due to neutrality constraints). In this case, the above 
argument fails and non-diagonal mat rices (j) become pos- 
sible candidates for the grouii d state ( Malekzadehl l2006t 
iRaiagopal and Schmittl . l2006[ ). 



D. Weak coupling solution of the gap equation 

We are now in a position to solve the QCD gap equa- 
tion (j22p for an order parameter with a given matrix 
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phase 


4>A 


Ai 


Aa 






CFL 


Si 


4 (1-fold) 


1 (8-fold) 






2SC 




1 (4-fold) 


(5-fold) 


(5c«,a-'5a3'5/33)(5'^-<5''5^') 





TABLE I Color- flavor structure of CFL and 2SC phases: Order parameters (f)^, eigenvalues Ar of the matrix MJvl\ and 
corresponding projectors Vr, derived from Eq. (|24|l . Color (flavor) indices are denoted a, P, {i, j). 



structure M. The matrix structure of the gap equation 
([22)) is handled by multiplying both sides of the equation 
by 7WtAj^+^ and taking the trace over color, flavor, and 
Dirac indices. 

The gap equation is sensitive to gluon modes with 
small momentum (p <^ nig) and even smaller energy 
(po ^ p^/rrv^ ^ p), meaning that medium effects in 
the gluon propagator have to be taken into account. In 
the low momentum limit, the gluon propagator takes 
on the st andard hard-dense loop approximation form 
(jBraaten a nd Pisarski, 1992), which we shall give below 
in Eqs. ([55]) and ([57)1 upon simplifying it as appropriate 
for Po ^ P- In order to obtain log(A//x) to order 
it suffices to keep only the leading terms in the propa- 
gator in the po p limit. We shall work in Coulomb 
gauge. Gauge independence of th e gap in a generalized 
Coul omb gauge was established in ( Pisarski and Rischkd . 
I2OO2D . an d a more formal proo f of gauge invariance wa s 
given in (|Gerhold and Rebhanl . l2003t iHou et all |2004[) . 
The gap equation reads 

A,, = ^ /-^ y Z{e,,)^ tanhf^) 

X 7^,(p)rot,(k,q)+Ab,eg,s,efc,.)r/(k,q)] ,(34) 

where we have abbreviated P = K — Q and have de- 
noted the gap function on the quasiparticle mass shell 
by Ak,r = A(efc_j.,k). We have denoted the traces over 
color, flavor, and Dirac space by 



Tr 



^;.(k,q) 



<t a(+) 



Tr 



(35) 

and T/(k,q) = ~{S''^ - p'p' ) (k, q) . The two terms 
inside the square bracket in Eq. (p4)) correspond to the 
contributions from electric and magnetic gluons. The 
dominant contribution comes from almost static gluons 
with Po <^ p. The static electric and almost static mag- 
netic gluon propagator give 



Di{p) 



M4(e + e')2 



(36) 
-e'), (37) 



where Mg = 



(37r/4)r 



With the gap equation now 



We can solve ([M)) for the zero temperature gap A on 
the Fermi surface. Or, we can solve for T in the A ^ 
limit, thus obtaining the critical temperature Tc- Solving 
for A, wc find that it has a weak coupling expansion of 
the form 



— - &olog(.g) - bo 
9 



(38) 



In our treatment of the fcrmion propagator, the gluon 
propagator, and in our truncation of the self-energy in 
Fig. m to one loop (for example neglecting vertex renor- 
malization) we have been careful to keep all effects that 
contribute to 60 1 but we have neglected many that con- 
tribute at order g\-Ogg and g. The formalism that we 
have presented can be used to evaluate ho and &o, 
and we shall describe the results in Sec. IIV.EI 

Before turning to quantitative results, it is worth high- 
lighting the origin and the importance of the leading 
— l/g behavior in ([35]), namely the fact that (A//x) '-^ 
exp(— constant/g). If in the gap equation of Fig. [H we 
were to replace the exchanged gluon by a contact inter- 
action, we would obtain a gap equation of the form 



A (X 5M 



A 



VF+A^ 



(39) 



with ^ = k — fi. This always has the solution A = 0; to 
seek nonzero solutions, we cancel A from both sides of 
the equation. Then, if A were 0, the remaining integral 
would diverge logarithmically at small ^. Therefore, wc 
find a nonzero A for any positive nonzero g no matter 
how small, with A oc exp(— constant/g^). This is the 
original ECS argument for superconductivity as a con- 
sequence of an attractive interaction at a Fermi surface. 
However, once we restore the gluon propagator the ar- 
gument is modified. The crucial point is that magnetic 
gluon exchange is an unscreened long-range interaction, 
meaning that the angular integral will diverge logarithmi- 
cally at forward scattering in the absence of any mech- 
anism that screens the magnetic interaction. The gap 
equation therefore takes the form 



A cx 



de- 



(40) 



stated fully explicitly, all that remains is to solve it. 



where 9 is the angle between the external momentum k 
and the loop momentum q and where 6 is some quantity 
with the dimensions of mass that cuts off the logarith- 
mic coUinear divergence of the angular integral. In the 
superconducting phase this divergence will at the least 
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be cut off by the Meissner effect, which screens gluon 
modes with p < A (since the Cooper pairs have size 
1/A) giving S A. This yields A g^A{\ogA)^ and 
hence a nonzero gap A ~ exp(— constant/g). This con- 
sequence of the long-range nature of the magnetic gluq n 
exchange was first discovered by Barrois (|Barroislll979t ). 
However, pursuing the argument as just stated yields 
the wrong value of the constant it was Son who 

realized that the coUinear divergence is cut off by Lan- 
dau damping at a larger value of the angle 9 than that 
at which the Meissner effect does so. Loosely speaking. 
Landau damping leads to (5 ~ (Am^)^/'^ ^ A. Son was 
then able to calculate the coe fficients of the 1/g term 
and of the logarithm in §^ |So3, Il999t ). The calcu 



lation of the constant bp was in it iated in (iHong et 



1 20001: IHsu and Schwetl. 120001: iPisarski and Rischke 
boOOallbl: IS chafcr and Wilczek'. h999S] an d com pleted in 
(jBrown ai. 2000b; Wang and RischkcJ, [IMl • Higher 
order terms are expected to be of order glogg, order g, 
and at higher order still may contain fractional powers 
and logarithms of g, see Sec. IV. A. 21 

The (A//z) DC exp(— constant/^) behavior means that 
the color superconducting gap is parametrically larger 
at /i — > oo than it would be for any four-fermion inter- 
action. Furthermore, asymptotic freedom ensures that 
l/g^nY increases logarithmically with /.t, which means 
that exp[— constant/(7(/.t)] decreases more slowly than 1/ ^ 
at large ^. We can therefore conclude that A increases 
with increasing /i at asymptotically large fi, although of 
course A/ ^ decreases. 

We conclude this subsection with a derivation of the 
correct value of the coefficient namely the constant 
in (A/fi) oc exp(— constant/g). This coefficient turns out 
to be independent of the spin-color- flavor structure M, 
and it is therefore simplest to present its derivation in 
the 2SC phase, in which there is only one gap parame- 
ter Afe = Ak^r=i, efc = £k,r=i- The leading behavior of 
the gap is completely determined by magnetic gluon ex- 
changes. We can also approximate the trace term by its 
value in the forward direction 7^(k, q) ~ 7^(k, k) = 2/3 
and set the wave function renormalization Z{qQ) = 1 (in 
the forward limit we also find 7oo(k, q) ~ '?t(k, q)). Car- 
rying out the angular integrals in the gap equation gives 



Afc = 



IStt^ 



aq — - log 



eg 2 



(41) 



Son observed that at this order we can replace the log- 
arithm by max{log(/^/efc), log(/.t/eg)}. Introducing loga- 
rithmic variables x = log[2/i/(^fc -f- ej,)] with ^j. = |fc — 
the integral equation (|¥T|) can be written as a differential 
equation 



A"(x) 



IStt- 



:A(:r), 



(42) 



with the boundary conditions A(0) = and A'(a;o) = 0. 
Here, xq ~ log(2/x/A) determines the gap on the Fermi 



surface. The solution is 

and thus 6_i = 3n^/V2. We conclude that in the weak- 
coupling limit the gap function is peaked near the Fermi 
surface, with a width that is much smaller than fi but 
much larger than A. Had we not set Z{qQ) = 1, the 
x-dependence of A{x) would be more complicated than 
the simple sinusoid in (1431), but the conclusion remains 
unchanged ( Wang and Rischki . l2002t ) . 



E. Gap and critical temperature at weak coupling 

The gap on the Fermi surface of a color superconductor 
at zero temperature can be written as 



Hg''"'e^'"' cxp 



V2gJ 



(44) 



to order g^ in the weak-coupling expansion of log(A//.i). 
We have derived the coefficient &_i in the exponent 
above, starting from a simplified version of the gap equa- 
tion p4p . with no wave function renormalization and a 
simplified gluon propagator. Upon restoring these effects, 
analysis of the gap equation yields 



5/2 



(45) 



In the following we shall define and explain the origin of 
each of the terms in this equation; we shall not present a 
complete derivation. 



• The factor g~^ and the numerical factor in Eq. ((45|) 
are due to large angle magnetic as well as electric 
gluon exchanges and are independent of the pattern 
of pairing in the color superconducting phase, i.e. 
independent of M.. 



The factor 



exp 



0.177 



(46) 



arises from the 
tion factor _Z(^ 



wave f unction reno r maliza- 
( Brown et aH l2000bl: 



IWang and Rischkg . l2002l ) and is also independent 
of M. and hence the same for all color supercon- 
ducting phases. 

• The factors that we have written as e~''e^'» are 
different in different color superconducting phases. 
The factor e~'^ is due to the angular structure of 
the gap. For the J = condensates whose gap 
equation we have derived, e~'^ = 1. Upon redoing 
the angular integrals for spin-1 con densates, we find 
that they are strongly suppressed I Schafed . l2000bl : 
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ISchmitd l2005l : ISchmitt et d\ . |2002| ). For spin-1 
pairing patterns in which quarks of the same chi- 
rahty form Cooper pairs, d = 6. A smaller suppres- 
sion occurs when quarks of opposite chirality pair, 
d = 4.5. Superpositions of these states yield values 
of d between these limits. Regardless, perturbative 
QCD predicts spin-1 gaps to be two to three orders 
of magnitude smaller than spin-0 gaps. 

• The factor depends on Ai, the color-flavor-spin 
matrix that describes the pattern of pairing in a 
particular color superconducting phase. In a phase 
in which MAi^ has two different eigenvalues Ai 
and A2, describing Tr[Pi] and Tr[7'2] quasiparticles 
respectively, we find 



c 



l(TT[Vi]Xi\ogXi 
2 



Tr[V2]X2 logAa) 



{TT[Vl]Xl+TT[r2]X2) 



(47) 



where the angular brackets denote an angular av- 
erage (trivial for J = phases). In the CFL phase, 
Ai = 1 and Tr[Pi] = 8 while A2 = 4 and Tr[7'2] = 1, 
meaning that there are 8 quasiparticles with gap A 
and 1 with gap 2 A. Evaluating Eg. (l47l) . we fin d 



2-1/3 in the CFL phase (jSchafeil . l2000at) . 



In the 2SC phase, e" 



1. Note that the ratio 



AnFT,/A2sn is also 2~-^/^ in an NJL model analysis 
( Raiagopal and Wilczekl . [2000[ ): this result depends 
on the structure of the condensates not on the na- 
ture of the interaction. From Acfl/A2sc = 2^^/'^ 
we can conclude the discussion begun in Sec. IIV.CI 
noting now that the condensation energy in the 
CFL phase is larger than that in the 2SC phase 
by a factor 3-2-2/3. 



• We can also determine the admixture of a color 
symmetric condensate in the CFL phase. In this 
case we have to use a two-parameter ansatz for 
the gap and solve two coupled gap equations. The 
color-symmetric gap parameter Ag is parametri- 
cally small compared to the color-antisymmetric 
gap A3, and the two gap equations decouple. We 
find Ag/Ag = V21og(2)g/(367r) which IS sup- 
pressed by both the coupling constant g and a large 
numerical factor. 



In evaluating (|^^ and (|l5l) . it suffices at present to 
evaluate g at the scale fi. The effect of choosing g{afi) 
with a either some purely numerical constant or some 
constant proportional to g or log A is order g, mean- 
ing that we cannot and need not determine a within 
our present calculation of log A to order g" . For a nu- 
merical estima te of the effects of a running g on A, sec 
(ISchafer j_ 2004 ]. The e ffects are not as large as envisioned 
in (jBeane et aZ.I . l200l[ ). 

We shall discuss a systematic approach to the cal- 
culation of corrections beyond 0{g^) in log(A//i) in 
Sec. IV. A. 21 There are a number of effects that have been 
considered, and were shown not to contribute at 0{g^), 



but for which the actual size of the 0{g) (or higher) cor- 
rection is not known. These include vertex corrections 
( Brown et a/.'.'2000aV the imaginary part of the gap func- 
tion IjEeng e t ai, 2006; Renter, 2006), and the modifica- 
tion of the glu on propagator due to the Meissner effect 
(lRischkel[200l . 

It is instructive to extrapolate the perturbative results 
to lower baryon densities for which the running coupling 
constant is not small. Taking fi ~ 400 — 500 MeV, and 
a strong coupling constant g ~ 3.5 (note that g ~ 3.56 
at /i = 400 MeV according to the two-loop QCD beta 
function, which of course should not be taken seriously 
at these low densities) one obtains A ~ 20 MeV. This is 
comparable to (but on the small side of) the range of typ- 
ical g aps A - (20 - 100) MeV (|Raiagopal and WilczeM . 
I2OOOII obtained using models in which the interaction be- 
tween quarks is described via a few model parameters 
whose values are chosen based upon consideration of zero- 
density physics, like the NJL models that we shall discuss 
in Sec. lVIIor numerical solut i ons o f the Dyson-S c hwin ger 
equations (jMarhauser et all . l2007t iNickel et all . [2006b). 
This qualitative agreement between two completely dif- 
ferent approaches, one based on using a model to ex- 
trapolate from ^ = to ^ = 400 - 500 MeV and the 
other based on applying a rigorous calculation that is 
valid for /i > 10^ MeV where the QCD coupling is weak 
at = 400 — 500 MeV gives us confidence that we un- 
derstand the magnitude of A, the fundamental energy 
scale that characterizes color superconductivity. Fur- 
thermore, the one nonperturbative interaction in QCD 
whose contribution to A has been evaluated reliably at 
high density, namely that due to the 't Hooft interac- 
tion induced by instantons, serves to increase A, bringing 
the high density computation into ev en better agr e ement 
with the model-based approaches (lAlford et all 1998 : 



Berges and Raia gopal'. '1999"; 'Carter and Diakonovl 1999 ; 



Rapp et all [T998i . ,200Q : .Sch^fer, .200i r 



Finally, we can use the gap equation (IMl) to extract 
the cr iti cal temperature Tr.. The resul t is ( Brown et 



2000bll3: iPisarski and Rischkd. l2000allb[ 
Schmitt et all l2002l : IWang and Rischkg . l2002t) 



Schmitl l200~ : 



Tc 
A 



(48) 



where 7 ~ 0.577 is the Euler-Mascheroni constant. This 
should be compared to the BCS result Tc/A = e"^ /tt ~ 
0.57. We observe that deviations from the BCS ratio 
occur in the case of two-gap structures and/or anisotropic 
gaps. Nevertheless, since is of order one, the critical 
temperature is always of the same order of magnitude as 
the zero-temperature gap. We see that for the 2SC phase 
Tc/A is as in BCS theory, whereas in the CFL phase this 
ratio is larger by a factor of 2^^^. It therefore turns out 
that Tc is the same in the CFL and 2SC phases. 

These estimates of Tc neglect gauge field fluctuations, 
making them valid only at asymptotic densities. We shall 
see in Scc. lV.Bj that including the gauge field fluctuations 
turns the second order phase transition that we find by 
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analyzing (j34p into a first order phase transition, and 
inereases Tc by a faetor 1 + 0(5), see Eq. ((74)) . 



F. Color and electromagnetic Meissner effect 

One of the characteristic properties of a superconduc- 
tor is the Meissner effect, the fact that an external mag- 
netic field does not penetrate into the superconductor. 
The external field is shielded by supercurrents near the 
interface between the normal phase and the supercon- 
ducting phase. The inverse penetration length defines a 
mass scale which can be viewed as an effective magnetic 
gauge boson mass. 

This effect can a lso be described as th e And erson- 
Higgs phenomenon ( Andersonl . Il963l : iHiggsl . [l963) ■ The 



difermion condensate acts as a composite Higgs field 
which breaks all or part of the gauge symmetry of the 
theory. The gauge fields acquire a mass from the Higgs 
vacuum expectation value, and the would-be Goldstone 
bosons become the longitudinal components of the gauge 
fields. A well known example in particle physics is pro- 
vided by the electrowcak sector of the standard model. 
The SU{2)l x gauge symmetry of the electroweak 

standard model is broken down to U{\)q by the expec- 
tation value of an SU (2) Higgs doublet which carries hy- 
percharge. There are three massive gauge bosons, the 
and the Z boson. The Z is a linear combination of 
the original and Y gauge bosons. The orthogonal lin- 
ear combination is the photon, which remains massless 
because the Higgs condensate is electrically neutral. 

The gauge symmetry in QCD is SU{^)c x J7(l)g. 
Different color superconducting order parameters re- 
alize different Higgs phases. The color gauge group 
may be partially or fully broken, and mixing be- 
tween diagonal gluons and photons can occur. In 
the following we shall concentrate on the 2SC and 
CFL phases and briefly mention other phases at the 
end of the section. Our starting point is the one- 



looD eauee boson Dolarization tensor (iLitim and Manuel 


2001: Rischke. 2000allbl: 


Rischke and Shovkovv, 2002 


Son and Steohanovl. 2000a 


), 



Km 



1 T 
2V 



Y^T^[f':,S{K)flS{K-P)\, (49) 



K 



where 



diag(.9 7^T, , -g 7^Tj) for a = 1 , . . . , 8 , 
diag(e7^Q,-e7^g) for a = 9. 

(50) 

Here, Q is the electric charge matrix Q ~ 
diag(2/3, — 1/3, — 1/3), and e the electromagnetic cou- 
pling constant. The polarization function can be defined 
as the second derivative of the thermodynamic potential 
with respect to an external gauge field. This quantity 
is equal to the derivative of the induced charge/current 



with respect to the gauge field. Electric charge screening 
is governed by the zero momentum limit of the Hoo com- 
ponent of the polarization tensor. The Meissner effect is 
related to a non-vanishing zero momentum limit of the 
spatial components H^ . We can define electric (Debye) 
and magnetic screening masses as 



ml,.,, = -limHS]O(0,p) 



(51a) 



ml,,,, EE - lim (5^-' - pV)n^,(0, p) . (51b) 

A calculation of the full momentum dependenc e of H,,,,., 
in th e 2SC and CFL phases can be found in ( Rischkd . 
|2001[ ) and (jMalekzadeh and Risch"k3 . l2006f) , respectively. 
One result that we shall need in Sec. IV. CI is the electric 
screening mass for gluons in the CFL phase, which is 
given by 



2 ,,2 



21 - 81og2g^/i 
36~ 



(52) 



The numerical factor 21 — 8 log 2 can be written as 15 4- 
(6 — 8 log 2), where the first term comes from diagrams in 
which the gluon couples to two octet quasiparticles and 
the second from coupling to one octet and one singlet 
quasiparticle. The log 2 factor is the log of the ratio of 
the singlet and octet gaps. 

In the following, we shall di scuss the Meissner masses. 
Results for the CFL phase (ICasalbuoni et aLl. 2001b : 



Rischkel l2000at ISon and Stephan ovl i2000al: IZarembo. 
2OOOD a nd the 2SC phase ((Casalbuoni et al, 2002b; 



Rischkel l2000d) are summarized in Table [TTl where we 
also lis t the screening masses for the single-flavor CSL 
phase ijSchmitt et aU boOSl . [200I . 

We observe that the chromomagnetic screening masses 
arc of order 17/1. This means that the screening length is 
much shorter than the coherence length ^ — 1/A, and 
color superconductivity is type I, see Sec. IV. Bl The fact 
that the screening masses are independent of the gap 
does not contradict the fact that there is no magnetic 
screening in the normal phase. Magnetic screening dis- 
appears for energies and momenta larger than the gap. 
Therefore, if the A ^ limit is taken before the limit 
p ^ then the magnetic screening vanishes, as expected. 
Of course, magnetic screening masses also vanish as the 
temperature approaches Tc. 



We also note that : 



is a 9 X 9 matrix, and the phys- 



ical masses are determined by the eigenvalues of this ma- 
trix. In the CFL phase, all magnetic gluons acquire the 
same nonzero mass, reflecting the residual SU{'i)c+L+R- 
In the 2SC phase, the Meissner masses of the gluons 1 
through 3 vanish, reflecting the symmetry breaking pat- 
tern SU{i)c SU{2)c- The unscreened gluons corre- 
spond to the generators of the unbroken SU{2)c, as they 
only see the first two colors, red and green. Cooper pairs 
are red-green singlets and so cannot screen these low mo- 
mentum gluons. 

In both 2SC and CFL phases, the off-diagonal masses 
vanish except for the eighth gluon and the photon. 
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M,78 



^ 0. The two-by-two part of the gauge 
boson mass matrices that describe the eighth gluon and 
the photon has one vanishing eigenvalue and one nonzero 
eigenvalue. The eigenvectors are characterized by a mix- 
ing angle 6*, given in the last column of Table [TTl This 
angle defines the new gauge fields, 



cos 6 y4° + sin , 



= - sin 6* A° + cos e Afj, , 



(53a) 
(53b) 



where and denote the fields for the eighth gluon 

and the photon, respectively. The A^ gauge boson feels 
a Meissner effect; it is the analogue of the massive Z- 
boson in the electroweak standard model. The Af^ gauge 
boson, on the other hand, experiences no Meissner effect 
because the diquark condensate is Q-neutral. This is the 
photon of the unbroken Abelian U{1)q gauge symme- 
try, consisting of simultaneous color and flavor (i.e. elec- 
tromagnetic) rotations. The field satisfies Maxwell's 
equations. Because g ^ e, the mixing angle is very small 
and the A^ photon contains only a small admixture of 
the original eighth gluon. 

In contrast, J = 1 color su perconductors show an elec - 
tromagnetic Meissner effect (jSchmitt et"all . |2003| . [200l . 
For example, in the CSL phase there is no mixing be- 
tween the gluons and the photon, as can be seen in the 
last row of Table UTl The photon acquires a mass since the 
electromagnetic group is spontaneously broken. Other 
candidate spin-1 phases, such as the polar, planar, or A 
phase involve mixing but also (except for a one-flavor 
system) exhibit an electromagnetic Meissner effect. This 
difference in phenomenology of spin-0 vs. spin-1 color su- 
perconductors m ay have consequences in compact stars 
(|Aguileral . [2?)07h . 



G. Chromomagnetic instability 



We have just seen in Sec. IIV.FI that color supercon- 
ductors have nonzero Meissner masses for some gluons 
and/or the photon, indicating a color or electromagnetic 
Meissner effect. However, as we discussed previously, in 
Sec. IIII.Bl if the CFL phase is stressed by a nonzero 
strange quark mass to the point that Cooper pairs break, 
the resulting gapless CFL (gCFL) phase foimd in anal- 
yses that presume a translationally in yariant condensate 
exhibi t s imaginary Meiss ner masses ( Casalbuoni et all . 
l2005bl : iFukushimal . |2005| ). This phenomenon was first 
discovered in the sim pler g apless 2SC (g2SC) phase 
(iHuang and Shovkovvl l2004allbl ) and can be understood 
in either the gCFL or g2SC context vi a a simplified anal- 
ysis i nvolving two quark species only ( Alford and Wari3 . 
I2005D that we introduced in Sec. lI.EI and shall pursue 
here. The negative Meissner mass squared implies that 
these phases are unstable toward the spontaneous gener- 
ation of currents, that break translation invariance. In 
this section we shall review the calculation of the Meiss- 
ner mass in a gapless color superconductor. 



We have seen in Sec. [TT] that the introduction of a 
nonzero strange quark mass, combined with the require- 
ment that matter be neutral and in beta equilibrium, 
serve to exert a stress on the CFL pairing that is con- 
trolled by the parameter ml/ (2^). This stress seeks to 
separate the bu and rs Fermi surfaces (and the bd and gs 
Fermi surfaces) but in the CFL phase they remain locked 
together in order to gain pairing energy cx A per Cooper 
pair. In the gCFL phase, on the other hand, there are 
unpaired bu and bd quarks in regions of momentum space 
in which the corresponding rs and gs states are empty 
— Cooper pairs have been broken yielding gapless exci- 
tations. We can describe the resulting chromomagnetic 
instability generically by picking one of these pairs, call- 
ing the quarks 1 and 2, and labelling their effective chem- 
ical potentials fii and /i2. The quasiparticle dispersion 
relations are 



(54) 



with the average chemical potential fi and the mismatch 
in chemical potentials 6fi as in Eq. (Note that the 
leading effect of the strange mass, oc ml/ {2^), is in- 
cluded in the effective chemical potential, meaning that 
we may use the massless dispersion relation of Eq. (|54p .) 
For fj,i = fi2 this yields identical dispersion relations for 
both degrees of freedom (and the same with a minus sign 
for the corresponding hole degrees of freedom which are 
omitted here). This is the usual situation of BCS su- 
perconductivity. For /ii 7^ /i2, however, one obtains two 
different quasiparticle excitations. A qualitative change 
appears at 5fi = A. For dfi > A the dispersion relations 
become gapless at the two momenta 



fc± = M ± v^p^^ 



(55) 



meaning that there are gapless quasiparticles on two 
spheres in momentum space. In the region of momentum 
space between these two spheres, the states of species 1 
are filled while those of sp ecies 2 are empty: the 1- 2 pair- 
ing has been "breached" ( Gubankova et "all. I2003D . (We 
have taken 5fi > 0.) This seems a natural way for the 
system to respond to the stress Sfj,, by reducing the num- 
ber of 2 particles relative to the number of 1 particles, 
albeit at the expense of lost pairing energy. In the larger 
gCFL context, such a response serves to alleviate the 
stress introduced by the requirements of neutrality and 
weak equilibrium. 

Gapless superconductivity ( Alford et aZ.! . [20b0al ) refers 
to the circumstance in which two species of fermions are 
paired in some regions of momentum space but in a shell 
(breach) in momentum space, bounded by two spherical 
effective Fermi surfaces, one finds unpaired fermions of 
just one of the two species. The term does not refer to 
situations in which some fermion species pair through- 
out momentum space while others do not pair at all, 
as for example in the 2SC phase. Nor does it apply to 
anisotropic superconductors in which the gap parameter 
vanishes in certain directions on the Fermi surface, as for 
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TABLE II Zero-temperature Meissner masses vtim, rotated Meissner masses ttia/, and gluon/photon mixing angle Q. The 
number a labels the gluons (a — 1,...,8) and the photon (a = 9). All masses are given in units of Nffj,^/{6TT^), where 
Nf = 3,2, 1 in the CFL, 2SC, CSL phases, respectively. We have abbreviated ?? = (21 - 81og2)/54, a = {3 + 41og2)/27, 
/3 = (6 — 41og2)/9. For the one-flavor CSL phase we denoted the quark electric charge by q. While the rotated photon in the 
CFL and 2SC phases is massless, the photon acquires a Meissner mass in the CSL phase. 



example in some single-flavor color superconductors or 
in the curCFL-A'*^ and crystalline color superconducting 
phases. The g2SC and gCFL phases are gapless super- 
conductors, in which the same quarks pair, yielding a 
nonzero order parameter A, while simultaneously featur- 
ing gapless excitations. Such phases turn out to suffer 
from the chromomagnetic instability as we now explain. 

The calculation of the Meissner mass can be done start- 
ing from Eq. (j49p . At zero temperature, in this simple 
context with two fermion species one finds 



m 



M 



1 



(56) 



where tuq is the Meissner mass obtained upon setting 
(5/i = 0, removing the stress. This expression shows that 
the Meissner mass becomes imaginary if and only if the 
spectrum is gapless. 

In essence, this is also what happens i n the CFL phase 
(|Casalbuoni et all l2005bl; iFukushimal . |2005[ ). In this 
case, of course, there are nine gauge bosons whose Meiss- 
ner masses were discussed in Sec. IIV.FI for the case of 
pairing in the absence of stress. The Meissner masses for 
pairing with mismatched Fermi surfaces are complicated 
and have to be computed numerically in general. How- 
ever, the reason for the negativity of mj^j is the same as 

in Eq. a negative term ex 5^-^ — appears 

for 5n > A. At the onset of gapless modes, Sfi = A, this 
term diverges and thus it dominates the Meissner masses 
at least for Sfi close to, but larger than. A. This is the 
story for the gluons Ai, A2, but it turns out that the 
Meissner masses for the gluons Aa, a = 4, 5, 6, 7, at first 
remain well-behaved for values of Sfi larger than A be- 
fore eventually also becoming imaginary for sufficiently 
large mismatches. The gluons ^3, Ag and the photon 
mix with each other. Two of the resulting new gauge 
bosons acquire an imaginary mass, just as the first two 
gluons. The third combination, Aq, remains massless, 
as expected from symmetry arguments. (The mixing be- 
tween these gauge bosons is a function of the mismatch 
and cannot be described by the mixing angle given in Ta- 
bleini) Although the details are clearly more complicated 
than in the simple two-species model, the conclusion re- 
mains that the chromomagnetic instability occurs if and 
only if there are gapless modes. 



This statement is not always correct, as the 
analysis of the gap less 2S C phase demonstrates 
(|Huang and Shovkovvi [2004allbl : iKirivamal l2006allbl ). In 
this phase, the gluons 1,2 and 3 are massless for arbitrary 
mismatches, reflecting the unbroken SU{2)c- One com- 
bination of the eighth gluon with the photon behaves as 
in Eq. (|56p while the other combination is massless. The 
Meissner masses for the gluons 4-7, however, are imagi- 
nary for > A as before but they are also imaginary in 
the parameter region A/y2 < V < A. Hence, the 2SC 
phase is unstable in a region where there are no gapless 
modes. Possible co nsequences of this int eresting behav- 
ior are discussed in (iGorbar aLl.l2006bl) and have been 
related to gluon condensation i Gorbar et a^.l . [2b06at) . 



We also know of an example where gapless pairing 
need not be accompanied by an instability. This is a 
two-species system where the coupling is allowed to grow 
so large that the gap becomes of the order of jl and 
even larger. In this case, a strong coupling regime has 
been identified where the gapless phase is free o f the 
chromomagnetic instabil ity ( Kitazawa et al\ . |2006[ ) . See 
( Gubankova et "oil. I2OO6D for a similar analysis in a non- 
relativistic system. The scenario in these examples can- 
not arise in QCD, since before /2 drops so low that A > /I, 
quark matter is replaced by nuclear matter. 

The chromomagnetic instability of the gCFL phase 
only demonstrates that this phase is unstable; it does not 
determine the nature of the stable phase. However, the 
nature of the instability suggests that the stable phase 
should feature currents, which must be counterpropagat- 
ing since in the ground state there can be no net current. 
Among the possible resolutions to the instability that we 
have enumerated in Sec. IIIIl two stand out by this ar- 
gument. In the meson supercurrent phase of Sec. IIII.Di 
that we shall discuss further in Sec. lV.Fi the phase of the 
CFL kaon condensate varies in space, yielding a current 
( Gerhold and Schafeil . l2006f ) . Ungapped fermion modes 
carry a counter-propagating current such that the total 
current vanishes. In the crystalline color superconduct- 
ing phases of Sec. Illl.Ci that we shall discuss further in 
Sec. I VII the diquark condensate varies in space in some 
crystalline pattern constructed as the sum of multiple 
plane waves. If the total current carried by the conden- 
sate is nonzero, it is cancelled by a counter-propagating 
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current carried by the ungappcd fcrmion modes that arc 
also found in the crystalline phases. It is important to 
note that in both these phases, the ungappcd modes 
have different Fermi surface topologies compared to that 
in the gCFL phase: they are anisotropic in momentum 
space, with unpaired fermions accommodated in one or 
many "caps" rather than in spherically symmetric shells. 
It turns out that in both these phases, the Meissner 
masses are real, meaning that these phases do not suffer 
from a chromomagnetic instabilit y. This was sh o wn in 
the meson supercurrent phase in ( Gerhold et al\ . l2007t ) 



and in the crystalline color superconducting phase in 
( Ciminale et We compare the free energies 

of these phases in Fig. [3l These two phases have to 
date been analyzed "in isolation" . It remains to be seen 
whether when they are analyzed in a sufficiently general 
framework that currents of either or both types are pos- 
sible they are distinct possibilities or different limits of 
the same more general inhomogeneous phase. 



V. EFFECTIVE THEORIES OF THE CFL PHASE 

At energies below the gap the response of supercon- 
ducting quark matter is carried by collective excitations 
of the superfluid condensate. The lightest of these excita- 
tions are Goldstone bosons associated with broken global 
symmetries. Effective theories for the Goldstone modes 
have a number of applications. They can be used to 
compute low temperature thermodynamic and transport 
properties, and to study the response to perturbations 
like nonzero quark masses and lepton chemical poten- 
tials. Other light degrees of freedom appear near special 
points in the phase diagram. Fermion modes become 
light near the CFL-gCFL transition, and fluctuations of 
the magnitude of the gap become light near Tc. 

Effective field theories can be constructed "top down" , 
by integrating out high energy degrees of freedom, or 
"bottom up" , by writing down the most general effective 
Lagrangian consistent with the symmetries of a given 
phase. In QCD at moderate or low density the micro- 
scopic theory is nonperturbative, and the top down ap- 
proach is not feasible. In this case the parameters of the 
effective Lagrangian can be estimated using dimensional 
analysis or models of QCD. If the density is very large 
then effective theories can be derived using the top down 
method. However, even in this case it is often easier to 
follow the bottom up approach, and determine the coef- 
ficients of the effective Lagrangian using matching argu- 
ments. Matching expresses the condition that low energy 
Green functions in the effective and fundamental theory 
have to agree. 

Quark matter at very high density is characterized by 
several energy scales. In the limit of massless quarks the 
most important scales are the chemical potential /i, the 
screening scale nig, and the pairing gap A. In the weak 
coupling limit we have /x ^ rUg ^ A. This hierarchy of 
scales can be exploited in order to simplify calculations 



of the properties of low energy degrees of freedom in the 
color superconducting phase. For this purpose we intro- 
duce an intermediate effective theory, the high density ef- 
fective theory (HDET), which describes quark and gluon 
degrees of freedom at energies below rUg . This theory will 
be described in Sees. EXTllYXSl Sees. EUEl are de- 
voted to effective theories of the CFL phase that allow 
us to determine the physics of the its low energy excita- 
tions. As we shall see in Sees. I VIll and I Villi these theories 
govern the phenomenology of the CFL phase even at den- 
sities not high enough for the weak coupling calculation 
of the gap parameter described in Sec. llVl to be reliable. 
In lV.GI we briefly mention effective field theories for some 
other color superconducting phases. 



A. High density effective theory 

The formalism discussed in Sec. lIVI can be extended to 
include higher orders in the coupling constants and the 
effects of nonzero quark masses. It can also be used to 
compute more complicated observables, like the disper- 
sion relations of collective excitations. In practice these 
calculations are quite difficult, because the number of 
possible gap structures quickly proliferate, and it is dif- 
ficult to estimate the relative importance of corrections 
due to the truncation of the Dyson-Schwinger equations, 
kinematic approximations, etc., a priori. 

There are two, related, strategies for addressing these 
issues: effective field theories and the renormalization 
group. Within the effective field theory approach we try 
to derive an effective Lagrangian for quasi-quarks and 
gluons near the Fermi surface, together with a power 
counting scheme that can be used to determine the mag- 
nitude of diagrams constructed from the propagators and 
interaction terms of the theory. This is the strategy that 
we will describe in Sees. IV.A.HIV.A.31 below. 

In the renormalization group approach we consider a 
general effective action for quarks and gluons at high 
baryon density, and study the evolution of the action 
as high energy degrees of freedom, energetic gluons and 
quar ks far away from the Ferni i surfa ce, are integrated 
out (Polchinsk i I1992I: IShankad . I1994D . This approach 
was applied to QCD with short range interaction s in 
(|Evans et oi.l . Il999allbl : ISchafer and Wilczekl . Il999hh . In 
this case one can show that for typical initial conditions 
the color antisymmetric, fiavor antisymmetric, J = 0, 
BCS interaction does indeed grow faster than all other 
terms, confirming in another way the arguments of Secs.|l] 
and ini that these are the channels in which the dominant 
diquark condensation occurs. In order to use the renor- 
malization group approach more quantitatively, one has 
to deal with the unscreened long range gluon exchanges, 
which is more difficult. Son studied the evolution of the 
BCS interaction usi ng the har d dense loop gluon prop- 
agator as an input (jSonl . [l99^ . The coupled evolution 
of static and dynamic screening and the BCS interaction 
has not been solved yet. A general scheme construct- 
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ing effective actions by integrati ng out hard modes was 
proposed in ( Renter et all . \206^ . 



1. Effective Lagrangian 

Consider the equation of motion of a free quark with 
a chemical potential /i 



(57) 



where 'ip± are eigenvectors of (a • p) with eigenvalue ±1 
and E± = — ^ ± p. If the quark momentum is near 
the Fermi momentum, p ^ pp ~ /i, then the solution 
1/'+ describes a low energy excitation ~ 0, whereas 
E_ —2^ corresponds to a high energy excitation. In 
order to construct an effective field theory based on this 
observation we de fine low and high-energy components 
of the quark field (|Hond . l2000al lbll 



1 ± a • irp 



(58) 



where vj? is the Fermi velocity and = (l,Vi?). The 
prefactor removes the rapid phase variation common to 
all fermions in some patch on the Fermi surface specified 
hyirp- We can insert the decomposition Eq. (|58p into the 
QCD Lagrangian and integrate out the ijj- field as well 
as hard gluon exchanges. This generates an expansion of 
the QCD Lagrangian in powers of 1/pp- At tree level, 
integrating out the fields is equivalent to solving their 
equation of motion 



1 /. 
2pp 



(59) 



(60) 



where 7|| = wp[-vp -7), ^ T^Ty ^'^'^ is the quark 
mass matrix. At 0{l/pp) the effective Lagrangian for 

tP+ is 

C = ^p\j^{iv'^D^)r^+.L - ^^'l,L [iPi.? + MAP] 
+ {L'^ R,M ^ AP) + .... 

The low energy ex pansion was studied in more detail 
m (jSchafed . l2003a[ ). There are a number of physical 
effects that have to be included in order to obtain a 
well-defined expansion. First, four-fermion operators 
have to be included. These operators naturally appear 
at 0(l/p|n) but their effects are enhanced by the large 
density of states N ^ pj^ on the Fermi surface. The 
most important of these operators is the BCS interac- 
tion [-i/;(v)-0(-v)][?/'t(v^)V;t(_v^)] The coefficient of the 
BCS operator was determined in (ISchafej . l2003ah . 

Because of the large density of states it is also nec- 
essary to resum quark loop insertions in gluon n-point 
functions. There is a simple generating functional for 
these effect s, known as the hard dense l oop (HDL) effec- 
tive action ( Braaten and Pisars 5 119921 



HDL 



— V r« 

9 



(61) 



This is a gauge invariant, but non-local, effective La- 
grangian. Expanding £hdl in powers of the gauge field 
produces 2,3, .. . gluon vertices. The quadratic term de- 
scribes dielectric screening of electric modes and Landau 
damping of magnetic modes. Higher order terms con- 
tain corrections to the gluon self interaction in a dense 
medium. 



2. Non-Fermi liquid effects and the gap equation 

In this section we shall analyz e the low energy expan- 
sion i n the regime A < ko < nig I Schafer and Schwenzeii 
I2OO6D . This energy range gives the dominant contribu- 



tion to the pairing gap and other low energy constants 
in the superconducting phase. Since electric fields are 
screened the interaction is dominated by the exchange of 
magnetic gluons. The transverse gauge boson propagator 
is 



k^ 



(62) 



(37r/4)mq and we have assumed that |fco| < 



where 

k. We observe that the propagator becomes large in the 
regime |fco| ^ k^/ni^. If the energy is small, |fco| "ig, 
then the typical energy is much smaller than the typical 
momentum. 



(m2jfco|)^/'»|fco| 



(63) 



This implies that the gluon is very far off its energy shell 
and not a propagating state. We can compute loop dia- 
grams containing quarks and transverse gluons by picking 
up the pole in the quark propagator, and then integrate 
over the cut in the gluon propagator using the kinematics 
dictated by Eq. ([551) . In order for a quark to absorb the 
large momentum carried by a gluon and stay close to the 
Fermi surface the gluon momentum has to be transverse 
to the momentum of the quark. This means that the term 
k'j_/{2^) in the quark propagator is relevant and has to be 
kept at leading order. Eq. ([55]) shows that k\/[2y) ^ ko 
as fco —> 0. This means that the pole of the quark prop- 
agator is governed by the condition fc|| ~ k']_/{2y). We 
find 



f/^^?/%'\ k,^g'/^^,^/^kl'\ 



(64) 



In this regime propagators and vertices can be simplified 
even further. The quark and gluon propagators are 



fcll - 



iesgn(fco) 



(65) 



(66) 



and the quark gluon vertex is gviiX^ /2). Higher order 
corrections can be found by expanding the quark and 
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gluon propagators as well as the HDL vertices in powers 
of the small parameter e = (fco/?n). 

The regime characterized by Eq. (|64[) is completely 
pertm'bative, i.e. graphs with ext ra loops are always sup- 
press ed by extra powers of e^/^ ( Schafer and Schwenzeii 
I2OO6D . The power of e can be found by using the fact that 



loop integrals scale as (^0^11^1^) ~ ^o^"^' fermion propa- 
gators scale as l/fc|| ^ l/fep^'^, gluon propagators scale 

as ~ l/k^'^ ^ and the quark-gluon vertex scales as 

a constant. Quark matter in the regime A < /cq < m is 
a non-Fermi liquid. The excitations are quasi-particles 
with the quantum numbers of quarks, but Green func- 
tions scale with fractional powers and lo garithm s of the 
energy and the coupling constant (iGerhold and Rebhanl . 
I2OO5HIPD et'ai\ . l2004l : ISchafer and Schwenzeii . l2004b^ . 

The corrections to Fermi liquid theory do not upset 
the logic that underlies the argument that leads to the 
BCS instability. For quark pairs with back-to-back mo- 
menta the basic one gluon exchange interaction has to be 
summed to al l orders, but all other i n teract ions remain 
perturbative ( Schafer and Schwenzeii . |2006[ ). The gap 
equation that sums the leading order transverse gluon 
exchange in the color-anti-symmetric channel is 



A(po) 



.25^ 



'2^ 



dkl 
(2^)2 

A(fco) 



kl 



fc2 



A(fco) 



■Po) 
(67) 



This equation is exactly equivalent to Eq. (|4ip. In par- 
ticular, all the kinematic approximations that were used 
to derive Eq. (|4T|) , like the low energy approximation to 
the HDL self energies and the forward approximation to 
the Dirac traces, are built into the effective field theory 
vertices and propagators. The effective theory can now 
be used to study corrections to the leading order result. 
Higher order corrections to the propagators and vertices 
of the effective theory modify the kernel of the integral 
equation Eq. (j67|l . The resulting correction to the gap 
function can be computed pcrturbatively, without hav- 
ing to solve the integral equation again, using a method 
that is similar to Rayleigh-Schrod i nger p erturbation the- 
ory (|Brown et all . l2000al; ISchafeil . l2003al) . 

The coefficients bo and bo introduced in Sec. lIV.Dl can 
be determined by matc hing the four- fe rmion operators 
in the effective theory ( Schaferl . l2003a[ l. The bo term 
also receives contributions from the fermion wave func- 
tion renormalization Z ^ log(fco). All other terms give 
corrections beyond 0{g'^) in log(Ao//^). Vertex correc- 
tions scale as F ~ p^^^ and are suppressed compared to 
the fermion wave function renormalization. The analo- 
gous statement in the case of phonon-induced electronic 
superconductors is known as Migdal's theorem. Gluon 
self energy insertions beyond the ko/k± term included 
in the leading order propagator are also suppressed by 
fractional powers of the coupling and the gluon energy. 



R R 



R L 
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FIG. 5 Mass terms in the high density effective theory. The 
first diagram shows a 0{MM^) term that arises from inte- 
grating out the field in the QCD Lagrangian. The second 
diagram shows a 0{M^) four- fermion operator which arises 
from integrating out and hard gluon exchanges. 



3. Mass terms 

A systematic determination of mass corrections to the 
high density effective theory is needed for calculations 
of the Goldstone boson masses in the CFL phase, and in 
order to understand the response of the CFL ground state 
to nonzero quark masses. Mass terms affect both the 
quark propagator and the quark-quark interaction. From 
Eq. ([5^ and (pU)) we see that integrating out the tp- field 
gives a correction to the energy of the 4'+ field of the form 
MM^ /{2pp). This term can be viewed as an effective 
flavor dependent chemical potential. We also note that 
this term is just the first in a tower of operators that arise 
from expanding out the energy of a free massive quark, 
E = [p^ + m^)^/^, for momenta near the Fermi surface. 
Higher order terms correspond to additional corrections 
to the chemical potential, the Fermi velocity, and to non- 
linear terms in the dispersion relation. 

There are no mass corrections to the quark-gluon ver- 
tex at 0{\/p\r). There are, however, mass corrections 
to the quark-quark interaction. In connection with color 
superconductivity we are mainly interested in the BCS 
interac tion. The diagram shown in Fig. Ogives ( Schafeij . 
l2002bD 



C = 



32p% 



[{\n\f'\MU{M),i] + {L^R,M^ Aft) . 



This is the leading interaction that couples the gap equa- 
tions for left and right handed fermions. We shall also 
see that the mass correction to the BCS interaction gives 
the leading contribution to the mass shift in the conden- 
sation energy, and the masses of the Goldstone bosons. 



B. Ginzburg-Landau theory 

At zero temperature fluctuations of the superconduct- 
ing state are dominated by fluctuations of the phase of 
the order parameter. Near the critical temperature the 
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gap becomes small and fluctuations of the magnitude of 
the gap are important, too. This regime can be described 
using the Ginzburg-Landau theory. Ginzburg and Lan- 
dau argued that in the vicinity of a second order phase 
transition the thermodynamic potential of the system can 
be expanded in powers of the order parameter and its 
derivatives. This method was used very successfully in 
the study of supcrfluid phases of ^He. 

The Ginzburg-Landau a pproach was firs t appl ied to 
color superconductivity in (iBailin and Lovd ll979l) . The 
problem was revisited by ( lida and Bavml . l200ll ) , who 
included the effects of unscreened gluon exchanges and 
charge neutrality. Consider the s-wave color anti-triplet 
condensate in QCD with three massless flavors. The or- 
der parameter can be written as 



(69) 



where (p^ is a matrix in color- flavor space. Note that here 
we have included the energy gap into (p^, in contrast to 
Eq. ([24|) . where (f>^ is dimensionless. We have fixed the 
orientations of left and right-handed condensates. Fluc- 
tuations in the relative color-flavor orientation of the left- 
and right-handed fermions correspond to the Goldstone 
modes related to chiral symmetry breaking, and will be 
considered in Sec. IV. CI Therefore, the ansatz ([55]) im- 
plies the assumption that chiral fluctuations near Tc are 
small compared to non-chiral gap fluctuations and fluctu- 
ations of the gauge field. The thermodynamic potential 
can be expanded as 

n = no + aTT{q}U) +l3i[Tr{(bU)]^ (70) 
+ (32 Tr ( [0^0] ^) + K Tr (V(/)V0^) + . . . 

The coefficients a,f3i,K can be treated as unknown pa- 
rameters, or determined in Q CD at weak coupling. The 
weak coupling QCD result is ( lida and Bavml . boOll ) 



a = m^-J^, 13, =(32 = N, (71) 



T 



7C(3) 
8^2 T2 



N 



(72) 



gauge field, and the transition is second order. In a type 
I superconductor the situation is reversed, and fluctua- 
tions of the gauge fie ld drive the transition first order 
(|Halperin et aZ.I . Il974[ ). 



In the weak coupling hmit, ^ ~ 1/A X ^ ^/{giA 
and color superconductivity is strongly type I. The role of 
gaug e field fluctuation s was studied in ((Bailin and Love , 
19841: iGiannakis et all |2004 iGiannakis and Renl . 12003 : 



Matsuura et al\ . l2004l : iNoronha et al\ . l2006l ). The con 



tribution to the thermodynamic potential is 



^lfi=^T 



(2^ 



log 1 



fc2 



fc2 



(73) 

where mA{k) is the gauge field screening mass. In QCD 
the momentum dependence of mA cannot be neglected. 
The contribution of the fluctuations O^; induces a cubic 
term oc ((>^ in the thermodynamic potential which drives 
the transition first order. The first order transition occurs 
at a critical temperature T* ([Giannakis et aL . i200j) 



T* - T, 



12V2^' 



(74) 



where Tc is the critical temperature of the second order 
transition obtained upon neglecting the cubic term. Al- 
though the result ([71]) cannot be trusted quantitatively 
at accessible densities, say fj, ~ 400 McV where g ~ 3.6, it 
does make it clear that the phase transition between the 
CFL (or 2SC) phase and the quar k-gluon plasma \yill be 
strongly first order. Noronha et al. ( Noronha et aZ.I . [20061 ) 
argue that Eq. (|74|) gives the complete 0(q) correction to 
the c ritical temperature (see, however, ( Matsuura et al\ . 
|2004| )V This implies that the transition to the color su- 
perconducting phase will occur at a critical temperature 
that is significantly elevated relative to the BCS estimate 
Tc = 0.57A that we obtained in Sec. IIV.EI The effects 
of gluon fluctuations are much more important here than 
those of photon fluctuations in a conventional type I su- 
perconductor. 



where N = ^"^ /{2tt'^) is the density of states on the Fermi 
surface. This result agrees with the BCS result. Using 
Eq. (|7T|) we can verify that the ground state is in the CFL 



phase 



We can also study many other issues. 



like the gluon screening lengths, the structure of vor- 
tices, the effects of electric and co lor neutrality, and the 
effects of nonze r o qua rk masses (|lida and Bavml 120021: 
llida et aLl . [2003 . 12005D . 

From the study of electronic superconductors, it is 
known that the nature of the finite temperature phase 
transition depends on the ratio k = A/^ of the screen- 
ing length A and the correlation length ^. If k > 1/^/2 
the superconductor is type II, fluctuations of the order 
parameter are more important than fluctuations of the 



C. Goldstone bosons in the CFL phase 

1. Effective Lagrangian 

In the CFL phase the pattern of chiral symmetry 
breaking is identical to the one at T = /i = 0. This im- 
plies that the effective Lagrangian has the same structure 
as chiral perturbation theory. The main difference is that 
Lorentz-invariance is broken and only rotational invari- 
ance is a good symmetry. The effective Lagrangian for 
the Goldstone modes is given by ( Casalbuoni and Gattol . 
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In terms of the new fields the Lagrangian takes the form 



£eff = ^Tr [9oS5oS^ - ^''5,E9,Et] 



(75) 



Tr 



+ [BTr(A/S^) + h.c] 

+ [^iTr(A/S^)Tr(AfS^) + A2Tt{MY) MY)) 

+A3Tr(ME^)Tr(il/^I]) + h.c] + .... 

Here E = eyiY){i(jf'\°' / f-n) is the chiral field, /^r is the 
pion decay constant and M is a complex mass matrix. 
The fields 0" describe the octet of Goldstone bosons 
(tt*, 7r°, _ftr*, 77). These Goldstone bosons are 

an octet under the unbroken SU{i)c+L+R symmetry of 
the CFL phase and their Q-charges under the unbroken 
gauge symmetry of the CFL phase are ±1 and as in- 
dicated by the superscripts, meaning that they have the 
same Q-charges as the Q-charges of the vacuum pseu- 
doscalar mesons. The chiral field and the mass matrix 
transform as E — > and M — *■ LMB) under chi- 

ral transformations {L,R) e SU{3)l x SU{3)b. For the 
present, we have suppressed the singlet fields associated 
with the breaking of the exact U{1)b and approximate 
U{1)a symmetries. We will give the effective Lagrangian 
for the massless Goldstone boson associated with super- 
fiuidity (i.e. from U{1)b breaking) below. 

The form of the effective Lagrangian follows from the 
symmetries of the CFL phase. It is nevertheless useful to 
understand how this Lagrangian arises upon integrating 
out high energy degrees of freedom. We start from the 
high density effective Lagrangian in the presence of a 
CFL gap term 
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+ - {Tr (XtV^i^t^i) - [Tr {X^l)]^ + h.c.} 

+ {L ^ R,X ^Y) . 

Here, tpL,B. are left and right-handed quark fields which 
transform as -0^ — > LiPlU'^ and tpR — > RiPrU'^ under 
chiral transformations {L,R) G SU{3)l x SU{3)r and 
color transformations U £ SU{3)c. We have suppressed 
the spinor indices and defined = ip^C^^ilj^ , where 
C is the charge conjugation matrix. The traces run over 
color or flavor indices and X, Y are fields that transform 
as X ^ LXU'^ and Y RYU"^. We wiU assume that 
the vacuum expectation value is {X) = (Y) = 1. This 
corresponds to the CFL gap term Aip^ipj ea/SA^^-'^- X, Y 
parametrize fluctuations around the CFL ground state. 
Note that fluctuations of the type X = Y correspond to 
the field (f> introduced in the previous section. 

For simplicity we have assumed that the gap term is 
completely anti-symmetric in flavor. We will derive the 
effective Lagrangian in the chiral limit M = M''' = 
and study mass terms later. We can redefine the fermion 
fields according to 



xUw^i9^)XL 



(78) 



— i Tr 
A 



f {Tr(xLXL)-[Tr(xL)]'} 

+ {L<-^ R,X ^Y). 

At energies below the gap we can integrate out 
the fermions. The fermion determinant generates 
a kinetic term for the chiral fields X and Y 
(|Casalbuoni and Gattol . Il999l ) 



f2 

C = - yTr [{X^DqX)^ 



liX^D^Xf] + {X 



Y) 
(79) 

For simplicity we have ignored the flavor singlet compo- 
nents of X and Y . 

The low energy constants and were calcu- 
lated by matching th e effective theory to weak co upling 
QCD calculations i n dSon and Stephanov . 2000al |b ), see 
also (iBeane all. 200C : iBedaaue and Schafeil . l2002t 
iRho et aLl . l2000l : IZarei^ . l2000[ ). The results are 



21 - 81og2 fi'^ 



18 



27r^ 



(80) 



The simplest way to derive these results, given the re- 
sults that we have already reviewed in Sec. IIV.FI is to 
recall that the gluon fleld acquires a magnetic mass due 
to the Higgs mechanism and an electric mass due to De- 
bye screening, and then to notice that Eq. (|79p shows that 
the electric mass is rn^ = g^f^, while the magnetic mass 
is m\j — v'^m'jj. This means that and are deter- 
mined by the Debye and Meissner masses for the gluons 
in the CFL phase that we have presented in Sec. lIV.Fi 
see Eq. ^ and Table |TT1 

Since the gluon is heavy, it can also be integrated out. 
Using Eq. ^ we get 



A. 



Y^d^Y) 



(81) 
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(77) 



This result can be substituted back into the effective La- 
grangian. The result is 

C,s^ ^TT[dom^^ -vld,m^^] , (82) 

where the Goldstone boson field is given by E = XY^ . 
This shows that the light degrees of freedom correspond 
to fluctuations of the color-flavor orientation of the left- 
handed CFL condensate relative to the right-handed one, 
as expected since these are the fluctuations associated 
with the spontaneously broken global symmetry. 



2. U{1)b modes and superfluid hydrodynamics 

Finally, we quickly summarize the effective theory for 
the U{1)b Goldstone mode. At order 0{{dip)^) we get 
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(|Son and Stephanovl . [2000al lbl) 



(83) 



where the low energy constants / and v are given by 



.2 „ 



(84) 



The field ip transforms d& tp — > Lp + a under U{X)b 
transformation of the quark fields ^ exp(iQ!)?/'- Be- 
cause U{1)b is an Abelian symmetry, the two-derivative 
terms do not contain any Goldstone boson self interac- 
tions. These terms are needed in order to compute trans- 
port properties of the CFL phase. Son noticed that self- 
interactions are constrained by Lorcntz invar iance (of th e 
microscopic theory) and U{1)b invariance (fS^ . I2002D . 
The analogous argument for non-r elativistic superfluids 
is described in (jGreiter et al 1 119891) . To leading order in 
g the effective theory of the U{1)b Goldstone boson can 
be written as 



(85) 



where the omitted terms are of the form d^ip'' with i > k. 
Expanding Eq. (j85p to second order in derivatives repro- 
duces Eq. ([83)1 . In addition to that, Eq. ([85)1 contains the 
leading three and four boson interactions. Using micro- 
scopic models one can obtain more detailed information 
on the properties of collective modes. A calculation of 
the spectral properties of the ip mo de in an NJL model at 
T = and T ^ can be found in ( Fukushima and 
I2005D . 

The spontaneous breaking of U {1)b is related to super- 
fluidity, and the U{1)b effecti ve theory ca n be interpreted 
as superfluid hydrodynamics (jSonl . [2003 l. We can define 
the fluid velocity as 



Da(p, 

Mo 



(86) 



where DqI/j = da(p+{n, 0) and jiQ = (DaipD^ipY/^ . Note 
that this definition ensures that the flow is irrotational, 
V X u = 0. The identification ((86|) is motivated by the 
fact that the equation of motion for the U{1) field ip can 
be written as a continuity equation 



d"{noUa) = 0, 



(87) 



where uq ~ S^q/tt^ is the superfluid number density. At 
T = the superfluid density is equal to the total density 
of the system, n = dP j dii\^=^„. The energy-momentum 
tensor has the ideal fluid form 

Tap = {(- + P)UaUp - Pgap, (88) 

and the conservation law d^Tap = corresponds to the 
relativistic Euler equation of ideal fluid dynamics. We 
conclude that the effective theory for the U{1)b Gold- 
stone mode accounts for the defining characteristics of 
a superfluid: irrotational, non-dissipative hydrodynamic 
flow. 
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FIG. 6 Contribution of the 0{M') BCS four-fermion opera- 
tor to the condensation energy in tlie CFL phase. The open 
squares correspond to insertions of the anomalous self energy 
A. 



3. Mass terms 

The structure of the mass terms in Eq. ((75|) is com- 
pletely determined by chiral symmetry. The coefficients 
B, Ai can be determined by repeating the steps discussed 
in the previous section, but keeping the mass terms in the 
high density effective theory. In practice it is somewhat 
easier to compute the coefficients of the chiral Lagrangian 
using matching arguments. For example, wc noticed that 
the easiest way to determine /^r is to compute the gluon 
screening mass in the microscopic theory. 

In Sec. IV. A. 31 wc showed that Xl = MA/V(2pi^) and 
X]i = M^M/{2pf) act as effective chemical potentials 
for left and right-handed fermions, respectively. For- 
mally, the effective Lagrangian has an SU{?>)l x SU{?>)ii 
gauge symmetry under which X^ fj transform as the tem- 
poral components of non- Abelian gauge fields. We can 
implement this approximate gauge symmetry in the CFL 
chiral theor y by promoting time derivat ives to covariant 
derivatives ( Bedaaue and Schafed . 1200^ . 



V 2PF J 



The mass dependent terms in the quark-quark interac- 
tion contribute to the gap and to the condensation en- 
ergy. In the chiral theory the shift in the condensation 
energy due to the quark masses is 



£ = -BTt{M) - Al [Tr{M)f - AsTr (M^) 
- y43Tr(M)Tr(Mt) + h.c. + .... 



(90) 



The contribution to the condensation energy from the 
mass correction to the BCS interaction is shown in Fig. [6] 
The diagram is proportional to the square of the conden- 
sate 



9 



(91) 



with the dependence on the mass matrix AI arising from 
the contraction of the BCS interaction with the CFL con- 
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densate. Wc get 



(92) 



where = A°/2. Wc note that the four-fermion op- 
erator is proportional to and the exp licit dep e ndence 
of the diagram on q cance ls. We find ( Schafed . l2002ht 
ISon and Stephanovl . [2000al lbl;i 



the effective theory is valid but the values of coefficients 
in it may not take on the values obtained by matching to 
high density calculations. 

In perturbation theory the coefficient B of the Tr(A/S) 
term is zero. B receives non-perturbative contributions 
from instantons. Instantons are semi-classical gauge con- 
figurations in the Euclidean time fun ctional integral that 
induce a fermion vertex of the form ( 't 



(97) 



£ = - 



3A2 

47r2 



Tr[i\/]) -Tr 



(93) 

This result can be matched against Eq. ((501) . We find 

B = and 



3A2 
4^ 



The result Ai ~ —A2 reflects the fact that the CFL order 
parameter is anti-symmetric in flavor (pure 3) to leading 
order in g. Using Eqs. (|89p and ([M)) we can compute the 
energies of the flavored Goldstone bosons 



Ek± 



where det/ denotes a determinant in flavor space. The 
't Hooft vertex ip7|) can be written as the product of the 
CFL condensate and its conjugate times (iptp), meaning 
that in the CFL phase Eq. (^T)) induces a nonzero quark 
condensate {ipip), as well as Goldstone boson masses 
™GB ^ "^('0V')//^- The instanton has gauge fleld 
A, A3 =0. (94) Af^ l/g, so its action is S" = %Tr^ I g^ ■ The effective 

coupling G is proportional to cxp(— S*) ~ cxp(— Stt^/^^), 
where g is the running coupling constant at a scale set 
by the instanton size p. 

In dense quark matter perturbativc gauge field screen- 
ing suppresses instantons of size p > l/^i, an d the ef- 
fective coupling G can be computed reliably (ISchafeil . 
l2002aD . Combined with the weak coupling result for 
{iptp), see Eq. (|9T|) . we get 
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, (95) 
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(98) 



where 
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(96) 



The mass matrix for the remaining neutral Goldstone 
bosons, which mix, can be found in ( Beane et al\ . l2000t 
ISon and Stephanovl . [2000allbl V We observe that the 0{m) 
terms lead to an inverted mass spectrum with the kaons 
being lighter than the pions. This can be understood 
from the microscopic derivation of the chiral Lagrangian. 
The Goldstone boson field is S = XF^, and a mode 
with the quantum number of the pion is given by 7r+ ~ 
^abc^ade^jb^^c^^^^dge^y rpj^^ structurc of the field opera- 
tors suggests that the mass is controlled by (mu+md)ms. 
By the same argument the mass of the is governed 
by {rriu + ms)mii^ and rriK < "m^,. We also note that 
the 0{m?') terms split the energies of different charge 
states. This can be understood from the fact that these 
terms act as an effective chemical potential for fiavor. 
Explicit calculations in an NJL model reproduce in 
(I5D ]) and the results (1551) . albeit with a d ifferent value of 
A ([Kleinhaus et fl/.l . i2007i:iRuggierill2007[ ). This serves as 
a reminder that in the CFL phase at moderate densities, 



where c = 0.155 and Kqcd is the QCD scale factor. 
In terms of B, (ipip) = —2B and the instanton contri- 
bution to the inass is 5m|-o — B{md + ms)/(2/^) 
(|Manuel and Tvtgatl . [2OO0I ). In the weak coupling limit, 
p ^ Aqcd, the instanton contribution is very small. 
However, because of the strong dependence on Aqcd 
the numerical value of B is quite uncertain. Using phe- 
nome nologica l constr aints on the instanton size distribu- 
tion (jSchafed . [2002al) concluded that the instanton con- 
tribution to the kaon mass a,t p = 500 MeV is of order 
10 MeV. 

Finally, we summarize the structure of the chiral ex- 
pansion in the CFL phase. Ignoring non-perturbative 
effects the effective Lagrangian has the form 

(99) 

Higher order vertices are suppressed by 9S/A whereas 
Goldstone boson loops are suppressed by powers of 
3I]/(47r/^). Since the pion decay constant scales as 
/tt ~ /i the effects of Goldstone boson loops can be ne- 
glected relative to higher order contact interactions. This 
is different from chiral perturbation theory at zero baryon 
density. We also note that the quark mass expansion 
contains two parameters, m? j p^ and m'^/(pA). Since 
A <C /i the chiral expansion breaks down if ~ pA. 



34 



m 




FIG. 7 (Color online) Phase structure of CFL matter as a 
function of the light quark mas s m a nd the strange quark 
mass iris, from (jKrvievski et all |2005| ). CFL denotes pure 
CFL matter, while K'^ and rj denote CFL phases with K° 
and/or r/ condensation. Solid lines are first order transitions, 
dashed lines are second order. Instanton effects have been 
neglected. 



This is the same scale at which BCS calculations find a 
transition from the CFL phase to a less symmetric state. 
We also note that the result for the Goldstone boson en- 
ergies given in Eq. ((95)) contains terms of 0{m? / ^'^) and 
0([m2/(^A)]2), but neglects corrections oi 0{[m'^ / 



The efi'ective Lagrangians ([75]) and ([55)1 describe the 
physics of the low momentum pseudo-Goldstone and 
Goldstone bosons of the CFL phase at any density. We 
have described the weak coupling computation of the co- 
efficients /tt, Wttj B, f and v as well as the 
/icff's in ([M)) . With the exception of B, all these results 
are expressed simply in terms of A, and the quark 
masses, with g not appearing anywhere. This suggests 
that the range of validity of these results, when viewed 
as a function of A, is bigger than the range of valid- 
ity of the weak coupling calculations on which they are 
based. As we decrease the density down from the very 
large densities at which the weak coupling calculation of 
A is under control, there is no indication that the rela- 
tions between the effective theory coefficients and A and 
fi that we have derived in this section break down. The 
only sense in which we lose control of our understanding 
of the CFL phase is that we must treat A as a parameter, 
in terms of which all the other effective theory coefficients 
are known. Since B introduces t/(l)A-brcaking physics 
that is not present at weak coupling and that does not 
enter the effective theory through any other coupling, it 
is not well constrained. 



D. Kaon condensation 

If the effective chemical potential in Eq. becomes 
larger than the corresponding mass term in Eq. (|95p . then 
the energy of a Goldstone boson can become negative. In 
the physically relevant case Wg ^ uid this applies 

in particular to the K'^ and the K~^. When the Gold- 
stone boson energy becomes negative the CFL ground 



state is reorganized and a Goldstone boson condensate is 
formed. The physical reason is that a nonzero rris dis- 
favors strange quarks relative to non-strange quarks. In 
normal quark matter the system responds to this stress 
by turning s quarks into (mostly) d quarks. In CFL mat- 
ter this is difficult, since all quarks are gapped. Instead, 
the system can respond by populating mesons that con- 
tain d quarks and s holes. 

The ground state can be determined from the effective 
potential 



•I 7T r 

4 



(100) 



Ai{[Tr(il/S]t)]'-Tr[(A/S]t)2]}, 



where Xl = il/Mt/(2pF), Xr = MHI/{2pf) and 
M ~ diag(TO„, m(i, m^) = diag(m, to, TOj). Here we only 
discuss the T = case. For nonzero temperature effects, 
in particular the calculat ion of the crit i cal tem perature of 
kaon condensation, see (|Alford et all l2007b[) . The first 
term on the right-hand side of Eq. (|100p contains the 
effective chemical potential 



2pF 



(101) 



and favors states with a deficit of strange quarks. The 
second term favors the neutral ground state S = 1. 
The lightest excitation with positive strangeness is the 
meson. We consider the ansatz S ~ cxp(iQ;A4) which 
allows the order parameter to rotate in the direction. 
The vacuum energy is 



V{a) = fl 



2pF 



+ TO^o(l — COS a) 



(102) 



where 



(4Ai//2)TO„(TOd + TO,) + B{md + 



TOs)/(2/^). Minimizing the vacuum energy we obtain 



cos(a) 



1 < TO^O 

2 

-4- pis > TOa'O 



(103) 



We conclude that there is a second order phase tran- 
sition to a kaon condensed state at fig = nij^o. The 
strange quark mass breaks the SU{3) flavor symme- 
try to SU{2)j X U{1)y ■ In the kaon condensed phase 
this symmetry is spontaneously broken to U{1)q. If 
niu = md , isospin is an exact sy mmetry and there are two 
exact Goldstone mo des (M iranskv and Shovkovvl l2002t 
ISchafer et all . l200ll) with zero energy gap, the and 
the K'^ . Isospin breaking leads to a small energy gap for 
the K+. 

Using the pcrturbative result for Ai, and neglecting in- 
stanton effects by setting i? = 0, we can get an estimate 
of the critical strange quark mass. The critical strange 
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quark mass scales as rrij'^/^^/^. Taking ~ 500 McV, 
A = 50 MeV, m„ = 4 MeV and = 7 MeV, 
we find m"'* ~ 68 McV, a result that corresponds to 
TO^o* = 5 MeV. If instanton contributions increase m/^o 
by 10 MeV, this would increase m™'* to 103 MeV, corre- 
sponding to the onset of kaon condensation depicted in 
Fig.H 

The difference in condensation energy between the 
CFL phase and the kaon condensed state is not neces- 
sarily small. In the limit /z^ ^ A we have sin a ^ 1 
and V{a) - /^AV2. Since fl is of order mVIStt^) 
this is an 0(1) correction to the pairing energy in 
the CFL phase. Microscopic NJL model calculations 
of the condensation e nergy in the kaon condensed 



phase can be found in (iBuballal. 2005b Forbei 



iKleinhaus et all l2007t iRuggieril. l2007l: IWarringa[|2006 
see also (jEbert and Klimenkol . j2007t lEbert et all 12007 



The CFL phase also contains a very light flavor neu- 
tral mode which can potentially become unstable. This 
mode is a linear combination of the 77 and 77' and its 
mass is proportional to muTrifi. Because this mode has 
zero strangeness it is not affected by the /ig term in 
the effective potential. However, since m„,TO(j ^ nig 
this state is sensitive to perturbative asTnl corrections 



( Krvievski et aLl . l2005l ). The resulting phase diagram is 
shown in Fig. [71 The precise value of the tetra-critical 
point (m*,m*) depends sensitively on the value of the 
coupling constant. At very high density m* is extremely 
small, but at moderate density m* can become as large 
as 5 MeV, comparable to the physical values of the up 
and down quark mass. 



E. Fermions in the CFL phase 

A single quark excitation with energy close to A is 
long-lived and interacts only weakly with the Goldstone 
modes in the CFL phase. This means that it is pos- 
sible to include quarks in the chiral Lagrangian. This 
Lagrangian not only controls the interaction of quarks 
with pions and kaons, but it also constrains the depen- 
dence of the gap in the fermionic quasiparticle spectrum 
on the quark masses. This is of interest in connection 
with the existence an d stability of the gapless CFL phase 
(lAlford et a^l2004bD . as we have discussed in Sees. II. El 
HLBl andimil 

The effective Lag rangian for fe rmions in the CFL phase 
is (iKrvievski andSchafer, . i2005i : iKrvievski and Yamadal . 
I2005h 



C = Tr(Ar^w^D^iV) -OTr(A^1'n'^75{^^,iV}(jl04) 
-FTr(7VV75[^^,7V]) 

{TY{NLNL)~[Tr{NL)Y 
-{L^ R) + h.c. 
Nl.r are left and right handed baryon fields in the adjoint 
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FIG. 8 (Color online) This figure shows the fermion spectrum 
in the CFL phase. For rus = there are eight fermions with 
gap A (set to 25 MeV as in Fig. [3} and one fermion with gap 
2 A (not shown). As discussed in Sec. IIIII the octet quasi- 
particles have the SC/(3) and '7(1)q quantum numbers of the 
octet baryons. Without kaon condensation gapless fermion 
modes appear at /is = A (dashed lines). With kaon conden- 
sation gapless modes appear at \is = 4A/3. (Note that the 
scale on the horizontal axis is 2/is.) 



representation of flavor 5/7(3). The baryo n fields origi- 
nate f r om quark-hadron compleme ntarity ( Afford et al\ . 
Il999al: ISchafer and Wilc^ . Il999a^ . We can think of iV 
as describing a quark which is surrounded by a diquark 
cloud, Nj^ ~ <lh{flL<lL)- The covariant derivative of the 



nucleon field is given by D^,N = ^^,N + i[V^,Af]. The 



vector and axial-vector currents are 



(105) 

where ^ is defined by = S. It follows that ^ transforms 
as C ^ L^C/^ = U^R^ with U € SU{3)v- The fermion 
field transforms as iV ^ UNU^. For pure SU{3) fiavor 
transformations L = R = V we have U = V . F and D 
are low energy constants that determine the baryon axial 
coup ling. In QCD at weak coupl ing, we find D = F = 
1/2 (|Krvievski and Schafeil . [2005h . 

The effective theory given in Eq. (|104p can be derived 
from QCD in the weak coupling limit. However, the 
structure of the theory is completely determined by chiral 
symmetry, even if the coupling is not weak. In particular, 
there are no free parameters in the baryon coupling to the 
vector current. Mass terms are also strongly constrained 
by chiral symmetry. The effective chemical potentials 
(Xl^Xh) appear as left and right-handed gauge poten- 
tials in the covariant derivative of the nucleon field. We 
have 



DoN 
Fo 



doN + i[To,N], 



(106) 



where Xl = MMy{2pF) and Xr = M'fM/{2pF) as 
before. (Xl,Xii) covariant derivatives also appears in 
the axial vector current given in Eq. (|105p . 
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We can now study how the fcrmion spectrum depends 
on the quark mass. In the CFL state we have ^ = 1. 
For = the baryon octet has an energy gap A and 
the singlet has gap 2 A. The leading correction to this 
result comes from the commutator term in Eq. (|106p . We 
find that the gap of the proton and neutron is lowered, 
Ap^„ = A — fis, while the gap of the cascade particles 
E!~,E!'' is increased, Ag = A + /ig. As a consequence 
we find gapless (p, n) excitations at = A. This result 
agrees with the spectrum discussed in Sec. IIII.BI if the 
identification p = (bu) and n = (bd) is made. 

The situation is more complicated when kaon con- 
densation is taken into account. In the kaon con- 
densed phase there is mixing in the (p, E"*", E~,5~) and 
(n, E°, S°, A®, A°) sector. For rriKO < < A the spec- 
trum is given by 



A± i^i, 

A, 
2A. 



(107) 

Numerical results for the eigenvalues are shown in Fig.[8l 
We observe that mixing within the charged and neutral 
baryon sectors leads to level repulsion. There are two 
modes that become light in the CFL window < 2A. 
One mode is a charged mode which is a linear com- 
bination of the proton and the E+, while the other 
mode is a linear combination of the neutral baryons 
(n, S'', A^, A°). The charged mode becomes gapless 
first, at fis = 4A/3. Corrections to t his resu l t wer e 
studied in the NJL model calculation of (lForbesl . [2005l) . 
which includes various subleading condensates and ob- 
tains Us = 1.22A at /i = 500 MeV. The neutral mode be- 
comes gapless only at fi^ = 2 A. The most important dif- 
ference as compared to the spectrum in the gapless CFL 
phase without kaon condensation is that for fig < 2A 
only the charged mode is gapless. 



of the vacuum energy on the kaon current jp^ = V(/)a'- 
The gradient term in the meson part of the effective La- 
grangian gives a positive contribution 

= Ivlfbl- (109) 

A negative contribution can arise from gapless fermions. 
In order to determine this contribution we have to cal- 
culate the fermion spectrum in the presence of a nonzero 
current. The relevant couplings are obtained from the 
covariant derivative of the fermion field in Eq. (jl06p and 
the D and F-terms in Eq. (|104p . The fcrmion spectrum is 
quite complicated. The dispersion relation of the lowest 
mode is approximately given by 



= A - 



(^-^o)^ 
2A 



4V -Ja, 



(110) 



where I = -v ■ p — pp and we have expanded uJi near its 
minimum = (Ms + v-Jif )/4. Eq. (|110p shows that there 
is a gapless mode if fig > 4A/3 — The contribution 
of the gapless mode to the vacuum energy is 



Oq — ^ 



(Ill) 



where c?f2 is an integral over the Fermi surface. The en- 
ergy functiona l Sm + £q was analyzed in ( KrvievskH . l2005t 
ISchafeil . I2006D . There is an instability near the point 
/ig = 4 A/3. The instability is resolved by the forma- 
tion of a Goldstone boson current. If electric charge neu- 
trality is enforced the magnitude of the current is very 
small, and there is no tendency towards the generation 
of multiple currents. It was also shown that a ll gluonic 
screening masses are real ( Gerhold et aLl . [20071 ). The sit- 
uation is more complicated if the neutral fermion mode 
becomes gapless, too. In this case the magnitude of the 
current is not small, and multiple currents may appear. 
This regime corresponds to the portion of the cmGFlj-K^ 
curve in Fig. [3] that is only slightly (invisibly) below the 
gCFL curve. 



F. Goldstone boson currents 



In Sec. IIV.GI we showed that gapless fermion modes 
lead to instabilities of the superfiuid phase. Here we will 
discuss how these instabilities arise, and how they can be 
resolved, in the context of low energy theories of the CFL 
state, by the formation of the meson supercurrent state 
introduced in Sec. IIII.DI The chromomagnetic instabil- 
ity is an instability towards the spontaneous generation 
of currents, that is to say the spontaneous generation of 
spatial variation in the phase of the diquark condensate. 
Consider a spatially varying U{1)y rotation of the neu- 
tral kaon condensate 



^(x) = [/(x)eA-C/nx), 



(108) 



where = exp(i7rA4) and [/(x) = exp(i0A'(x)A8). This 
state is characterized by nonzero vector and axial-vector 
currents, see Eq. (|105p . We shall study the dependence 



G. Other effective theories 

Effective Lagrangians have also been been constructed 
for color superconducting phases other than the CFL 
phase. The effective theory for the light singlet axia l 
mode in the 2SC phase can be found in ( Beane et al\ . 
I2OOOI ). The phonon effective theory in the crystalline 
color superconducting phase is discussed in Sec. IVI.CI 

It is also interesting to study effective theories in 
QCD-like theories at large density. Some of these the- 
ories do not have a sign problem and can be stud- 
ied on the lattice with algorithms that are available 
today. Of particular inte r est are QCD with Nr. 



2 colors (lAlles et ali 120061; iFukushima and lidal. 



Hands et al. , 2006 , 1999 , Kogut et all 1200 ll . Il999l I200C 



2OO2I: iNishida et all l2004l) and 



2007 



densitv (iKoeut and Sinclair! . 1200^ 
I2OOII: ISplittorff et al\ . I2OOIL 



) CD at finite isospin 
[ ISon and Stephanovl . 
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VI. NJL MODEL COMPARISONS AMONG CANDIDATE 
PHASES BELOW CFL DENSITIES 

As we have explained in Sec. [ill ^it sufficiently high 
densities, where the up, down and strange quarks can be 
treated on an equal footing and the disruptive effects of 
the strange quark; mass can be neglected, quark matter is 
in the CFL phase. At asymptotic densities, the CFL gap 
parameter Acfl and indeed any property of CFL quark 
matter can be calculated in full QCD, as described in 
Sec. IIVI At any density at which the CFL phase arises, 
its low energy excitations, and hence its properties and 
phenomenology, can be described by the effective field 
theory of Sec. |Vl whose form is known and whose pa- 
rameters can be systematically related to the CFL gap 
^CFL- If we knew that the only form of color super- 
conducting quark matter that arises in the QCD phase 
diagram were CFL, there would therefore be no need to 
resort to model analyses. However, as we have discussed 
in Sec. IIIIl / (fiAcFh) may not be small enough (at 
fJ- = Mnuc where the nuclear — > quark matter transition 
occurs) for the QCD phase diagram to be this simple. 

Even at the very center of a neutron star, /x cannot be 
larger than about 500 MeV, meaning that the (density 
dependent) strange quark mass Mg cannot be neglected. 
In concert with the requirement that bulk matter must 
be neutral and must be in weak equilibrium, a nonzero 
Ms favors separation of the Fermi momenta of the three 
different flavors of quarks, and thus disfavors the cross- 
species BCS pairing that characterizes the CFL phase. If 
CFL pairing is disrupted by the heaviness of the strange 
quark at a higher fi than that at which color supercon- 
ducting quark matter is superseded by baryonic matter, 
the CFL phase must be replaced by some phase of quark 
matter in which there is less, and less symmetric, pairing. 

Within a spatially homogeneous ansatz, the next phase 
down in density is the gapless CFL (gCFL) phase de- 
scribed in Sec. IIII.BI However, as we have described 
in Sec. lIV.Gi such gapless paired states suffer from a 
chromomagnetic instability: they can lower their energy 
by the formation of counter-propagating currents. It 
seems likely, therefore, that a ground state with counter- 
propagating currents is required. This could take the 
form of a crystalline color superconductor, that we have 
introduced in Sec. IIII.CI Or, given that the CFL phase 
itself is likely augmented by kaon condensation as de- 
scribed in Sees. III. CI and IV. Dl it could take the form of 
the phase we have described in Sec. IV. Fl in which a CFL 
kaon condensate carries a current in one direction bal- 
anced by a counter-propagating current in the opposite 
direction carried by gapless quark quasiparticles. 

Determining which phase or phases of quark matter oc- 
cupy the regime of density between hadronic matter and 
CFL quark matter in the QCD phase diagram, if there is 
such a regime, remains an outstanding challenge. Barring 
a major breakthrough that would allow lattice QCD cal- 
culations to be brought to bear despite the fermion sign 
problem, a from-first-principlcs determination seems out 



of reach. This leaves two possible paths forward. First, 
as we describe in this section, we can analyze and com- 
pare many of the possible phases within a simplified few 
parameter model, in so doing seeking qualitative insight 
into what phase(s) are favorable. Second, as we shall 
describe in Sec. IVIIII we can determine the observable 
consequences of the presence of various possible color su- 
perconducting phases in neutron stars, and then seek to 
use observational data to rule possibilities out or in. 

A. Model, pairing ansatz, and homogeneous phases 

We shall employ a Nambu-Jona-Lasinio (NJL) model 
in which the QCD interaction between quarks is replaced 
by a point-like four-quark interaction, with the quantum 
numbers of single-giuon exchange, analyzed in mean field 
theory. This is not a controlled approximation. However, 
it suffices for our purposes: because this model has at- 
traction in the same channels as in QCD, its high density 
phase is the CFL phase; and, the Fermi surface splitting 
effects whose qualitative consequences we wish to study 
can be built into the model. Note that we shall assume 
throughout that Acfl <C /i. This weak coupling assump- 
tion means that the pairing is dominated by modes near 
the Fermi surfaces. Quantitatively, this means that re- 
sults for the gaps and condensation energies of candidate 
phases are independent of the cutoff in the NJL model 
when expressed in terms of the CFL gap Acfl: if the cut- 
off is changed with the NJL coupling constant adjusted 
so that Acfl stays fixed, the gaps and condensation en- 
ergies for the candidate crystalline phases also stay fixed. 
This makes the NJL model valuable for making the com- 
parisons that are our goal. The NJL model has two pa- 
rameters: the CFL gap Acfl which parametrizes the 
strength of the interaction and M^/(4/x), the splitting 
between Fermi surfaces in neutral quark matter in the 
absence of pairing. The free energy of candidate patterns 
of pairing can be evaluated and compared as a function 
of these two parameters. 

As a rather general pairing ansatz, we shall consider 

(ud) - A3 ^ exp (2iq^ r) 

a 

(us) A2 ^ exp (2iq2 • r) 

a 

(ds) ~ Ai ^ exp (2iq^ r) . (112) 

a 

If we set all the wave vectors to zero, wc can use 
this ansatz to compare spatially homogeneous phases in- 
cluding the CFL phase (Ai = A2 = A3 = Acfl), the 
gCFL phase (A3 > A2 > Ai > 0) and the 2SC phase 
(A3 = A2SC; Ai = A2 = 0). Choosing different sets of 
wave vectors will allow us to analyze and compare dif- 
ferent crystalline color superconducting phases of quark 
matter. 

NJL models of varying degrees of complexity have 
been used for a variety of purposes beyond the scope 
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of this review. For example, whereas we treat Aqfl 
and quark masses as parameters and use the NJL model 
to compare different patterns of pairing at fixed values 
of these parameters and /i, it is possible instead to fix 
the NJL coupling or couplings and then self-consistently 



functions of u (lAbuki et al. 


. 12005: Abuki and Kimihiro 


20061: Blaschke et all 2005 


; iBuballa and Oerte' 




2002 


loDolito et al. 


. 20071 Mishra and Mishra, 2004 


2006 


Riister et al\. 


2005: Steiner et al. 20021 Warringa 


,|2006|) 



Doing so reintroduces sensitivity to the cutoff in the 
NJL model and so does not actually reduce the 
number of parameters. Also, these models tend 
to find rather larger values of Ms than in analyses 
that go beyond NJL models, for exampl e the analy- 



sis us ing Dyson-Schwinger equations in (jNickel et al\ . 
l2006aD . There have also been many investigations 
of the phase diagram in the ^-T plane in NJL 



or solved for self-consistentlvl (1 Abuki and Kunihiro. 


20061: Barducci et all 2004 Berges and Rai 


agODal, 199£ 


Fukushima et all 2005: He et al. 2007 


; lida et al. 


2005.: .Kashiwa et al. 


. 20071: Mishra and Mishral. 


2005 


Riister et al. |2004 


20051. l2006l: ISchwarz et al\. 


1998 


Warringa et al. 2005 


). Although many of their features 



are sensitive to the cutoff as well as the chosen couplings, 
these NJL phase diagrams indicate how rich the QCD 
phase diagram may turn out to be, as different conden- 
sates vanish at different temperatures. One result that 
has been obtained using the Ginzburg-Landau approxi- 
mation as well as in NJL models and so is of more general 
validity is that upon heating the CFL phase at nonzero 
but small M^/^, as T increases A2 vanishe s first, then 
Al , and then A3 (IFukushima et all l2005t llida et all 



^1 , anu men ^3 iiruKuamma ei ut.i . izuum iiiua ei 
l2004l : lRiister et aal2004l ). However, it remains to be seen 
how this conclusion is modified by the effects of gauge- 
field fluctuations, which for Ms = turn the mean-field 
Ginzburg-Landau second order transition into a strong 
first order phase transition at a significantly elevated 
temperature, see Sec. IV. Bl and Eq. (|74p . 

We shall analyze quark matter containing massless u 
and d quarks and s quarks with an effective mass Ms. 
The Lagrangian density describing this system in the ab- 
sence of interactions is given by 



5"^**^- - M^f 



Q/3 



(113) 



where i, j = 1,2,3 are flavor indices and a, (3 = 1,2,3 are 
color indices and we have suppressed the Dirac indices, 
where M"-^ = (5"'^ diag(0, 0, Ms)y is the mass matrix, 
and where the quark chemical potential matrix is given 
by 



Q/3 



rriCtQ . rTiQ/3 



(114) 

with Qij = diag(2/3, — 1/3, — l/3)y the quark electric- 
charge matrix and T3 and Tg the diagonal color gen- 



erators. In QCD, /ie, /i3 and /ig are the zeroth com- 
ponents of electromagnetic and color gauge fields, and 
the gauge field dynamics ensure t hat they take on values 
such that the matter is neutral (lAlford and Raiagopal , 



2002 



2003 



Dietrich and Rischkd . I200I iGerhold and Rebhanl. 



Krvievskil . I2OO3I ). satisfying the neutrality condi 



tions ([3|). In the NJL model, quarks interact via four- 
fermion interactions and there are no gauge fields, so we 
introduce /ig, /is and by hand, and choose them to 
satisfy the neutrality constraints The assumption of 
weak equilibrium is built into the calculation via the fact 
that the only flavor-dependent chemical potential is /ie, 
ensuring for example that the chemical potentials of d 
and s quarks with the same color must be equal. Be- 
cause the strange quarks have greater mass, the equality 
of their chemical potentials implies that the s quarks have 
smaller Fermi momenta than the d quarks in the absence 
of BCS pairing. In the absence of pairing, then, because 
weak equilibrium drives the massive strange quarks to be 
less numerous than the down quarks, electrical neutrality 
requires a /ie > 0, which makes the up quarks less numer- 
ous than the down quarks and introduces some electrons 
into the system. In the absence of pairing, color neutral- 
ity is obtained with fJ-3 ~ fJ-s = 0. 

The Fermi momenta of the quarks and electrons in 
quark matter that is electrically and color neutral and in 
weak equilibrium are given in the absence of pairing by 





= /i ^ 


fJ-e 

3 


Pf 


= p- 


2/ie 

3 


Pf 




^ , Pe 


Pf 


= Me 


1 



(115) 



where we have simplified pp by working to linear order 
in /te and M^. To this order, electric neutrality requires 
/ie = M2/(4/i), yielding 



Pf 



Pf 
Pf 
Pf 



= A* 



Ml_ 

12/f 

Ml 

6/i 

5A/? 



12/j 



Pf 



Pf 



Ml 
4/i 

Ml 
"4^7 



Ml 
4/i ' 



(116) 



as illustrated in Fig. |2| We see from (|115|) that to leading 
order in and /ie , the effect of the strange quark mass 
on unpaired quark matter is as if instead one reduced 
the strange quark chemical potential by Mg/(2/i). We 
shall make this approximation throughout. Upon making 
this assumption, we need no longer be careful about the 
distinction between pF^s and /t's, as we can simply think 
of the three flavors of quarks as if they have chemical 
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potentials 



= Mm 

Mm = Pf 



with 



Ml 
8fi 



5ii , 



(117) 



(118) 



where the choice of subscripts indicates that 25^2 is the 
sphtting between the Fermi surfaces for quarks 1 and 3 
and 2(5/i3 is that between the Fermi surfaces for quarks 1 
and 2, identifying u, d, s with 1, 2, 3. 



As 



described 



2001 



Alford and Ra i agopaj 20021 : iRajagopal and Wilczek 



Alford et all 



20041 : 



Steiner g/J . 12002^ ). BCS pairing introduces 
quahtative changes into the analysis of neutrality. For 
example, in the CFL phase /ig = and fi^ is nonzero 
and of order M^/jj,. This arises because wherever BCS 
pairing occurs between fermions whose Fermi surface 
would be split in the absence of pairing, the Fermi 
momenta of these fermions are locked together. This 
maximizes the pairing energy gain while at the same 
time exacting a kinetic energy price and changing the 
relation between the chemical potentials and the particle 
numbers. This means that the /i's required for neutrality 
can change qualitatively as happens in the CFL example. 

The NJL interaction term with the quantum numbers 
of singie-gluon exchange that we add to the Lagrangian 
(fna is 



-'interaction 



(119) 



where we have suppressed the color and flavor indices 
that we showed explicitly in ()113p . and have continued 
to suppress the Dirac indices. The full expression for 
pA. jg (r^'^)„,,^, = 7"(7'^)a/3%. The NJL coupling 
constant A has dimension -2, meaning that an ultraviolet 
cutoff A must be introduced as a second parameter in 
order to fully specify the interaction. We shall define A 
as restricting the momentum integrals to a shell around 
the Fermi surface, fi — A < \p\ < fi + A. 

In the mean-field approximation, the interaction La- 
grangian (|119p takes on the form 

interaction = ^i^A{x)i^'^ + i^'^A(a;)V', (120) 

where A(a::) is related to the diquark condensate by the 
relations 



A{x) ^ ^Ar^^(vV'^)(r^,)^ 
A{x) = ^A(r^^)^(^^^)r^, 

= 7°A^(x)7° . 



(121) 



The ansatz (|112p can now be made precise: we take 

A{x) = Acf{x) (E) Cj^ , (122) 
with the color-flavor part 

3 

AcF{x)c.^M = ^ ^ A(q?)e2«i?-e,„^e/,j . (123) 

1=1 q? 

We have introduced notation that allows for the possibil- 
ity of gap parameters A(qj) with different magnitudes 
for different / and for different a. In fact, we shall only 
consider circumstances in which A(qJ) = A/, as in (|112p . 

The full Lagrangian, given by the sum of (|113|1 and 
p20[) . is then quadratic and can be written very sim- 
ply upon introducing the two component Nambu-Gorkov 
spinors (|12|) in terms of which 



c=hr^+^ A(.) 

2 i A{x) ■ 



(124) 



Here, /i = /i7o and /i is the matrix (jll4p . 

The propagator corresponding to the Lagrangian (jl24|) 
is given by 



{^{x)^{x')) 



i^PixMx')) {^{x)r {x')) 
{ij^{x)^{x')) (^^(x)V'^(x')) 

iG{x,x') iF(x,x') \ 
iF{x,x') iG(x,x') ) ' 



(125) 



where G and G are the "normal" components of the prop- 
agator and F and F are the "anomalous" components. 
They satisfy the coupled differential equations 



i^ + l^ A{x) \/G(x,x') F{x,x') 
A{x) {i^- iJiY )\ F{x,x') G{x,x') 



1 
1 



(126) 



(5W(.T-X') 



We can now rewrite (|12ip as 



A{x)^jXT^-'F{x,x){Ta.)' 



3i 



(127) 



A{x)^j\{T^yFix,x)TA. , 



either one of which is the self-consistency equation, or 
gap equation, that we must solve. 

Without further approximation, (|127p is not tractable. 
It yields an infinite set of coupled gap equations, one 
for each A(qj), because without further approximation 
it is not consistent to choose finite sets {q/}. When 
several plane waves are present in the condensate, they 
induc e an infinite tower of higher momentum conden- 
sates ( Bowers and Raiagopall [20021 ). In the next subsec- 
tion, we shall make a Ginzburg-Landau (i.e. small-A) 
approximation which eliminates these higher harmonics. 
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Of course, an even more dramatic simplification is ob- 
tained if we set all the wave vectors qf to zero. Still, even 
in this case obtaining the general solution with M,, ^ 
and A i ^ A2 ^ A3 is somewhat involv ed ( Alford et al\ . 
l2004bl . l2005d; iFukushima et all |2005| ). We shall not 
present the resulting analysis of the CFL^gCFL transi- 
tion and the gCFL phase here. The free energies of these 
phases are depicted in Fig. [31 and their gap parameters 
are depicted below in Fig. [101 

If we simplify even further, by setting Mg = and 
Ai = A2 = A3 = AcFL, the gap equation determin- 
ing the CFL gap para meter Aqfl can then be evalu ated 
analytically, yielding ( Bowers and Raiagopali . |2003 ) 



AcFL = 23 Aexp 



2/i2A 



(128) 



We shall see below that in the limit in which A <C 
AcFLji^M ^ results for the myriad possible crys- 

talline phases can be expressed in terms of Aqfl; nei- 
ther A nor A shall appear. This reflects the fact that 
in this limit the physics of interest is dominated by 
quarks near the Fermi surfaces, not near A, and so once 
AcFL is used as the parameter describing the strength 
of the attraction between quarks, A is no longer visi- 
ble; the cutoff A only appears in the relation between 
AcFL and A, not in any comparison among different pos- 
sible paired phases. We are using the NJL model in 
a specific, limited, fashion in which it serves as a two 
parameter model allowing the comparison among differ- 
ent possible paired phases at a given Acfl and Mg. 
NJL models have also been employed to estirnate the 
value of Aqft, at a given /x dXlford et al\. Il998l. IT999b : 



Berges and Raiagopal . Il999l: Carter and Diakonovl 1998 : 



Raiagopal and Wilczekl . 12005 " Rapp tt aLn'l998[ ): doing 
so requires normalizing the four-fermion interaction by 
calculating some zero density quantity like the vacuum 
chiral condensate, and in so doing introduces a depen- 
dence on the cutoff A. Such mean-field NJL analyses 
are important complements to extrapolation down from 
an analysis that is rigorous at high density and hence 
weak coupling, described in Sec. IIVI and give us con- 
fidence that we understand the magnitude of Acfl ~ 
10 — 100 MeV. This estimate rece ives further support 
from the lattice-NJL calculation of ( Hands and Waltersl . 
|2004[ ) which finds diquark condensation and a ~ 60 MeV 
gap in an NJL model whose parameters are normalized 
via calculation of and a constituent quark mass 

in vacuum. With these as inputs, A is then calculated on 
the lattice, i.e. without making a mean-field approxima- 
tion. With an understanding of its magnitude in hand, 
we shall treat Acfl as a parameter, thus making our 
results insensitive to A. 

We shall focus below on the use of the NJL model that 
we have introduced to analyze and compare different pos- 
sible crystalline phases, comparing their free energies to 
that of the CFL phase as a benchmark. The free energy 
of the 2SC phase is easily calculable in the same model, 
and the free energies of the unstable gapless CFL and 



gaples s 2SC phases can also be obtained ( Alford et al\ . 
l2005cl) . These free energies are all shown in Fig. [3l The 
free energies of phases with various patterns of single- 
flavor pairing have also been calculated in the same 
model I Alford et all . l2003l) . The NJL model is not a nat- 
ural starting point for an analysis of the kaon condensate 
in the CFL-A'° phase, but with considerable effort this 



has been accomplished in (iBuballal. [2005bl : iForbesl . [20051: 
iKleinhaus et al. L l2007l : IWarringal . l2006D . The curCFL-AfO 



phase of Secs. lIII.DI and lVJl in which the if ^-condensate 
carries a current, has not been analyzed in an NJL model. 
But, because both the CFL-A'° and curCFL-A'*^ phases 
are continuously connected to the CFL phase, they can 
both be analyzed in a model-independent fashion using 
the effective field theory described in Sec. [Vl The CFL- 
A'° and curCFL-AT" curves in Fig. [31 were obtained as 
described in Sec. [Vl It remains a challenge for future 
work to do a calculation in which both curCFL-AT'^ and 
crystalline phases are possible, allowing a direct compari- 
son of their free energies within a single calculation and a 
study of whether they are distinct as current results seem 
to suggest or are instead different limits of some more 
general inhomogeneous color superconducting phase. 



B. Crystalline phases 



Crystalline color superconductivitv (lAlford et al. 
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Leibovich et all 
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Mannarelfi et all l2006bl) 


naturally permits pairing be- 



tween quarks living at split Fermi surfaces by allowing 
Cooper pairs with nonzero net momentum. In three- 
flavor quark matter, this allows pairing to occur even 
with the Fermi surfaces split in the free-energetically op- 
timal way as in the absence of pairing, meaning that neu- 
tral crystalline phases are obtained in three-fiavor quark 
matter with the chemical potential matrix (|114p simpli- 
fied to fi = 3°"^ (g) diag {nu, iidi l^s) wit h the fiavor chem- 
ical potentials given sirnply by ( 1171) llCasalbuoni et al. , 
2005a : [Mannarelfi et all l2006bl : [Raiagopal and Sharmal . 



2006b[ ). up to hig her order corrections that have been 
investigated in (|Casalbuoni et all . [20061 ) . This is the ori- 
gin of the advantage that crystalline color superconduct- 
ing phases have over the CFL and gCFL phases at large 
values of the splitting 6fi. For example, by allowing u 
quarks with momentum p + q3 to pair with d quarks 
with momentum — p -|- qg, for any p, we can pair u and 
d quarks along rings on their respective Fermi surfaces. 
In coordinate space, this corresponds to a condensate of 
the form {ud) ~ A3exp(2zq3 • r). The net free energy 
gained due to pairing is then a balance between increas- 
ing I q3 1 yielding pairing on larger rings while exacting a 
greater kinetic energy cost. The optimum choice turns 
out to be Iqsl = ?/<5m3 with = 1.1997, corresponding 
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to pai ring rings on t he Fer mi surfaces with opening angle 
67.1° (lAlford a^.l ■ l2001a^ . Pairing with only a single qa 
is disadvantaged because the only quarks on each Fermi 
surface that can then pair are those lying on a single 
ring. This disadvantage can be overcome in two ways. 
First, increasing A widens the pairing rings on the Fermi 
surfaces into pairing ban ds which fill in, forming p airing 
caps, at large enough A (iMannarelh aLl . l2006bf ). Sec- 
ond, it is possible to cover larger areas of the Fermi sur- 
faces by allowing Cooper pairs with the same |q3 1 but var- 
ious qs, yielding (ud) ~ A3 exp(2 i q^-r) with the qg 

chosen from some specified set {q3,q3,q3, . . .} = {q3}. 
This is a condensate modulated in position space in some 
crystalline pattern, with the crystal structure defined by 
{qg}. In this two-flavor context, a Ginzburg-Landau 
analysis reveals that the best {qg} contains eight vectors 
pointing at the corners of a cube, say in the (±1, ±1, ±1) 
directions in momentum space, yielding a face-centered 
cubic structure in position space ( Bowers and Raiagopall 
I2OO2D . 

This subsection desc ribes the analysis of th r ee-flavo r 
crystalline phases in ( Raiagopal and Sharmal . l2006a[ ). 
We use the ansatz given by (|122p and (|123p for the three- 
flavor crystalline color superconducting condensate. This 
is antisymmetric in color {a, 13), spin, and flavor (i, j) in- 
dices and is a generalization of the CFL condensate to 
crystalline color superconductivity. We set Ai = 0, ne- 
glecting (ds) pairing because the d and s Fermi surfaces 
are twice as far apart from each other as each is from the 
intervening u Fermi surface. Hence, / can be taken to run 
over 2 and 3 only. {q2} and {q3} deflne the crystal struc- 
tures of the (us) and (ud) condensates respectively. We 
only consider crystal structures in which all the vectors 
in {q2} are equivalent to each other in the sense that any 
one can be transformed into any other by a symmetry op- 
eration of {qj} and same for {q3}. This justifies our sim- 
plifying assumption that the (us) and (ud) condensates 
are each specified by a single gap parameter (A2 and 
A3 respectively), avoiding having to introduce one gap 
parameter per q. We furthermore only consider crystal 
structures which are exchange symmetric, meaning that 
{qj} and {q3} can be exchanged by some combination 
of rigid rotations and reflections applied simultaneously 
to all the vectors in both sets. This simplification, to- 
gether with 6fi2 = SfJ-s (an approximation corrected only 
at order M^/^^), guarantees that we find solutions with 
A2 = A3. 

We analyze and compare candidate crystal structures 
by evaluating the free energy f2(A2,A3) for each crys- 
tal structure in a Ginzburg-Landau expansion in pow- 
ers of the A's. This approximation is controlled if 
A2, A3 ^ AcFL, <^Mj with AcFL the gap parameter in the 
CFL phase at M"^ / fi = 0. The terms in the Ginzburg- 
Landau expansion must respect the global U{\) symme- 
try for each fiavor, meaning that each A/ can only ap- 
pear in the combination |A/p. (The U{1) symmetries 
are spontaneously broken by the condensate, but not ex- 
plicitly broken.) Therefore, il(A2, A3) is given to sextic 



order by 

f^(A2,A3) 



/'2a2|A2p + P3a3|A3|2 



^(/32|A2|4 + /33|A3|4 + /332|A2nA3 

^(72|A2r + 73|A3r 



7322 I A3 n A2 I'' + 7233 1 A3 h A2 1 



(129) 



where we have chosen not ation consistent with that use d 
in the two flavor study of ( Bowers and Raiagopall . [2003 ). 
which arises as a special case of (|129p if we take A2 or 
A3 to be zero. Pi is the number of vectors in the set 
{q^^}. The form of the Ginzburg-Landau expansion (|129p 
is model- independent, whereas the expressions for the 
coefficients a/, /?/, /3/j, 7/, and 7/jj for a specific crystal 
structure are model-dependent. We calculate them in 
the NJL model described in Sec. IVI.Al For exchange 
symmetric crystal structures, a2 = a3 = a, f32 ~ P3 = f3, 
72 = 73 = 7 and 7233 = 7322 . 

Because setting one of the A/ to zero reduces the prob- 
lem to one with two-flavor pairing only, we can obtain a, 
P and 7 via applying the two-fl avor analysis described 
in ( Bowers and Raiagopall . [2002| ) to either {q2} or {q3} 
separately. Using a as an example, we learn that 



ai = a{qi,S^ii) 



2(7/ \qi-dfii 

2'°H4(g|-<^A*?) 



.(130) 



Here, qi = |q/| and A2SC is the gap parameter for 
the 2SC (2-fiavor, 2-color) BCS pairing obtained with 
6fij = and A/ nonzero with the other two gap param- 
eters set to zero. Assuming that Acfl ^ the 2SC 
gap p arameter is given by A2SC = 23Acfl ( Schafeil . 
l2000a[ ). see Sec. lIVI In the Ginzburg-Landau approxima- 
tion, in which the A/ are assumed small, we must first 
minimize the quadratic contribution to the free energy, 
and only then investigate the quartic and sextic contri- 
butions. Minimizing ai fixes the length of all the vectors 
in the set {q/}, and eliminates the possibility of waves 
at higher harmonics, yielding qj = rjSiJi wi th tj = 1.1997 



the so lution to ^ log [(r? + 1)/{t] - 1)] = 1 (|Alford et 
l2001al) . Upon setting qj = rjSfij, p30p becomes 



■log 



A^ 

^2SC 



4(5Aif (7]2 - 1) 



(131) 



Once the qi have been fixed, the only dimensionful 
quantities on which the quartic and sextic coefficients 
can depend are the Su r (iBowers and Raiagopall l2002t 
iRaiagopal and Sharmal . l2006bl ). meaning that for ex- 
change symmetric crystal structures and with Sfi2 — 
= we have /3 = jS/S^i'^, 13^2 = %2l^\i^, 7 = 
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^/S^JL^ and 7322 = 7322/1^^^ where the barred quanti- 
ties are dimensionless numbers which depend only on 
{€[2} and {€[3} that must be evaluated for each crystal 
structure. Doing so requires evaluating one-loop Feyn- 
man diagrams with 4 or 6 insertions of A/'s. Each 
insertion of A/ (AJ) adds (subtracts) momentum 2qj 
for some a. The vector sum of all these external mo- 
menta inserted into a given one-loop diagram must van- 
ish, meaning that the calculation consists of a book- 
keeping task (determining which combinations of 4 or 6 
q°'s selected from the sets {q^} satisfy this momentum- 
conservation constraint) that grows rapidly in complex- 
ity with the complexity of the crystal structure, and 
a loop integration that is nontrivial because the mo- 
me ntum in the propagator c hanges after each insertion. 
In ( Raiagopal and Sharmal . l2006bf ). this calculation is 
carried out explicitly for 11 crystal structures in the 
mean-field NJL model of Sec. IVI.AI upon making the 
weak coupling (AcpL and (5/i both much less than /i) 
approximation. Note that in this approximation neither 
the NJL cutoff nor the NJL coupling constant appear in 
any quartic or higher Ginzburg-Landau coefficient, and 
as we have seen above they appear in a only within Aqfl- 
Hence, the details of the model do not matter as long as 
one thinks of Aqfl as a parameter, kept <|; /i. 

It is easy to show that for exchange symmetric crys- 
tal structures any extrema of ri(A2,A3) in (A2,A3)- 
space must either have A2 = A3 = A, or have one o f 
A2 and A3 vanishing ( Raiagopal and Sharmal . l2006b[ ). 
It is also possible to show that the three-flavor crys- 
talline phases with A2 = A3 = A are electrically neutral 
whereas two-fiavor solutions in which only one of th e A's 
is nonzero are not ( Raiagopal and Sharmal . [2006tJ ) . We 
therefore analyze only solutions with A2 = A3 = A. 
We find that 0(A, A) is positive for large A for all 
the c rystal structures that have be en investigated to 
date ( Raiagopal and Sharmal . l2006b[ )."^ This allows us 
to minimize ri(A, A) with respect to A, thus evaluating 
A and fi. 

We begin with the simplest three-flavor "crystal" 
structure in which {q2} and {q3} each contain only 
a single vector, making t he {us) and {ud) conden sates 
each a single plane wave ([Casalbuoni et al liooli)- We 
call this the 2PW phase. Unlike in the more realis- 
tic crystalline phases we describe below, in this "crys- 
tal" the magnitude of the (ud) and (us) condensates 
are unmodulated. This simple condensate neverthe- 
less yields a qualitative lesson which proves helpful 
in winnowing the space of multiple plane w ave crys- 
tal structures ( Raiagopal and Sharmal . l2006bV For this 
simple "crystal" structure, all the coefficients in the 



2005aLlMannarelfi et aLl . l2006bt 
2006br ). The terms that occur 



Ginzb urg-Landau free energy can be evaluated analyti- 
callv (ICasalbuoni et al\ . 
iRaiagopal and Sharmal . 
in the three-fiavor case but not in the two-flavor case, 
namely /332 and 7322, describe the interaction between 
the two condensates, and depend on the angle </) be- 
tween q2 and qg. For any angle 0, both P32 and 7322 
are positive. And, both increase monotonically with 
4> and diverge as (j) tt This divergence tells us that 
choosing qj and qg precisely antiparallel exacts an infi- 
nite free energy price in the combined Ginzburg-Landau 
and weak-coupling limit in which A <^ 6fi, Acfl ^ 
meaning that in this limit if we chose = tt we find 
A = 0. Away from the Ginzburg-Landau limit, when 
the pairing rings on the Fermi surfaces widen into bands, 
choosing (f> = t: exacts a finite price meaning that A is 
nonzero but smaller than that for any other choice of 
(j). The high cost of choosing qj and qg precisely an- 
tiparallel can be understood qualitatively as arising from 
the fact that in this case the ring of states on the u- 
quark Fermi surface that "want to" pair with d-quarks 
coincides precis ely with the ring that " wants to" pair 
with s-quarks ( Mannarelli et all l2006bl ). This simple 
two plane wave ansatz has been analyzed upon making 
the weak-coupling approximation but without makin, 
the G inzburg-Landau approximation ( Mannarelli et al 
l2006bl ). All the qualitative lessons learned from the 
Ginzburg-Landau approximation remain valid and we 
learn further that the Gi nzburg-Landau approxirn ation 
always underestimates A ( Mannarelli et aUl2006tJ ). 



* This is in marked contrast with what happens with only two 
flavors (and upon ignoring the requirement of neutrality.) in 
that context, many crystal structures have negative 7 and 
henc e sextic order free energies t hat arc unbounded from be- 
low l lBowers and Raiagopall . l2002t) . 



The analysis of the simple two plane wave "crys- 
tal" structure, together with the observation that in 
more complicated crystal structures with more than 
one vector in {q2} and {q3} the Ginzburg-Landau co- 
efficient P32 (7322) is given in whole (in part) by a 
sum of many two plane wave contributions, yields one 
of two rules for constructing favorable crystal struc- 
tures fo^_thTee;^avp£_crYKte superconductiv- 
ity ( Raiagopal and Sharmal . l2006bt ): {q2} and {q3} 
should be rotated with respect to each other in a way that 
best keeps vectors in one set away from the antipodes of 
vectors in the other set. The second rule is that the sets 
{q2} and {q3} should each be chosen to yield crystal 
structures which, seen as separate two-flavor crystalline 
phases, are as favo rable as possible. The 11 crystal struc- 
tures analyzed in ( Raiagopal and Sharmal . l2006bl ) allow 
one to make several pairwise comparisons that test these 
two rules. There are instances of two structures which 
differ only in the relative orientation of {q2} and {q3} 
and in these cases the structure in which vectors in {q2} 
get closer to the antipodes of vectors in {q3} are disfa- 
vored. And, there are instances where the smallest angle 
between a vector in {q2} and the antipodes of a vector in 
{qg} are the same for two different crystal structures, and 
in these cases the one with the more favorable two-flavor 
structure is more favorable. These considerations, to- 
gether with explicit calculations, indicate that two struc- 
tures, which we denote "2Cube45z" and "CubeX" , are 



particularly favorable. 

Ill the 2Cube45z crystal, {q2} and {q^} each contain 
eight vectors pointing at the corners of a cube. If we ori- 
ent {q2} so that its vectors point in the (±1,±1,±1) 
directions in momentum space, then {q3} is rotated 
relative to {q2} by 45° about the z-axis. In this 
crystal structure, the (ud) and (us) condensates are 
each given by the most favored two -flavor crystal struc- 
ture ( Bowers and Raiagopall |2002[ ). The relative rota- 
tion maximizes the separation between any vector in {q2} 
and the nearest antipodes of a vector in {q^}. 

We arrive at the CubeX structure by reducing the 
number of vectors in {q2} and {q^^}- This worsens the 
two-flavor free energy of each condensate separately, but 
allows vectors in {q2} to be kept farther away from the 
antipodes of vectors in {q^}. We have not analyzed all 
structures obtainable in this way, but we have found one 
and only one which has a condensation energy compa- 
rable to that of the 2Cube45z structure. In the CubeX 
structure, {q2} and {qa} each contain four vectors form- 
ing a rectangle. The eight vectors together point toward 
the corners of a cube. The 2 rectangles intersect to look 
like an "X" if viewed end-on. The color, flavor and posi- 
tion space dependence of the CubeX condensate is given 
by 
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where a = VStt/q = 4. 536/(5 /i = 36.29/x/Afg is the lattice 
spacing. For example, with M^/fJ- = 100, 150, 200 McV 
the lattice spacing is a = 72,48,36 fm. We depict this 
condensate in Fig. [HI 

In Figs. [To] and [3l we plot A and versus / fi for 
the most favorable crystal structures that we have found, 
namely the CubeX and 2Cube45z structures described 
above. We have taken the CFL gap parameter Aqfl = 
25 MeV in these figures, but they can easily be re s caled t o 
any value of Aqfl M ( Raiagopal and Sharmal . [2"006bt ): 
if the A and / ji axes are rescaled by Acfl and the en- 
ergy axis is rescaled by A^pL. Fig. [TOl shows that the gap 
parameters are large enough that the Ginzburg-Landau 
approximation is at the edge of its domain of reliability. 
However, results obtained for the simpler 2PW crystal 
structures suggest that the Ginzburg-Landau calculation 
underestimates A and the condensation energy and that, 
even when it breaks down, it is a good qual i tative guide 
to the favorable structure ( Mannarelli et al ]. l2006bD . We 
therefore trust the result, evident in Fig. [31 that these 
crystalline phases are both impressively robust, with one 
or other of them favored over a wide swath of M'^/fi and 
hence density. We do not trust the Ginzburg-Landau 
calculation to discriminate between these two structures, 
particularly given that although we have a qualitative 
understanding of why these two are favorable we have 
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FIG. 9 The CubeX crystal structure of Eq. ([132)) . The fig- 
ure extends from to a/2 in the x, y and z directions. Both 
A2(r) and A3(r) vanish at the horizontal plane. A2(r) van- 
ishes on the darker vertical planes, and A3(r) vanishes on the 
lighter vertical planes. On the upper (lower) dark cylinders 
and the lower (upper) two small corners of dark cylinders, 
A2(r) = 4-3. 3A (A2(r) = -3.3A). On the upper (lower) 
lighter cylinders and the lower (upper) two small corners of 
lighter cylinders, A3(r) = -3.3A (A3(r) = -1-3. 3A). The 
largest value of |A/(r)| is 4A, occurring along lines at the 
centers of the cylinders. The lattice spacing is a when one 
takes into account the signs of the condensates; if one looks 
only at |A/(r)|, the lattice spacing is a/2. 



no qualitative argument for why one should be favored 
over the other. We are confident that 2Cube45z is the 
most favorable structure obtained by rotating one cube 
relative to another. We are not as confident that CubeX 
is the best possible structure with fewer than 8-1-8 vec- 
tors. Regardless, the 2Cube45z and CubeX crystalline 
phases together make the case that three-flavor crys- 
talline color superconducting phases are the ground state 
of cold quark matter over a wide range of densities. If 
even better crystal structures can be found, this will only 
further strengthen this case. 

Fig.[3|shows that over most of the range of //i where 
it was once considered a possibility, the gCFL phase can 
be replaced by a much more favorable three-flavor crys- 
talline color superconducting phase. We find that the 
two most favorable crystal structures have large conden- 
sation energies, easily 1/3 to 1/2 of that in the CFL phase 
with Ms = 0, which is 3A^pL/i^/7r^. This is at first sur- 
prising, given that the only quarks that pair are those 
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FIG. 10 Gap parameter A versus M'^/fi for: the CFL gap 
parameter (set to 25 MeV at / fJ. — 0), the three gap 
parameters Ai < A2 < A3 describing (ds), (us) and {ud) 
pairing in the gCFL phase, and the gap parameters in the 
crystaUine color superconducting phases with CubeX and 
2Gube45z crystal structures. Increasing M^/jj. corresponds 
to decreasing density. 



or other is favored over the CFL, gCFL and unpaired 
phases for 2.9Acfl < M^/f^ < 10.4Acfl, as shown 
in Fig. El For Aqfl = 25 MeV and Ms = 250 MeV, 
this translates to 240MeV < ^ < 847MeV. With these 
choices of parameters, the lower part of this range of fi 
(higher part of the range of / fi) is certainly super- 
seded by nuclear matter. And, the high end of this range 
extends beyond the ^ ~ 500 MeV characteristic of the 
quark matter at the densities expected at the center of 
neutron sta rs. This qualitative feature persists in the 
analysis of ( Ippolito et all l2007t ) in which Mg is solved 
for rather than taken as a parameter. If neutron stars do 
have quark matter cores, then, it is reasonable to include 
the possibility that the entire quark matter core could 
be in a crystalline color superconducting phase on the 
menu of options that must ultimately be winnowed by 
confrontation with astrophysical observations. (Recall, 
that if AcFL is larger, say ^ 100 MeV, the entire quark 
matter core could be in the CFL phase.) As we shall see 
in the next subsection, crystalline color superconducting 
quark matter is rigid, with a very large shear modulus, 
while at the same time being superfluid. This provides a 
possible origin for pulsar glitches, as we shall discuss in 
Sec. IVIILFI 



lying on rings on the Fermi surfaces, whereas in the CFL 
phase with Ms ~ pairing occurs over the entire u, d 
and s Fermi surfaces. It can to a degree be understood 
qualitatively once we recall that there are in fact many 
rings, and note that as A increases, the pairing rings 
spread into bands on the Fermi surfaces, and for A as 
large as that we find to be favored these bands have ex- 
panded and fille d in, becoming many "pol ar caps" on the 
Fermi surfaces ( Mannarelli et al I l2006bD . In addition 
to being free-energetically favorable, these crystalline 
phases are, as far as we know, stable: t hey do not suf- 
fer from the chromomagnetic instability fCi minale et ai . 
2006; Gatto and Ruggicrj, ,2007: Giannakis et all 12005 : 



Giannakis and Renl . l2005bl ) an d they are also s table with 



respect to kaon condensation ( Anglani et all [20071 ). In 
simplified analogue contexts, it has even been possible 
to trace the path in configuration space from the unsta- 
ble gapless phase (analogue of gCFL) downward i n free 
energy to the stable crysta lline phase (|Fukushimal . 120061 : 
IFukushima and lidal . [20071 ) . 

Fig. [3] also shows that it is hard to find a crystalline 
phase with lower free energy than the gCFL phase at 
the lower values of M^//i (highest densities) within the 
"gCFL window" . At these densities, however, the calcu- 
lations described in Sec. |V| demonstrate that the gCFL 
phase is superseded by the stable CFL-K^ and curCFL- 
phases, as shown in Fig. [31 

The three-flavor crystalline color superconducting 
phases with CubeX and 2Cube45z crystal structures are 
the lowest free energy phases that we know of, and hence 
candidates for the ground state of QCD, over a wide 
range of densities. Within the Ginzburg-Landau approx- 
imation to the NJL model that we have employed, one 



C. Rigidity of crystalline color superconducting quark 
matter 

The crystalline phases of color superconducting quark 
matter that we have described in the previous subsec- 
tion are unique among all forms of dense matter that 
may arise within neutron star cor es in one respect: they 
are rigid ( Mannarelli et all l2007l ) . They are not solids 
in the usual sense: the quarks are not fixed in place 
at the vertices of some crystal structure. Instead, in 
fact, these phases are superfluid since the condensates all 
spontaneously break the U{1)b symmetry corresponding 
to quark number. We shall always write the conden- 
sates as real. This choice of overall phase breaks U{1)b, 
and spatial gradients in this phase correspond to super- 
currents. And yet, we shall see that crystalline color 
superconductors are rigid solids with large shear mod- 
uli. The diquark condensate, altho ugh spatially inhomo- 



geneous, can carry sup ercurrents ( Alford et all l2001al : 
iMannarelli et a/.l . l2007t ). It is the spatial modulation of 
the gap parameter that breaks translation invariance, as 
depicted for the CubeX phase in Fig. [21 and it is this 
pattern of modulation that is rigid. ^ This novel form of 
rigidity may sound tenuous upon first hearing, but we 
shall present the effective Lagrangian that describes the 



' Supersolids jAndreey and Lifshitl. Il969l: IChesteJ. Il970l : 
iKim and Chaiil . l2004al lU iLeggettJ . Il97d Isonl . |2005| ) are another 
example of rigid superfluids, but they differ from crystalline color 
superconductors in that they are rigid by virtue of the presence 
of an underlying lattice of atoms. 
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phonons in the CubeX and 2Cube45z crystalline phases, 
whose lowest order coefficients have been calculated in 



the N JL model that we are employing (jMannarelli et ali . 

We shall then extract the shear moduli from the 
phonon effective action, quantifying the rigidity and in- 
dicating the presence of transverse phonons. The fact 
that the crystalline phases are simultaneously rigid and 
superfluid means that their presence within neutron star 
cores has potentially observable consequences, as we shall 
describe in Sec. IVlll.FI 

The shear moduli of a crystal may be extracted 
from the effective Lagrangian that describes phonons 
in the crystal, namely space- and time-varying dis- 
placements of the crystalline pattern. Phonons in 
two-flavor crystall i ne pha ses were first investigated in 
( Casalbuoni et al\ . I2002al l3). In the present context, we 
introduce displacement fields for the (ud), (us) and (ds) 
condensates by making the replacement 



AiY^e''^'- 



,2iqJ-(r-uj(r)) 



(133) 



in p23|l . One way to obtain the effective action describ- 
ing the dynamics of the displacement fields u/(r), in- 
cluding both its form and the values of its coefficients 
within the NJL model that we are employing, is to be- 
gin with the NJL model of Sec. IVI.AI but with and 
integ rate out the fermion fiel ds. After a lengthy calcula- 
tion ( Mannarelli et all l2007l ) . this yields 



5[u] = i J d^xY^Kj 



(134) 



where m, n, v and w are spatial indices running over x, 
y and z and where we have defined 
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Upon setting Ai = and A2 = A3 = A 
2Ai2|A|V 



(135) 
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7r2(r/2 - 1) 



0.664 /i^lA^I . (136) 



S[u] is the low energy effective action for phonons in any 
crystalline color superconducting phase, valid to second 
order in derivatives, to second order in the gap param- 
eters A J and to second order in the phonon fields u/. 
Because we are interested in long wavelength, small am- 
plitude, phonon excitations, expanding to second order in 
derivatives and in the phonon fields is satisfactory. More 
complicated terms will arise at higher order, for exam- 
ple terms that couple the different u/ 's, but it is legit- 
imatc to neglect these complications ( Mannarelli et all . 



I2OO7I ). Extending this calculation to higher order in the 
Ginzburg-Landau approximation would be worthwhile, 
however, since as we saw in Sec. IVI.BI this approxima- 
tion is at the edge of its domain of reliability. 

In order to extract the shear moduli, we need to com- 
pare the phonon effective action to the gen eral theory 
of elastic media ( Landau and Lifshitd . Il98lh . which re- 
quires introducing the strain tensor 



(137) 



2 V dx"" dx"" 
We then wish to compare the action (|134p to 

S[n\ = I fd^x(Y.Y.PTi9ouTKdouT) 

which is the general form of the action in the case in 
which the effective action is quadratic in displacements 
and which defines the elastic modulus tensor A™"^"™ for 
this case. In this case, the stress tensor (in general the 
derivative of the potential energy with respect to s™") is 
given by 



= A 



mvnw ^nw 



(139) 



The diagonal components of a are proportional to the 
compression exerted on the system and arc therefore re- 
lated to the bulk modulus of the crystalline color super- 
conducting quark matter. Since unpaired quark matter 
has a pressure ~ /i^, it gives a contribution to the bulk 
modulus that completely overwhelms the contribution 
from the condensation into a crystalline phase, which is 
of order fi^A^. We shall therefore not calculate the bulk 
modulus. On the other hand, the response to shear stress 
arises only because of the presence of the crystalline con- 
densate. The shear modulus is defined as follows. Imag- 
ine exerting a static external stress aj having only an 
off-diagonal component, meaning cr™^' 7^ for a pair of 
space directions m =/= v, and all the other components of 
a are zero. The system will respond with a strain s"™. 
The shear modulus in the mv plane is then 



mil 
mv / 



^ \ mvmv 

2^' 



(140) 



where the indices m and v are not summed. For a general 
quadratic potential with cr™" given by (|139p , i^}™ simpli- 
fies partially but the full simplification given by the last 
equality in (|140|1 only arises for special cases in which the 
only nonzero entries in A"*""™ with m ^ v are the A™^™" 
entries, as is the case for all the crystal structures that 
we consider. 

For a given crystal structure, upon evaluating the sums 
in (|134p and then using the definition p37p to compare 
P34p to (|138p , wc can extract expressions for the A tensor 
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and thence f or the shear moduh. Thi s analysis, described 
in detail in ( Mannarelli et al ],[2003), shows that in the 
CubeX phase 
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while in the 2Cube45z phase 
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We shall see in Sec. I VIII. F I that it is relevant to check 
that both these crystals have enough nonzero entries in 
their shear moduli that if there are rotational vortices 
are pinned within them, a force seeking to move such a 
vortex is opposed by the rigidity of the crystal structure 
described by one or more of the nonzero entries in the 
z^/. This is demonstrated in ([Man narel li et a/.„ ,20071 ). 

We see that all the nonzero shear moduli of both the 
CubeX and 2Cube45z crystalline color superconducting 
phases turn out to take on the same value. 



16 



(143) 



with K defined by (|136p . Evaluating k yields 
l^CQM = 1.18 /i^A^ 



2.47 



MeV 



fm^ V 10 MeV 



400 MeV 



(144) 



From (|144p we first of all see that the shear modu- 
lus is in no way suppressed relative to the scale ^^A^ 
that could have been guessed on dimensional grounds. 
And, second, we discover that a quark matter core in 
a crystalline color superconducting phase is 20 to 1000 
times more rigid than t h e crus t of a conventional n eutro n 
star (iMannarelh et all 1200 7*: 'Stro hmaver et all Il99l[ ). 
Finally, see ( Mannarelli et a l.. 200^ for the extraction 
of the phonon dispersion relations from the effective ac- 
tion (|134p . The transverse phonons, whose restoring 
force is provided by the shear modulus and which cor- 
respond to propagating ripples in a condensation pat- 
tern like that in Fig. [SI turn out to have direction- 
dependent velocities that are typically a substantial frac- 
tion of th e speed of light, in the specific instances eval- 
uated in (jMannarelli et all |2007[ ) being given by a/I/S 
and \/2/3. This is yet a third way of seeing that this 
superfluid phase of matter is rigid indeed. 



VII. TRANSPORT PROPERTIES AND NEUTRINO 
PROCESSES 



In Sec. IVIIII we shall discuss how the observation of 
neutron star properties constrains the phase structure of 



dense quark matter. A crucial ingredient in these analy- 
ses are the transport properties as well as neutrino emis- 
sivities and opacities of different phases of quark matter. 

Using the methods introduced in Sec. |V] it is possi- 
ble to perform rigorous calculations of transport proper- 
ties of the CFL phase. The results are parameter free 
predictions of QCD at asymptotically large density, and 
rigorous consequences of QCD expressed in terms of a 
few phcnomenological parameters (/^r, • ■ •) at lower 
density. 

In the case of other color superconducting phases we 
perform calculations using perturbative QCD or models 
of QCD. For many quantities the results depend mainly 
on the spectrum of quark modes, and not on details of 
the quark-quark interaction. 



A. Viscosity and thermal conductivity 

Viscosity and thermal conductivity determine the dis- 
sipated energy E 'm & fluid with nonzero gradients of the 
velocity v and the temperature T, 



E 



T 



d^x{d,T)\ (145) 



The transport coefficients 77, C, and ft are the shear and 
bulk viscosity and the thermal conductivity, respectively. 
Eq. p45p is strictly valid only for non-relativistic fluids. 
In the case of relativistic fluids there is an extra contri- 
bution to the dissipated energy which is p roportional to 
K and the gradient of fj. (|Landau and Lifs'h itz. 1987.) • In 
terms of its hydrodynamic properties a superfluid can be 
viewed as a mixture of a normal and a superfluid com- 
ponent characterized by separate flow velocities. The 
shear viscosity is entirely due to the normal component, 
but there are contributions to the bulk viscosity which 
are related to st resses in the superfluid flow r elative to 
the normal one ( Andersson and Comeij . l2006t iGusakovl 
I2OO7I: iKhalatnikovL Il989ir In the following we shall ne- 
glect these effects and interpret Vi in Eq. (|145[) as the 
normal fluid velocity. 

In neutron stars an important contribution to the 
bulk viscosity arises from electroweak effects. In a bulk 
compression mode the density changes periodically and 
electroweak interactions may not be sufficiently fast to 
reestablish weak equilibrium. Weak effects occur on the 
same time scale as the oscillation period of the neutron 
star and the frequency dependence of the bulk viscosity 
is important. We define 



C(^) = 2(i^ 



SVo 



(146) 



where u; is the oscillation frequency, (...) is a time aver- 
age, and SVq/Vo is the fractional change in the volume. 
The coefficient in Eq. p45p is the u ^ limit of ({uj). 



47 



10 
10 



30 



H = 400 MeV M = 120 MeV t = 1 ms 



28 



10 



26 



24 




22 



10 
10 



-is 20 

10 



E 
o 



-2 



10 
10 
10^ 
10 
10 



16 



12 



10 



- - - " / ICF kaonsi / - 



0.01 




T (MeV) 



FIG. 11 (Color online) Bulk viscosities as functions of tem- 
perature for an oscillation period t = l-Kjuj = \ ms. CFL 
phase: contribution from the process ^ ip + {f for different 
values of 5m = m^o — tJ'K'> contribution from </p <^ <^ + (^, 
see Eq. (|150|l . 2SC phase and unpaired quark matter: contri- 
bution from the process -I- d <-> u-\- s. 



If a single weak process is responsible for reestablishing 
chemical equilibrium, the frequency dependent bulk vis- 
cosity can be written in the form 



C(^) = G 



7 



7^ + CJ^ 



(147) 



The prefactor C accounts for the dependence of the equi- 
librium densities (e.g., the net difference between the den- 
sity of strange and non-strange quarks if the weak process 
changes strangeness) on the respective chemical poten- 
tials, and 7 is the characteristic inverse time scale of the 
flavor changing process. Eq. p47p shows that, for a given 
w, C has a maximum at 7 = a;. At this point the time 
scale of the microscopic process matches the one of the 
external oscillation. If more than one weak process con- 
tributes t o reequilibration, Eg. fl47l) becomes more com- 
plicated (lAlford and Schmittl . [2007 : iHaensel et al. I. I2000I: 
ISa'd et all l2007bD 



1. CFL phase 

The normal fluid is composed of quasi-particle exci- 
tations. In the CFL phase all quark modes are gapped 
and the relevant excitations are Goldstone bosons. At 
very low temperature, transport properties are domi- 
nated by the massless Goldstone boson (p associated with 
the breaking of the U{1)b symmetry. Using the results 
in Sec. IV.C.2[ we can compute the mean free path l^p 
of the (p due to <-> if + (p and f + if <-> <^ + V? scat- 
tering. Small angle scattering contributions give rise to 



cx fi*/T^ (|Manuel et all l2005f) and ~ 1 km at 
T = 0.1 MeV, while large angle sca ttering contributions 
yield an even longer cx ^ /T^ (IShovkovv and Ellis! . 
I2OO2). The thermal conductiv ity n due to the (p is given 
bv (IShovkovv and Ellij ■ l2002^ 



(148) 



where is the ip mean-free path between large an- 
gle scatterings and v is the ip velocity from Eqs. ((83|) 
and ([84]) . For temperatures below ^ 1 MeV the ther- 
mal conductivity is very large and macroscopic amounts 
of CFL matter are expected to be isothermal. The 
electric conductivity in CFL matter is dominated by 
thermal electrons and po sitrons and was estimated in 
(IShovkovv and EUisl . [200I . 

At low temperatures, the shear viscosity of the CFL 
phase is dom inated by the ip co ntribution, which was 
computed in ( Manuel et all l2005l ) and is given by 



77 = 1.3 X 10-4^ 



(149) 



The bulk viscosity C vanishes in an exactly scale in- 
variant system. For realistic quark masses the dom- 
inant source of scale breaking is the strange quark 
m a,ss. The contribution from the p rocess (p (p + (p 
is (|Manuel and Llanes-Estradal [20071 ) 



^^0.011^ 



(150) 



We show this contribution in Fig. [TTJ The other contri- 
bution to the CFL bulk viscosity presented in the figure 
comes from the process ip + ip and was studied for 

arbitrary lu in jAlford et all l2007al) . We observe that 
at T ~ (1 - 10) MeV the bulk viscosity of CFL mat- 
ter is comparable to that of unpaired quark matter. For 
T < 1 MeV, C is strongly suppressed. Depending on the 
poorly known value for Sm = 771^0 — fJ-KO (here assumed 
to be positive, a negative value corresponds to kaon con- 
densation), the pure ip contribution given in Eq. p50p 
may dominate over the K'^ ^ ip + ip reaction at low 
enough temperatures. However, for T < 0.1 MeV the ip 
mean free path is on the order of the size of the star, i.e., 
the system is in the coUisionless rather than in the hydro- 
dynamic regime, and the result ceases to be meaningful. 

Thermal conductivity and viscosities for the CFL-/v ° 
phase have not yet been computed. The existence of 
a gapless K'^ Goldstone mode in this phase will intro- 
duce new contributions. However, since the CFL results 
for K and 77 are already dominated by a gapless mode, 
namely the ip, the modifications to these quantities are 
not expected to be significant. The modification to C will 
be more significant, since the kaon contribution to this 
quantity is already important in the CFL phase. 
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FIG. 12 Contributions to the process u + d-^u + s in the 
2SC phase. A gapped fermion is marked with the gap A at 
the respective hne. (We have omitted (smaU) contributions 
from anomalous propagators.) 



2. Other phases 

For unpaired, ultrarelativistic three-flavor quark mat- 
ter, thermal and electric conductivity as well as shear vis- 
cosity have been computed in (Hcisel berg and Pethickl . 
Il993 ). In the low-temperature limit (in particular 



T ^ mo with the electric 
A^/gVV(27i"^)) they are 



screening mass mjj = 



0.5- 



M. 

v2 



Ccl 



0.01 



2 2 2/3 



a 



2215/3 



and 



77 



4.4 X 10 



4 2/3 
-3 "^g 
0,2715/3 



(151) 



(152) 



These quantities have not yet been computed for partially 
gapped color superconductors such as the 2SC phase. 
The presence of ungapped modes, however, suggests that 
the results only differ by a numerical factor from the un- 
paired phase results. 

The dominant flavor changing process that contributes 
to the bulk viscosity in unpaired quark rn atter is the re- 
action (jAnand et all |2000| - iMadsenl . Il992l) 



u + s . 



(153) 



Other relevant processes are the semi -leptonic processes 
u + e d + and u + e ^ s + v,. ( Dong etoll . l2007t 



ISa'd et aLl . r2007b[ ). 

In a partially gapped phase the bulk viscosity is also 
dominated by the process p53p . In the 2SC phase of 
three flavor quark matter, the number of d-quarks pro- 
duced per unit time and volume, F, due to (|153p can 
be computed from the diagrams shown in Fig. 1121 The 
combinatorical factors in front of the diagrams are ob- 
tained upon counting color degrees of freedom: one can 
attach one of three colors to each of the two weak ver- 
tices, giving rise to 9 possibilities. In the 2SC phase all 
blue quarks and all strange quarks are unpaired while 



all other modes are paired, see Table H] Consequently, 
4 of the 9 possibilities contain three gapped modes (red 
or green for both vertices), 2 contain two gapped modes 
(red or green for one, blue for the other vertex), 2 con- 
tain one gapped mode (blue for one, red or green for 
the other vertex), and one contains only unpaired modes 
(blue for both vertices). Therefore, at very low tempera- 
ture, r ^ A, where the contributions of gapped quarks 
arc exponentially suppressed, F i s to a very good approx- 
imation given by ( Madsenl . l2000t ) 



F 



2SC 



_J_-punp 



for T < A . 



(154) 



since only the one reaction containing only unpaired 
modes contrib utes. The rate F^"p was computed in 
(IMadsenl . [T993I) . 

For larger temperatures, the contribution from gapped 
modes cannot be neglected. Each diagram yields a contri- 
bution which (for one direction of the process) schemat- 
ically reads 



F cx 



/ J^(5(ei£i + e2£2 - 6363 - 6424 + (5/^) 



X /(eiei)/(e2e2)/(-e3e3)/(-e4e4) 



(155) 



Here, £i are the quasiparticlc energies, 6fj, = iJ,s — fJ-d, and 
/ is the Fermi distribution function. ^ is a function of 
the momenta Pi and the signs ~ ±1. The sum over the 
signs Ci is very important in a paired system: the process 
u + d u + s not only receives contributions from 2 — > 2 
processes, but also from 3^1 and 1^3 reactions 
involving pairs created or absorbed by the condensate. 

From the net production rate of d quarks F one obtains 
the characteristic inverse time scale 7 needed for the bulk 
viscosity in Eq. (|147p . For small external volume oscil- 
lations SVo/Vq, F is linear in the resulting oscillation in 
chemical potentials, F = XSfi. Then, 7 = BX, where 
B depends only on the equilibrium flavor densities. The 
resulting bulk viscosity as a function of temperature for 
a typical oscillation frequency uj/{27t) = lms~^ is shown 
in Fig. [TlJ A critical temperature of Tc = 30 MeV is 
assumed. For low temperatures, the time scale of the 
nonlcptonic process is much smaller than the oscillation 
frequency 7 ^ oj, implying ^ cx 7. Consequently, from 
Eq. (|154p we conclude C2SC — Cunp/9. For large tem- 
peratures, however, we have 7^0; and thus ^ oc I/7. 
Consequently, the superconducting phase, which has the 
slower rate, has the larger bulk viscosity. 

The bulk viscosity has also been comput ed for two- 
flavor quark matter with single- flavor pairing ( Sa'd et 'all . 
l2007al) . In this case there are also ungapped modes and 
thus the result is similar to the one of the 2SC phase. 
The main difference is the lower critical temperature for 
single-flavor pairing. As a consequence, these phases are 
unlikely to exist for temperatures larger than that at 
which the bulk viscosity of the unpaired phase is maxi- 
mal. Therefore, the bulk viscosity cannot be larger than 
that of the unpaired phase. 
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B. Neutrino emissivity and specific heat 

Neutrino emissivity determines the rate at which quark 
matter can loose heat via neutrino emission. For the 
purpose of studying how neutron stars with ages rang- 
ing from tens of seconds to miUions of years cool, as we 
shall discuss in Sec. IVIII.Cl it is appropriate to treat the 
matter as completely transparent to the neutrinos that 
it emits. 



1. CFL phase 

In CFL quark matter, all quasifermion modes are 
gapped and neutrino emissivity is dominated by reac- 
tions involving (pseudo)-Goldstone modes such as 



± 



(f + (f ip - 



Vp + Vf 



(156a) 
(156b) 
(156c) 



These processes wer e studied in ( Jaikumar et aZ.1 . 120021 : 
iReddv et al. 1, l2003bf ). The decay rates of the massive 
mesons tt^, K^, and tt'^ are proportional to their num- 
ber densities and are suppressed by Boltzmann factors 
exp(— i?/T), where E is the energy gap of the meson. 
Since the pseudo-Goldstone boson energy gaps are on 
the order of a few MeV, the emissivities are strongly 
suppressed as compared to unpaired quark matter for 
temperatures below this scale. Neutrino emission from 
processes involving the ip is not exponentially suppressed, 
but it involves a very large power of T, 



2, ,4 



(157) 



and is numerically very small. Reddy et al. also studied 
the neutrino mean free path l^. For T ^ 30 MeV the 
mean free p ath is on the order of 1 m, but for T < 1 MeV, 
> 10 km (|Reddv et aZ.l . l2003bl ). In the CFL-K° phase, 
is almost the same as in the CFL phase, wh ile the 
neutrino emissivity is larger ( Reddv et aZ.I . l2003al) . 

The specific heat of CFL matter is also dominated by 
the if, yielding 



Cv 



2n' 
15^ 



(158) 



This is much smaller than the specific heat of any phase 
containing unpaired quarks, as we shall see below. 



2. Other phases 

The density of thermally excited ungapped fermions 
is proportional to /x^T while that of ungapped bosons is 

. This means that in any degenerate system (T <C /i) 
ungapped fermion modes, if they exist, will dominate 



the neutrino rates. In unpaired quark matter neutrino 
emissivity is dominated by the direct Urea processes 



- e 
d 



(electron capture) , (159a) 
(/3-decay) . (159b) 



The radiated e nergy per unit of time and volume is 
(llwamotol [T980l ) 



457 
630 



(160) 



Note that this result is proportional to the strong cou- 
pling constant a^. The tree- level processes for mass- 
less quarks are approximately coUinear and the weak 
matrix element vanishes in the forward direction. A 
nonzero emissivity arises from strong interaction cor- 
rections, which depress quark Fermi momenta relative 
to their chemical potentials. Because they do not at 
the same time depress the electron Fermi momentum, 
this opens up phase space for the reactions p59p . A 
nonzero emissivity can also arise from quark mass ef- 
fects, or higher order corrections in T/j^i. Since strange 
quark decays are Cabbibo suppressed and T/ii is small 
the dominant contribution is likely to be that propor- 
tional to as, namely (|160p . Note that we have not 
i ncluded non-Fermi liquid cor rections of 0{aslog{T)) 
(jSchafer and Schweiize^ . l2004al) . 

In order to determine the rate at which neutron stars 
cool we also need to know the specific heat. In unpaired 
quark matter 



NrN 



Cv 



(161) 



where we have again neglected terms of 0(aslog(r)) 
( Idp et gTI . |2004| ) and assumed the flavor chemical po- 
tentials to be equal. We see that the specific heat (|158p 
in the CFL phase, whose excitations are bosonic, is much 
smaller than that in unpaired quark matter. 

In the case of 2SC matter, the neutrino emissivity at 
low temperature is 1/3 of that of unpaired quark matter. 
The 2SC emissivity for a r bitrary temperatures can be 
found in ( Jaikumar et all l2006b[ ). In addition to the 
direct Urea process neutrino pair production 



q + q^q + q + iye + i'e 



(162) 



(q is any quark flavor and £ denotes neutrino flavor ) 
has also been studied ( Jaikumar and Prakashl l200l[ ). 
The rate of this process is parametrically smaller than 
the direct Urea process for very small temperatures 
(exp(— 2A/r) vs. exp(— A/T)), but it may play a signif- 
icant role for temperatures close to the superconducting 
phase transition temperature Tc- 

For the LOFF phase similar arguments apply. The 
presence of ungapped m odes renders its specific heat 
( Casalbuoni et all . |2003( ) and its neutrino emissivity due 



to direct Urea process es virtually ind i stingu ishable from 
the unpaired phase ( Anglani et all . l2006f) . However, 
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interesting effects of crystalline structures may be ex- 
pected for other cooling mechanisms. This is not un- 
like effects in the crust of a conventional neutron star, 
where for instance electron-phonon scattering as well as 
Bragg diffraction of electrons lea d to neutrino emission 
via bremsstrahlung processes, see ( Yakovlev et aZj . boOlf ) 
and references therein. 

The direct Urea processes have also been considered for 
the gapless CFL phase. A distinctive feature of this phase 
is the fact that the energy of one of the quark modes 
is approximately quadratic in momentum. This implies 
a strong enhancement in the specific heat, which leads 
to very slow cooling at very small temperatures when 
phot on emission f r om the surface dominates the energy 
loss ( Alford et al\ . l2005lJ ). However, the instability of 
this phase at small temperatures suggests that this result 
is most probably of no relevance for astrophysics. 

Finally, it is interesting to consider neutrino emission 
from single flavor paired matter. Single flavor spin-one 
pairing involves small gaps, as well as nodes in the gap 
parameter, and the emissivity is expected to be larger 
than that of matter with spin zero pairing. The emissiv- 
ity of two-flavor quark matter with (uu) and (dd) pairing 
was studied for different spin-one orde r parameters in 
( Schmitt et al ]. l2006HWang et al. iiooi). The result can 
be written as 



457 2 ^6 
bdv 



i + |5(A/r) 



(163) 



where the M-quark and d-quark gaps are assumed to be 
identical. All spin - one p hases analyzed in ( Schmitt et al\ . 
l2006l : IWang etal\ . \200^ and described by Eq. con- 
tain ungapped modes similar to the 2SC phase. There- 
fore, the emissivity at low temperatures is simply one 
third of that of unpaired quark matter. (In the case of 
color-spin locking, all excitations become gappe d if one 
takes i nto account no n zero q uark masses (.Aguilera et all 
120051 : ISchmitt et all l2006t ) and/or more complicated 
struc tures of the order parameter ( Marhauser et all . 
l2007l ).) The contribution to p63)) that arises from the 
paired quark s, is described by th e nontrivial function 
g{A/T); see (jSchmitt et al\ . I2006D for the explicit form 
and numerical evaluation of this function for arbitrary 
temperatures. In Table IIIIl we present the behavior of 
this function for temperatures much smaller than the 
gap, r ^ A, for various single-flavor spin-one color 
superconducting phases. Although this contribution is 
small compared to the contribution of the ungapped 
modes, we can use it to show the effect of different 
(anisotropic) gap structures on the parametric behav- 
ior of the neutrino emissivity. We see that, while fully 
gapped modes lead to an exponential suppression of the 
emissivity, nodes in the gap weaken this suppression to 
a power law. The power law depends on the behavior of 
the gap in the vicinity of the nodes. 
The specific heat can be written as 



We show the suppression function IC{A/T) for the spe- 
cific heat in Table IIIIl We see that an exponential sup- 
pression of the emissivity goes along with an exponential 
suppression of the specific heat. 



VIII. COLOR SUPERCONDUCTIVITY IN NEUTRON 
STARS 

Neutron stars are the densest material objects in the 
universe, with masses of order that of the sun (Mq) and 
radii of order ten km. Depending on their mass and on 
the stiffness of the equation of state of the material of 
which they are composed, their central density lies be- 
tween ^ 3 and ~ 12 tim es nuclear saturation density 
(no = 0.16 nucleons/fm^) i Lattimer and Prakashl . [2OOII 
I2007D . Neutron stars consist of an outer crust made of 
a rigid lattice of positive ions embedded within a fluid 
of electrons an d (in the inner layer of the c rust) super- 
fluid neutrons ( Negele and Vautherinl . [l973[ ). Inside this 
crust, one finds a fluid "mantle" consisting of neutrons 
and protons, both likely superfluid, and electrons. De- 
termining the composition of neutron star cores, namely 
of the densest matter in the universe, remains an out- 
standing challenge.^ If the nuclear equation of state is 
stiff enough, neutron stars are made of neutrons, protons 
and electrons all the way down to their centers. If higher 
densities are reached, other phases of baryonic matter 
(including either a pion condensate (iBahcall and WoH , 



1965 



1972 



Bavm and Campbelll . Il978l: iMigdall. Il97ll: ISawver 
ScalapinoT 1972), a kaon condensate ( Brown . 1995 



Kaplan and Nelsonl.ll986l). or a nonzer o density of one or 



several hyperons ( Glendennind . [l985l )) may result. Or, 
neutron star cores may be made of color superconducting 
quark matter. 

The density at which the transition from baryonic mat- 
ter to quark matter occurs is not known; this depends on 
a comparison between the equations of state for both, 
which is not well-determined for either. Very roughly, we 
expect this transition to occur when the density exceeds 
one nucleon per nucleon volume, a criterion which sug- 
gests a transition to quark matter at densities > 3no. 
The question we shall pose in this section is how as- 
trophysical observation of neutron stars could determine 
whether they do or do not contain quark matter within 
their cores. We have seen throughout the earlier sections 
of this review that quark matter at potentially accessi- 
ble densities may be in the CFL phase, with all quarks 
paired, or may be in one of a number of possible phases 
in which there are some unpaired quarks, some of which 
are spatially inhomogeneous. If quark matter does exist 
within neutron stars, with their temperatures far below 



For review articles on neutron stars as laboratories for under- 



cv = T{fj.l + fij) 



i + ^/C(A/T) 



(164) 
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phase 


gap structure 


oc 


ic{4>) oc 


CSL 


isotropic (no nodes) 


0exp(— '/20) 


0"/2exp(-^0) 


planar 


anisotropic (no nodes) 


0^/2 exp(-0) 


(^2exp(-0) 


polar 


point nodes (linear) 






A 


point nodes (quadratic) 







TABLE III Suppression function C/(0) for neutrino emissivity in direct Urea processes and suppression function /C(</)) for specific 
heat for four spin-one color-superconducting phases (abbreviating (j) = ^/T, and everything in the limit (j) — » oo). While fully 
gapped modes yield exponential suppression, nodes in the gap yield power law suppressions. The gap functions in the polar 
and A phases differ in the angular direction in the vicinity of the point nodes. A linear behavior leads to a stronger suppression 
than a quadratic behavior. 



the critical temperatures for these paired phases, it will 
be in some color superconducting phase. We shall see 
in this section that these different phases have different 
observational consequences, making it possible for a com- 
bination of different types of observational data to cast 
light upon the question of which phase of color super- 
conducting quark matter is favored in the QCD phase 
diagram, if in fact neutron stars do feature quark matter 
cores. 

Before turning to the signatures of quark matter in 
neutron star cores, we mention here the more radical 
possibility that nuclear matter in bulk is metastable at 
zero pressure, with the true ground state of strongly in- 
teracting matter in the infinite volume limit being color 
superconducting three- flavor quark matte r. According to 



this "strange qu ark n iatter hypothe sis" ( Bodmerl . Il97ll : 
iFarhi and Jaffd . 119841 : IWitteni . ll984) . ordinary nuclei are 
either stabilized by virtue of their small size or are 
metastable with lifetimes vastly exceeding the age of 
the universe. If this hypothesis is correct, some of the 
stars that we think are neutron stars m ay be strange 
stars , made entirely of quark matter (lAlcock et al. . 
19861: lAlcock and Olintol . Il988l: iFarhi and Jaffd . Il984 : 



Haensel et all Il986[) . Strange stars may have a thin 



crust (of order 100 meters thick) of positive ions sus- 
pend ed above the quark matter surface by an electric 
field (jAlcock et all . Il986l ). or they may have a compa- 
rably thin crust of positive ions embedded within the 
(neg atively charged) ou t er layer of the q uark matter it- 
self (|Alford et all l2006l : iJaikumar et~d{ . l2006al ). They 
cannot, however, have a conventional, km thick, crust. 
And, there are many indications that neutron stars in 
fact do have conventional crusts. For example, the rich 
phenomenology of X-ray bursts is well-understood only 
within this setting. More recent evidence comes from the 
analysis of the quasi-periodic oscillations with frequencies 
in the tens o f Hz detecte d in the aftermath of magnetar 
superbursts (llsrael et al . 20051: IStrohmavcr and Wattsl . 
[2005, 200l; IWatts and Strohmav cr. 2006, 2^007|), which 
can be understood as sei smic oscillations of a conven- 
tiona l neutron star crust (IStrohmaver and Wattsl . 120051 . 
l2006l : IWatts and Strohmaveii l2007l ) whereas the thin 
crusts of a strange star would oscillate at much higher 
frequencies ( Watts and Reddvl . [2007i ) . Even if most com- 
pact stars are neutron stars not strange stars, it remains a 



logical possibility that some strange stars exist, meaning 
that all ordinary neutron stars are metastable. Although 
possible this scenario is unlikely, given that merger events 
in which strange stars in an inspiralling binary are tidally 
disrupted would litter the universe with small chunks of 
quark matter ( "strangclcts" ) and one must then under- 
stand why these have not catalyzed the conversion of all 
neutron stars to strange stars ([Fri edman and Caldwcl|. 
Il99l[ ). We shall devote the remainder of this section to 
the more challenging task of using observational data to 
constrain the more conservative scenario that quark mat- 
ter exists only above some nonzero transition pressure, 
namely within the cores of conventional neutron stars. 



A. Mass-radius relation 

It has long been a central goal of neutron star astro- 
physics to measure the masses M and radii R of many 
neutron stars to a reasonable accuracy. Mapping out the 
curve in the mass-radius plane along which neutron stars 
are found would yield a strong constraint on the equa- 
tion of state of dense matter. As this program represents 
such a large fraction of the effort to use observations of 
neutron stars to constrain dense matter physics, we begin 
by considering its implications for the presence of quark 
matter within neutron star cores. 

The larger the maximum mass that can be attained 
by a neutron star, the stiffer the equation of state of 
dense matter, and if stars with masses close to 2 M© 
are found then the existence of phases with a soft equa- 
tion of state, such as baryonic matter with kaon or 
pion condensation, can be ruled out. However, although 
the quark matter equation of state is not known from 
first principles, it may easily be as stiff as the stiffer 
equations of state posited for ordinary nuclear mat- 
ter, and neutron stars with quark matter cores can in 
fact reach masses of order 2 (lAlford et all 12005. 



20011 : iRiister and Rischkd . |2004| ) 



Baldo et all l2003t iBlaschke et all l2007l : iFraga et all 



The equation of state for CFL quark matter can be 
parametrized to a good approximation as ( Alford et all 
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l2005ari 



n = -P = -- 



47^2 



(l-c)^4 + ^(A/2-4AV+Sc 



(165) 

If c were zero, the /i^ term would be that for nonin- 
ter acting quarks; c parametrizes the leading effect of 
interactions, modifying the relation between pp and 
/i. At high densities, c = 2 as /tt to leading order 
in the strong coupling constant (jBavm and Chinl . Il976t 
iFreedman and McLerranl . Il977[ ). Analysis of higher or- 
der co rrections suggests t hat c > 0.3 at accessible den- 
sities ( Fraga et all 200lt ). BeS can be thought of as 
parametrizing our ignorance of the fi at which the nuclear 
matter to quark matter transition occurs. The M'^jj? 
term is the leading effect of the strange quark mass, 
and is common to all quark matter phases. The pres- 
sure of a color superconducting phase with less pairing 
than in the CFL phase would have a smaller coefficient 
of the A^/x^ term, and would also differ at order M^, 
here lumped into a change in -Beff- Because pairing is a 
Fermi surface phenomenon, it only modifies the term, 
leaving the larger term untouched. However, it can 
nevertheless be important because at accessible densi- 
ties the term is largely can celled by i?cff, enha.ncing 
the importance of the ter m ( Afford and Reddvl . l2003t 
iLugones and Horvathl . l2002f ) . Remarkably, and perhaps 
coincidentally, if we make the (reasonable) parameter 
choices c = 0.3, = 275 MeV and A = 100 MeV 
and choose -Boff such that nuclear matter gives way to 
CFL quark matter at the relatively low density 1.5 ng, 
then over the entire range of higher densities relevant to 
neutron stars the quark matter equation of state (|165[) 
is almost indistinguishable from the nuclear equation 
of state due to Akmal, Pa ndharipande and Ravcnhall 
(APR) (jAkmal et all Il998f) that is one of the stiffest 
nuclear equations of state in the c ompendium found in 
(ILattimer and Prakashl[2003 . l2007^ . Neutron stars made 
entirely of nuclear matter with the APR equation of 
state and neutron stars with a quark matter core with 
the equation of state (|165p with the parameters just de- 
scribed fall along almost indistinguishable curves on a 
mass vs. radius plot, with the most significant difference 
being that the APR equation of state admits neutron 
stars with maximum mass 2.3M0, whereas the introduc- 
tion of a quark ma tter core reduces the maximum mass 
slightly, to 2.OM0 (|Alford et aZ.I . l2005a^ . 

The similarity between a representative quark matter 
equation of state and a representative nuclear equa- 
tion of state makes clear that it will be very hard to 
use a future determination of the equation of state 
to discern the presence of quark matter. However, 
although the numbers in the above paragraph should 
be taken as indicative rather than definitive, they do 
suggest that the existence of a neutron star whose mass 
was reliably determined to be > 2Mq would make it 
hard to envision such a star (and hence any lighter 
stars) having a quark matter core of any appreciable 
size. Now that the mass of PSR J0751+1807 has been 



revised downward from (2.1 zfc O.2)M0 llNice et al\ . 
I2OO5I) to (1.261:^2)^0 (|Nice and Stairl l2007^ ■ 
the heaviest known neutron star orbited by a 
white dwarf is PSR J0621-I-1002, whose mass is 



1.69t;^;ig)M0 (iNice and Staird . l2007l : ISplaver et ai\ . 



I2OO2D . Also, one of the two pulsars Ter5I and Ter5J (in 
a globular cluster) must have a mass that is > I .68M0 
at the 95% confidence level ( Ransom et al\ . l2005t ). 
and the mass of the X-ray pulsar Vela X-1 is above 
1.6 Mq (|Barziv et al\ . \2Qm . 

Given our lack of knowledge of the equations of state 
for nuclear and quark matter, measuring neutron star 
masses and radii alone do not allow us to reach our goals. 



B. Signatures of the compactness of neutron stars 

If we could detect gravity waves from neutron stars 
spiraling into black holes in binary systems, the gravi- 
tational wave form during the last few orbits, when the 
neutron star is being tidally disrupted, will encode infor- 
mation about the density profile of the neutron star. For 
example, upon assuming a conventional density profile, 
the gravi ty wave form encodes information about the ra- 
tio M/ R ((Faber et aLl . l200^ . essentially via encoding the 
value of the orbital frequency at which tidal deformation 
becomes significant. This suggests a scenario in which 
the presence of an interface separating a denser quark 
core from a less-dense nuclear mantle could manifest it- 
self via the existence of two orbital frequency scales in 
the wave form, the first being that at which the outer 
layers are deformed while the denser quark core remains 
spherical and the second being the time at wh ich even the 
quark core is disrupted (|Alford et aZ.I . l2001b[ ). This idea 
must be tested in numerical relativity calculations, and 
it may turn out to be better formulated in some other 
way. For example, perhaps the gravity wave form can be 
used to constrain the first few moments of the density 
profile, and this information can then be used to con- 
trast neutron stars with standard density profiles char- 
acterized by a single length scale R with those which are 
anomalously compact because they have a "step" in their 
density profile. Whatever the best formulation turns out 
to be, it seems clear that if LIGO sees events in which 
the tidal disruption of a neutron star occurs within the 
LIGO band-width, the gravity wave data will constrain 
the "compactness" of the neutron star, providing infor- 
mation about the density profile that is complementary 
to that obtained from a mass-radius relation. 

If there is a "step" in the density profile at an inter- 
face, LIGO gravity waves may provide evidence for its 
presence. But, should a density step be expected if color 
superconducting quark matter is found in the core of a 
neutron star? There are two qualitatively distinct possi- 
bilities for the density profile, depending on the surface 
tension of the quark matter/nuclear matter interface a. 
If a is large enough, there will be a stable, sharp, inter- 
face between two phases having different densities (but 



53 



the same chemical potential). If a is small enough, it 
becomes favorable instead to form a macroscopic vol- 
ume filled with a net-neutral mixture of droplets of neg- 
atively charged quark matter and positively charged nu- 
clear matter, see Sec. IIII.Hl which allows a continuous 
density profile. The distinction between these two sce- 
narios has been analyzed quantitatively for the case of a 
first order ph ase transition from n uclear matter to CFL 
quark matter ( Alford et aLl . [2001bl ). This is the simplest 
possible phase diagram of QCD, with a single transition 
between the phases known to exist at nuclear density and 
at asymptotically high density. We have seen earlier in 
this review that this simple QCD phase diagram is ob- 
tained if AcFL is large enough, allowing CFL pairing to 
fend off stresses that seek to split Fermi surfaces, all the 
way down in density until the nuclear matter takes over 
from quark matter. A sharp interface between the (elec- 
trically insulating) CFL phase and (electrically conduct- 
ing) nuclear matter features charged boundary layers on 
either side of the interface, which play an important role 
in determining the <j above which this st ep in the den- 
sity profile is stable (jAlford et al The critical 
a is about 40 MeV/fm'^, lower than dimensional analy- 
sis would indicate should be expected, meaning that the 
sharp interface with a density step is more likely than a 
mixture of charged components. The increase in the den- 
sity at the interface can easily be by a factor of two. The 
critical a above which a sharp interface is favored has not 
been evaluated for the case of a first order phase transi- 
tion between nuclear matter and color superconducting 
phases other than the CFL phase. 

It is also possible that the long term analysis of the 
binar y double pu lsar PSR J0737-3039A (jBurgav et all 
l2003l : iLvne et al\ . [2Qo3) may yield a measurement of 
the moment of inert ia of t his 1. 3 4 solar mass neutron 
2004 l2006t iLattimer and Schutzl . 
This could be another route 



star (IKramer et al 



20041 ) 



l2005l : lMorrison et al\ . 
to constraining the compactness of a neutron star, and 
perhaps gaining evidence for or against a step in the den- 
sity profile of this star. 



C. Cooling 

The avenues of investigation that we have described so 
far may constrain the possible existence of quark mat- 
ter within neutron star cores, but they are not sensitive 
to the differences among different color superconducting 
phases of quark matter. We turn now to the first of 
three observational signatures that have the potential to 
differentiate between CFL quark matter and other color 
superconducting phases. 

Within less than a minute of its birth in a supernova, 
a neutron star cools below about 1 MeV and becomes 
transparent to neutrinos. For the next million years or 
so it cools mainly via neutrino emission from its inte- 
rior. Photon emission from the surface becomes domi- 
nant only later than that. This means that information 



about properties of the interior, in particular its neutrino 
emissivity and heat capacity, can be inferred from mea- 
surements of the temperature and age of neutron stars. 
Because all forms of dense i natte r are good heat con- 
ductors (|Shovkovv and Ellisl . [2003 ). neutron star interi- 
ors are isothermal and the rate at which they cool is de- 
termined by the volume integrals over the entire interior 
of the local emissivity and the local specific heat. This 
means that the cooling tends to be dominated by the 
properties of whichever phase has the highest neutrino 
emissivity and whichever phase has the highest specific 
heat. 

Different forms of dense matter fall into three cate- 
gories, ordered by decreasing neutrino emissivity. The 
first category includes any phase of matter that can emit 
neutrinos via direct Urea processes, yielding an emissiv- 
ity cx T^. Examples include unpaired quark matter, 
phases of quark matter with some unpaired quarks in- 
cluding the crystalline phases and the phases with single 
flavor pairing in Table III, baryonic matter containing 
hyperons, nucleonic matter augmented by either a pion 
or a kaon condensate, and even ordinary nuclear matter 
at sufficiently high densities that the proton fraction ex- 
ceeds about 0.1. For the specific case of unpaired quark 
ma tter, the emissivity is given by (fleOl) (|lwamotal . Il980l 
which can be written as 



ei, ~ (4 X lO^^erg cm ^ s ^) x 



a. 



0.5 V 100 MeV / V500 MeV 



T 



109 K 



(166) 



where we have taken ~ / (4/^) , appropriate for neu- 
tral unpaired quark matter. (Note that as ~ 0.5 is com- 
parable to the value c ^ 0.3 that we used in Sec. IVIII.Al 
according to the lowest order relation c = 2as/Tr. The 
q;s/0.5 factor in (|166p could be replaced by c/0.3.) The 
emissivity of other phases of quark matter in which 
only some quarks are unpaired, including the crystalline 
phases, is reduced relative to (|166|) . but only by factors 
of order unity. 

Ordinary nuclear matter at densities not too far above 
no, where the proton fraction is less than 0.1, falls into 
a second category in which there is no phase space for 
direct Urea processes and neutrino emission occurs only 
via modified Urea processes like n + X^p + X + e + D 
with X some spectator nucleon, giving the much lower 
emissivity 



(1.2 X lO^Vg cm"'^ s"') — 



no 



2/3 



T 



109 K 



(167) 

Neutron stars whose interiors emit neutrinos at this 
rate, perhaps modified by effects of nucleon superfluid- 
ity. cool following a family of standard cooling curves (see 
(|Page et aZ.I . |2004| : lYakovle^ et al. I. I2OOII . I2OO5D and refer- 
ences therein), taking 10^ to 10^ years to cool below 10® 
K. 

CFL quark matter constitutes a third category. As 
we have seen in Sec. IVIIi it is unique among all phases 
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of dense matter in having an emissivity cx T^^ that is 
many orders of magnitude smaUcr than (|167p . Further- 
more, whereas all other phases of dense matter have a 
specific heat oc T, in the CFL phase the specific heat is 
controlled by bosonic excitations making it oc T^. This 
means that if a neutron star has a CFL core, the total 
neutrino emissivity and the total heat capacity of the star 
are both utterly dominated by the contributions of the 
outer layers, whether these are made of nuclear matter 
or of some phase that admits direct Urea reactions. The 
CFL core holds little heat, and emits few neutrinos, but 
is a good conductor and so stays at the same tempera- 
ture as the rest of the star. The rest of the star controls 
how the star cools. 

Finally, the single-flavor color superconducting phases 
are interesting because they represent a potential transi- 
tion from the first to the thir d categories ( Aguilera et all . 
I2OO5I : iGrigorian et all I2OO5I ): their critical temperatures 
are so low that if some quarks can only pair in spin-one 
channels, they will not pair until after the star has cooled 
through an initi al epoch of direct Urea emission; and , 
in ce r tain cases (lAguilera et all l2005l : iMarhauser et all 
I2OO7I : ISchmitt et all I2006D alfguarks can be gapped be- 
low the critical temperature for color-spin locked pairing, 
meaning that these phases ultimately become like CFL 
quark matter, playing no role in the cooling of the star 
which at late times will be controlled by the modified 
Urea processes in the nuclear matter mantle. 

We can now describe a possible future path to the dis- 
covery of CFL quark matter cores within neutron stars. 
Suppose that LIGO detects the gravitational waves from 
the tidal disruption of a neutron star with some known 
mass spiralling into a black hole and, as we discussed in 
Sec. IVIILBl suppose that evidence is found that the den- 
sity profile of the neutron star has a denser core within 
a less dense mantle, consistent with the existence of a 
step in the density profile. Suppose furthermore that 
it was understood by then that neutron stars with that 
mass cool following one of the family of standard cooling 
curves, meaning that there can be no component of their 
interior within which direct Urea processes are allowed at 
any time. This combination of observations would rule 
out the possibility that the dense core, inside the density 
step, contained any of the color superconducting phases 
that we have discussed except CFL. 

The scenario above may be unlikely, because there are 
a growing number of lines of evidence that although the 
cooling of many neutron stars is broadly consistent with 
the standard cooling curves, some fraction of neutron 
stars cool much more quickly. Examples of neutron stars 
that are too cold for their age includ e those in the su- 
pernova remnants 3C5 8 and CTAl (|Page et all |2004 
iPage and Reddvl . |2006[ ). A second, less direct, piece of 
evidence is provided by an unsuccessful search for the 
X-ray emission from a cooliiig neutron star in 15 other 
super nova remnants (|KaDlanl . l2007l : iKaplan et al\ . l2006l 
|2004| ). Ahhough some of these supernovae may have been 
Type I A supernovae which do not produce neutron stars, 



and although some may have produced black holes, it is 
likely that many of these supernovae remnants do contain 
neutron stars. Their nonobservation results in an upper 
limit on their temperatures, and in all cases this upper 
limit falls below the standard cooling curves. A third 
line of evidence comes from neutron stars that undergo 
transient bouts of accretion ( Brown et aZ.I . [l998[ ) . X-ray 
observations of one of these, SAX J1808. 4-3658, during 
its quiescent phase yield an u pper limit on the th ermal 
luminosity of the neutron star ( Heinke et aZ.I . [20071 ) . The 
mean accretion rate averaged over many transient ac- 
cretion episodes is known, meaning that the average ac- 
cretion heating of the star is known. The fact that the 
thermal luminosity is as low as it is means that the accre- 
tion heating of the star must be balanced by cooling by 
neutrino emission at a rate that far exceeds (|167p . The 
emissivity for unpaired quark matter (jl60[) is consistent 
with the data, as are the direct Urea rates for sufficiently 
dense nuclear matter and for hyperon matter. Pion con- 
densation or kaon condensation yield emissivities that are 
proportional to but with prefactors that are about two 
orders of magnitude smaller than that in (|160p . and are 
ruled out as explan ations for the ability of SAX J1808.4- 
3658 to keep cool (|Heinke et aZ.I . l2007i) . Similar conclu- 
sions can also be inferred from the (even lower) limit on 
the quiesc ent luminosity of the soft X-ray transient IH 
1905-hOOO (|Chakrabartvi . l2007l : lJonk~ et al. l,[2003,[2Q0i), 
although in this instance the time-averaged accretion rate 
is not as well known. 

By now it certainly seems clear that some neutron 
stars cool much faster than others. It is then reason- 
able to speculate that lighter neutron stars cool following 
the standard cooling curve and arc composed of nuclear 
matter throughout whereas, based on the three lines of 
evidence above, heavier neutron stars cool faster because 
they contain some form of dense matter that can radi- 
ate neutrinos via the direct Urea process. This could be 
quark matter in one of the non-CFL color superconduct- 
ing phases, but there are other, baryonic, possibilities. 
If this speculation is correct, then if neutron stars con- 
tain CFL cores they must be "inner cores" , within an 
outer core made of whatever is responsible for the rapid 
neutrino emission. 



D. r-modes limiting pulsar spins 

A rapidly spinning neutron star will quickly slow 
down if it is unstable with respect to bulk flows 
known as Rossby modes, or r-modes, whose restoring 
force is due to the Coriolis effect and which trans- 
fer the star's angular momentum into gravitatio n al ra- 
diation (lAnderssonl 19981 : Andersson and Comerl . 20011 : 



lAndersson erami999a llb[ |Friedman and Morsinkl . ll998[ ). 

For any given interior composition and temperature, 
above some critical spin frequency there is an insta- 
bility which leads to an exponentially growing r-mode. 
This means that as a neutron star is spun up by ac- 
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crction. its spin will be limited by a value very slightly 
above this critical frequency, at which the accretion 
torque is bala nced by gravitat ional radiati on from the 
r-mode flows (Andersson et ai . 1999a.b : .Bildstenl . Il998t 
ILindblom et aLl . ll998tlOwen et aLl . ll998[) . From a micro- 
physical point of view, the r-mode instability is limited 
by viscous damping: the greater the damping, the higher 
the critical spin above which r-modes become unstable. 
The critical frequency is controlled by the shear viscos- 
ity in some regimes of temperature (typically lower) and 
by the bulk viscosity in others (typically higher). This 
means that the existence of pulsars with a given spin, as 
well as any observational evidence for an upper limit on 
pulsar spins, can yield constraints on the viscosities of 
neutron star interiors. 

There is observational evidence for a physical limit on 
pulsar spins. The fastest known pulsar is a recently dis- 
cover ed radio pulsar spinning at 716 Hz ( Hessels et all 
hood ). However, it is not easy to draw inferences from 
the distribution of spins of the many known radio pul- 
sars as to whether 716 Hz is close to some physical limit 
on the spin frequency because there are significant obser- 
vational biases that make it harder to find faster radio 
pulsars. The most rapid pulsars are "recycled", meaning 
that they were spun-up during an episode of accretion 
from a binary companion. During such accretion, a neu- 
tron star may be visible as an X-ray pulsar. The spin 
frequencies of the 13 known millisecond X-ray pulsars 
lie between 270 and 619 Hz. What makes this signif- 
icant is first of all that the episodes of accretion have 
long enough durations that they could easily spin a neu- 
tron star up beyond 1000 Hz, and second of all that 
there are no selection biases that preclude the discov- 
ery of X-ray pulsars with frequencies as l arge as 2000 
Hz (|Chakrabartvl |2004 IChakrabartv et~ai\ . l2003f) . Anal- 
ysis of the observed distribution of X-ray pulsar spin 
frequencies leads to two conclusions: first, the distri- 
bution is consistent with being uniform;^ and, second, 
there is some physical effect that sets a limit on the al- 
lowed spin of a p ulsar which (with 95% confidence) is at 
730 H z or lower ( ChakrabartvL l2004l : IChakrabartv et all . 
|2003( ). It is unlikely that a spin-limit in this vicinity 



These data thus rule out a prop osal for how small quark matte r 
cores could have been detected l|Glendenning and Webeil . l200lh . 
If slowly-rotating neutron stars just barely reach quark-matter 
densities in their center, then rapidly spinning oblate neutron 
stars, which have slightly lower central density, will not con- 
tain quark matter. This "spinning out" of a quark matter core 
could be detected either by ano malies in braking indices o f radio 
pulsars that are slowing down l lGlendenning et all . Il997f) or by 
anomalous populatio n statistics of X-ray pulsars tha t are being 
spun up by accretion l lGlendennine and Weberl . l20pi | '). The data 
on X-ray pulsars show no sign of such an effect jchakrabartvl . 
|2004| : IChakrabartv et all |2003| ) indicating that, if quark matter 
is present, spinning the star and making it oblate does not get 
rid of it. If neutron stars do have quark matter cores, therefore, 
the quark matter must occupy a reasonable fraction of the star. 



can be attributed to centrifugal break-up of the spin- 
ning neutron star: unless neutron star radii are larger 
than anticipated, this "mass shedding limit" is signifi- 
cantly higher, above 1 kHz. On the other hand, if the 
observed limit on pulsar spin frequencies is attributed to 
the onset of the r-mode instability, the resulting con- 
straint on the viscosities of neutron star interiors is 
broadly consistent with the viscosities of nuclear matter, 
although this consistency is somewhat loose given the 
uncertainties in neutron star densities and in their tem- 



peratures while being spun u p ([Andersson et "^■ Il999bl : 
iLevin and Ushomirskvl . l200l[ ). 



The physics of the r-mode instability definitively rules 
out the possibility that accreting X-ray binary pul- 
sars are strange star s that are com posed of CFL quark 
matter throughout ( Madsenl . l200Cl( ). From the results 
of Sec. IVIIi we can conclude that CFL quark matter 
has negligible shear damping, and signif icantly smaller 
bulk viscosity than nu c lear m atter. (See ( Haensel et all 
l200nl . l20nial : ISawved . I1989D . fo r calculations of bulk 
viscosity in nuclear matt er, ( Haensel et all l2001bl : 



ILindblom and Owenl . |2002|) for baryonic matter con- 
tainin g hyperons, and ( Chatterjee and BandvopadhvavL 
I2OO7I ) for baryonic matter containing a kaon condensate.) 
A CFL strange star would therefore have a critical fre- 
quency at which the r-mode instability sets in measured 
in Hz or fractions of Hz, in gross disagreement with the 
data on spin frequencies of both X-ray and radio pulsars.^ 



It is a very interesting question, at present unresolved, 
whether the presence of a CFL quark matter core within 
an ordinary neutron star introduces unstable r-modes at 
low spin frequencies. If there is a density step at the 
nuclear/CFL interface, there may be oscillation modes 
localized near that interface. The question is whether 
there arc r-modes that are sufficiently well localized on 
the CFL side of the interface that they are undamped, 
or whether the tails of the mode wave functions that ex- 
tend into the nuclear matter side of the interface result 
in enough damping to prevent the modes from becom- 
ing unstable. Nobody has solved for the r-mode wave 
functions for a rotating star whose density profile has a 
step at an interface, with viscous dissipation occurring 
on one side of the interface only.^ If it were to turn out 
that a star with a CFL core is even close to as unstable 
with respect to r-modes as a star that is made entirely of 
CFL matter, the existence of pulsars spinning with hun- 
dreds of Hz frequencies would immediately rule out the 
possibility that these neutron stars have CFL cores. 



* Strange stars made of unpaired quark matter o r of 2 SC quark 
matter can be consistent with the data l lMadsenl . l2000l'l . 

^ Certain other oscillation modes ("f-modes" and "g- modes") of a 
nonrotating neutron star whose density profile includes a density 
step have been computed l lSotani et aLl . [2OO2I) . 
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E. Supernova neutrinos 

The only time when a neutron star emits enough neu- 
trinos to be detectable on earth as a neutrino source is 
during the first few seconds after the supernova explosion. 
The time-of-arrival distribution of supernova neutrinos 
could teach us about possible phase transitions to CFL 
quark matter (ICarter and Reddvl. 2000l: Jaikuinar et all . 
' " 120041 : ll 



l2002l : iKundu and Reddvl . l2004l : iReddv et al\ . l2003al lb^. 

All phases of quark matter and nuclear matter except 
CFL have short enough mean free paths that the neutri- 
nos detected from a supernova are emitted from a sur- 
face of last scattering called the neutrinosphere, inside of 
which they were diffusing. This surface of last scattering 
moves inward to higher densities during the first seconds 
after the supernova, as the protoneutron star cools. Sup- 
pose that a volume in the core of the protoneutron star 
has made a transition into the CFL phase, in which neu- 
trinos scatter only off Goldstone bosons which are less 
numerous (number density oc rather than oc /x^T for 
ungapped quark excitations). As this core cools, the neu- 
trino mean free path within a CFL core becomes longer 
than in any phase of matter in which there arc unpaired 
quarks (or nucleons) off which the neutrinos can scat- 
ter. The last supernova neutrinos to arrive could carry 
information about conditions when the neutrinosphere 
reaches the CFL core. Perhaps there may even be en- 
hanced neutrino luminosity at the end of an otherwise 
dropping timc-of-arrival distribution, as all those neutri- 
nos that were previously trapped within the transparent 



nos tnat were previously trapped witmn the transp arent 
core fly out unimpeded ( Carter and Reddvl . 1200(11 ). De- 
termining whether this proposed signature can arise re- 
quires implementing the transition to a CFL core, with 
its long neutrino mean free paths, within a full-fledged 
simulation of neutrino transport during a supernova. 



F. Rigid quark matter and pulsar glitches 

The existence of a rigid crystalline color superconduct- 
ing core within neutron stars may have a variety of ob- 
servable consequences. For example, if some agency (like 
magnetic fields not aligned with the rotation axis) could 
maintain the rigid core in a shape that has a nonzero 
quadrupole moment, gravity waves would be emitted. 
The LIGO non-detec tion of such gravity waves from 
nearby neutron stars ( Abbott et all . l2007l ) already lim- 
its the possibility that they have rigid cores that are de- 
formed to the maximum extent allowed by the shear mod- 
ulus (|168p . upon assuming a range of possible breaking 
strains, an d this constraint will tigh t en as LIGO contin- 
ues to run (jHaskell et am2007l: [Lh] . l2007t ). (The analo- 
gous constraint o n strange sta rs that are rigid throughout 
was obtained in ( Owenl . l2005l ).) Perhaps the most excit- 
ing implication of a rigid core, however, is the possibility 
that (some) pulsar "glitches" could originate deep within 
a neutron star, in its quark matter core. 

A spinning neutron star observed as a pulsar gradually 



spins down as it loses rotational energy to electromag- 
netic radiation. But, every once in a while the angular 
velocity at the crust of the star is observed to increase 
suddenly in a dra matic event called a glitch. The stan- 
dard e x planation (lAlDarl.ll977t[Alpar eit adll984al.ll996 . 
1984bl : lAnderson and Itohl . Il975l: "Epstein and Bavm 



1992: Jonesl. Il997l: iLink and Epsteinl . Il996l: iLink et 

1993 : iPines and Alpar . 19851 ) requires the presence of a 
superfluid in some region of the star which also features 
a rigid array of spatial inhomogeneities which can pin the 
vortices in the rotating superfluid. In the standard ex- 
planation of pulsar glitches, these conditions are met in 
the inner crust of a neutron star which features a neu- 
tron superfluid coexisting with a rigid array of positively 
charged nuclei that may serve as vortex pinning sites. 
We shall see below that a rigid core made of crystalline 
color superconducting quark matter also meets the basic 
requirements. 

The viability of the standard scenario for the origin of 
pulsar glitches in neutron star crusts has recently been 
questioned (|Link . .2007ai) . Explaining the issue requires 
understanding how the basic requirements come into play 
in the generation of a glitch. As a spinning pulsar slowly 
loses angular momentum over years, since the angular 
momentum of any superfluid component of the star is 
proportional to the density of vortices the vortices "want" 
to move apart. However, if within the inner crust the vor- 
tices are pinned to a rigid structure, these vortices do not 
move and after a time this superfluid component of the 
star is spinning faster than the rest of the star. When 
the "tension" built up in the array of pinned vortices 
reaches a critical value, there is a sudden "avalanche" in 
which vortices unpin, move outwards reducing the angu- 
lar momentum of the superfluid, and then re-pin. As this 
superfluid suddenly loses angular momentum, the rest of 
the star, including in particular the surface whose angu- 
lar velocity is observed, speeds up — a glitch. We see 
that this scenario requires superfluidity coexisting with 
a rigid structure to which vortices can pin that does not 
easily deform when vortices pinned to it are under ten- 
sion. In very recent work. Link has questioned whether 
this scenario is viable because once neutron vortices are 
moving through the inner crust, as must happen during 
a glitch, they are so resist ant to bending that they may 
never re-pm (|Linkl . l2007ah . Link concludes that we do 
not have an understanding of any dynamics that could 
lead to the re-pinning of moving vortices, and hence that 
we do not currently understand glitches as a crustal phe- 
nomenon. 

We have seen in Sec. IVI.BI that if neutron star cores 
arc made of quark matter but Acfl is not large enough 
for this quark matter to be in the CFL phase, then all of 
the quark matter core — and hence a significant fraction 
of the moment of inertia of the star — may be in one of 
the crystaUine phases described in Sec. IVI.BI By virtue 
of being simultaneously supcrfluids and rigid solids, the 
crystalline phases of quark matter provide all the neces- 
sary conditions to be the locus in which (some) pulsar 
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glitches originate. Their shear moduli (|144p . namely 
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with A the gap parameter in the crystalline phase as 
in Fig. [ini make this form of quark matter 20 to 
1000 times more rigi d than the crust of a n eutro n 
star (jMannarelli et al\ . l2007l : IStrohmaver tt aA Il99lt) . 
and hence more than rigid enough for glitches to origi- 
nate within them. The crystalline phases are at the same 
time superfluid, and it is reasonable to expect that the 
superfluid vortices that will result when a neutron star 
with such a core rotates will have lower free energy if 
they are centered along the intersections of the nodal 
planes of the underlying crystal structure, i.e. along 
lines along which the condensate already vanishes even 
in the absence of a rotational vortex. A crude estimate 
of the pinning force on vortices within crystalline color 
superco nducting quark matter indicates that it is suf- 
ficient ( Mannarelli et ad l2007l ). So, the basic require- 
ments for superfluid vortices pinning to a rigid structure 
are all present. The central questions that remain to 
be addressed are the explicit construction of vortices in 
the crystalline phase and the calculation of their pin- 
ning force, as well as the calculation of the timescale 
over which sudden changes in the angular momentum 
of the core are communicated to the (observed) surface, 
presumably either via the common electron fluid or via 
magnetic stresses. 



Much theoretical work remains before the hypoth- 
esis that pulsar glitches originate within a crystalline 
color superconducting neutron star core is developed fully 
enough to allow it to confront data on the magnitudes, 
relaxation timescales, and repeat rates that characterize 
glitches. Nevertheless, this hypothesis offers one imme- 
diate advantage over the conventional scenario that re- 
lied on vortex pinning in the neutron star crust. It is 
impossible for a neutron star anywhere within which ro- 
tationa l vort i ces are pinned to precess on ^ year time 
scales (|Linkl . l2006l l2007bl : ISedrakian eTd\ . Il999[ ). and 
yet ther e is now evidenc e that several pulsars are pre- 
cessing (IChukwude et al. 1, 120031; IShabanova et all l200ll: 
IStairs et all l2000l ). Since all neutron stars have crusts, 
the precession of any pulsar is inconsistent with the pin- 
ning of vortices within the crust, a requirement in the 
standard explanation of glitches. On the other hand, per- 
haps not all neutron stars have crystalline quark matter 
cores — for example, perhaps the lightest neutron stars 
have nuclear matter cores. Then, if vortices are never 
pinned in the crust but are pinned within a crystalline 
quark matter core, those neutron stars that do have a 
crystalline quark matter core can glitch but cannot pre- 
cess while those that don't can process but cannot glitch. 
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